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Recent Developments in Aircraft Control 


ROBERT C. SEAMANS, JR.,* FRANK A. BARNES,* THEODORE B. GARBER, } 
AND VINCENT W. HOWARD** 


Massachusetts Institute of Technology 


ABSTRACT 


Inhabited aircraft can be controlled directly by the pilot or 
indirectly through the use of instruments, electronics, and servo 
equipment. Uninhabited aircraft can be controlled by means 
of programmers, radio links, or homing systems. In all case 
human beings must select and supervise the mission of the 
aircraft and must make command decisions when new informa- 
tion becomes available. It is felt that human beings can exercise 
judgment more effectively when they are relieved of routine 
operations, especially actions involving muscular fatigue and 
manual dexterity. 

However, the present trend toward supersonic speeds is 
introducing more violent changes in the aerodynamic loading 
of aircraft, especially in maneuvering flight. As a result, large 
hinge moment variations are expected. In addition, the center 
of pressure shifts position with Mach Number, causing significant 
changes in static stability. Consequently, equipment to position 
the control surfaces and to stabilize the aircraft in the maneuver 
desired by the human operator appears desirable in both in 
habited and uninhabited aircraft 

In this paper, the mechanization of a pitch control system is 
discussed and its design is analyzed in nondimensional form 
In the analysis, consideration is given to changes in static 
stability, to the effect of fuselage bending, and to the location 
of gyro elements in the aircraft. In addition, the ability of the 
control system to minimize the effect of gust interference is 


discussed. 


INTRODUCTION 


1 poy PRINCIPAL FUNCTION OF COMMERCIAL AVIATION 
is the transportation of passengers and freight. 
Regardless of whether a particular operation is sched- 
uled or nonscheduled, painful losses in revenue result 
whenever flights must be delayed or canceled. The 
military services also utilize aircraft for transportation 
and, in addition, use aircraft for bombing, tactical 
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support, and interception of enemy aircraft. At a 
time when such great reliance is being placed on com- 
mercial and military aircraft, flight operation is rapidly 
becoming more difficult as a result of the increasing 
numbers of aircraft in use and the higher speeds of 
flight now obtained. Accordingly, much attention 
is being paid toward means of improving flight opera- 
tions. 

Many factors may interfere with the conduct of 
a particular commercial or military flight. Bad 
weather may cause the flight to be delayed, rerouted, 
or canceled. In addition, the aircraft may be difficult 
to fly, the efficiency of pilot and crew may be reduced, 
and the auxiliary equipment in the aircraft may be 
inaccurate and unreliable. An ideal aircraft systemtt 
would suppress all interfering effects and allow the 
flights to operate on schedule, but aircraft systems, 
unfortunately, are not ideal. However, recent de- 
velopments are yielding important improvements. 
For example, improved navigation systems can muini- 
mize the effect of bad weather during cruise, in the 
traffic pattern, and during the subsequent landing. 
In addition, better understanding of both atmospheric 
turbulence and aircraft structures may lead to control 
systems that alleviate the effects of horizontal and 
vertical gusts. Human error is closely coupled with 
the difficulty of flying an aircraft —flight crews burdened 
with unnecessary activities caused by aircraft de- 
ficiencies are likely to make costly mistakes. Conse 
quently, considerable effort is being expended on 
improvements to air-frame and cockpit design. How- 
ever, more fundamental changes will be required before 
completely reliable all-weather operation is achieved. 

The discussion in this paper concerns itself with those 
improvements that are being made to the control of the 
angular orientation of an aircraft. An aircraft is 
difficult to maintain on its proper heading when it is 


+? The complete aircraft system includes the aircraft itself, all 
auxiliary equipment in the aircraft, the flight crew, and all 
supporting aids located on the ground 
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improperly damped in any of its modes of oscillation. 
These modes include the long- and short-period longi- 
tudinal oscillations and the ‘Dutch roll’ lateral 
oscillation. Other handling problems may arise from 
changes in static stability, from a decrease in the 
effectiveness of the control surfaces, or from excessive 
changes of any other stability parameters. Moreover, 
hinge moment variations can place inordinate demands 
on the pilot. 

Recently developed integrating gyro units can be 
used to sense orientational deviations of the aircraft 
resulting from marginal stability, atmospheric turbu- 
lence, and other interfering effects. These gyro units 
are effective when the aircraft is in any attitude, and 
consequently they are ideally suited to aircraft that 
require complete maneuverability. The integrating 
gyro units were developed in the Instrumentation 
Laboratory at M.I.T. under the direction of Dr. C. S. 
Draper, the work being done under the auspices of the 
Armament Laboratory of the Wright Air Development 
Center (USAF). Although the authors have applied 
integrating gyro units to a number of aircraft guidance 
and control systems, they were not directly involved in 
the engineering development of these gyro units. 
Hence, this paper discusses only the application of the 
integrating gyro units to aircraft guidance and control 
systems. 

The control systems discussed in this paper are 
based on developments for the Armament Laboratory 
of the Wright Air Development Center (USAF) and 
the Bureau of Ordnance of the U.S. Navy. The 
analysis is nondimensionalized in order to broaden the 
application to a large variety of aircraft systems. 
Various types of aircraft control systems are discussed 
in Sections I and II. The response of a system for 
controlling pitch orientation is presented in Section 
III, and the effects of gusts and fuselage bending on 
this system are included in Sections IV and V. The 
details of the mathematical analysis of the aircraft 
and control systems are contained in the appendixes. 


(I) ArIRCRAFT GUIDANCE AND CONTROL SYSTEMS 


On the second day of July in 1914, Lawrence Sperry 
won a $10,000 safety prize when he utilized an 
automatic pilot to guide an aircraft by a reviewing 
stand in Paris, France. Since that date, autopilots 
have come into widespread use. In operation, the 
pilot places his aircraft on course at the desired speed 
and altitude and then engages the autopilot. Two- 
degree-of-freedom gyro units in the autopilot sense 
deviations of the aircraft from trim in pitch, yaw, 
and roll. Electrical signals proportional to these 
deviations are transmitted to servomechanisms that 
displace the control surfaces and return the aircraft 
to the desired orientation. Rate gyro units are some- 
times included to improve the damping characteristics 
of the aircraft, and a pressure-sensing element is used 
to indicate changes in static pressure, thus providing 
a means for altitude correction. In addition, aircraft 
have been flown automatically along ILS (Instrument 
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Landing System) beams by introducing the ILS signal 
into the autopilot. 

The purpose of the autopilot is to assist the pilot 
in following prescribed flight paths. Consequently, the 
combination of radio aids, pilot, and autopilot provide 
an aircraft with a navigation system. Other types of 
aircraft guidance and control systems are orientational] 
control systems with manual guidance, aircraft automatic 
tracking systems, and missile homing systems. De- 
scriptions of these systems, together with discussions 
of their applications, are included in this section. 
The section concludes with a discussion of the operation 
of a generalized guidance and control system. 


(1.1) Navigation System 


The input to a navigation system is the correct 
flight path. The system output is the aircraft position. 
For example, the correct flight path is the center of 
the ILS beam for a blind-landing navigation system. 
When the aircraft is not at the center of the beam, a 
command signal is provided by the ILS receiver and a 
coupling unit. 

Prior to entry into the landing beam, the aircraft 
must follow a traffic pattern such that a maximum 
number of aircraft are maintained in the vicinity of 
the terminal airport consistent with aircraft safety. 
Currently, a block system is used with aircraft stacked 
in 1,000-ft. layers. When one airplane lands, each 
aircraft is advised to descend to the next lower level. 
The correct flight altitude is designated by the control 
tower and the actual altitude is measured by a pressure 
altimeter. Any deviations are noted by the pilot, 
who repositions the aircraft accordingly. The limited 
capacity of data-processing equipment on the ground 
and of data relays from the ground to the aircraft has 
kept the traffic-control system from becoming auto- 
matic. However, systems now under development 
for the dispatch of interceptor aircraft can be adapted 
to airport traffic control at some future data. 

For navigating between airports, the flight plan 
establishes the correct flight path. Deviations from 
the flight plan are detected by the pilot and the navi- 
gator from visual observations of landmarks, from 
the use of radio aids, and from celestial observations. 
When deviations occur, the pilot alters the flight 
direction in order to return to the desired course. 
Fully automatic systems must provide a measure of 
latitude, longitude, and altitude, either by incorporat- 
ing radio links with the ground, by dead reckoning, or 
by use of inertial and celestial systems. 


(1.2) Orientational Control Systems with Manual 

Guidance 

When an aircraft is flying with the autopilot engaged, 
the pilot can change the angular orientation of the 
aircraft by disengaging the autopilot while the changes 
are being made, or, when small corrections are required, 
he can make corrections by use of turn and trim adjust- 
ments. Consequently, an autopilot can be used to 
provide limited orientational control but is unsuited 
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to the maneuverability requirements of certain military 
aircralt. 

The input to an orientational control system with 
manual guidance is the orientation established by the 
pilot. The system output is the actual orientation. 
For example, a pilot may decide to head toward a 
ground target that he has identified. To accomplish 
this, he must introduce corrections to the pitch and 
yaw control systems until the correct orientation is 
reached. The pilot introduces the commands to the 
control system by means of a control stick. The 
function of the control system is to maintain an angular 
velocity proportional to the command and to minimize 
all interfering effects caused by atmospheric turbulence. 


(1.3) Aircraft Automatic Tracking Systems 


The orientational control system with manual 
guidance can be used in conjunction with sighting 
equipment to aim the aircraft for gunnery and rocketry 
under visual conditions. However, radar is required for 
target detections in all-weather operations, and, since 
automatic tracking is desirable in fast-closing attacks, 
it appears that aircraft automatic tracking systems 
will be used whenever possible in all-weather inter- 
ceptors. 

The input to an aircraft automatic tracking system 
is the line from interceptor to target, commonly called 
the line of sight. The system output is the aircraft 
orientation. The correct orientation is offset from 
the line of sight to account for the ballistics of the 
projectile and the motion of the target during the time 
of flight of the projectile. Consequently, the correct 
orientation depends upon the type of armament used 
for interception and must be computed continuously 
as the target aspect varies during the firing run. 


(1.4) Missile Homing Systems 


In general, missile guidance can be divided into three 
phases: (1) launching, (2) mid-course, and (3) terminal 
guidance. Terminal guidance can be effected by 
transmitting commands from the ground or from 
inhabited aircraft, in which cases the target position 
is determined by equipment external to the missile. 
Terminal guidance can also be effected with a missile 
homing system, which requires a target detector in 
the missile. A missile homing system requires fewer 
data links than other terminal guidance systems but 
at a cost of greater missile complexity. 

The purpose of the missile homing system is to 
maneuver the missile into collision with the target. 
The input to the system is the target flight direction. 
The system output is the missile flight direction. The 
missile flight direction differs from the missile orienta- 
tion by the amount of the angle of attack and the angle 
of sideslip. When the target takes evasive action, its 
flight direction changes, and consequently the missile 
flight direction must be altered continuously. The 
difference between the correct flight direction and the 
actual flight direction must be determined, and a 
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TABLE 1] 
Aircraft Guidance and Control Systems 


Type of Input to the Output of the Output of 
Guidance Guidance and Guidance and the Control 
and Control Control Svstem, Control System, System, 

Sy stem Q(gex)yin q t q ut 
Navigation Correct flight Aircraft posi Aircraft flight 
system path tion direction 
Orientational Angular Aircraft angu Aircraft pitch 


orientation 
established tion 


lar orienta ind yaw 


control sys 
ingular 


tem with 


manual by the pilot velocities 
guidance 
Aircraft Line of sight Aircraft Aircraft pitch 


and yaw 
angular 
velocities 


angular 
orientation 


automatic to target 


tracking 


system 

Missile Target flight Missile flight Missile 
homing direction direction normal and 
system lateral ac- 


celerations 


command signal must be introduced into the control 


system. 


(1.5) Operation of a Generalized Guidance and Control 

System 

The four guidance and control systems discussed 
in this section are listed in Table | along with input, 
Giges)iny to each of the systems, and the output, Qies)outs 
of each of the systems. Each of the systems can be 
represented by the functional diagram of Fig. 1. It 
can be seen from Fig. 1 that the signal modifier supplies 
the control system with a command signal. Then the 
control system maneuvers the aircraft in order to 
change the guidance and control system output until 
the command signal is minimized. The command 
signal reaches a minimum value when the output from 
the sensing element matches the input to the guidance 
and control system. The relationship between system 
performance and the performance of its components is 
indicated in Fig. 1. The mathematical form used 
to express this relationship is called the relating 
function. The relating function of a component or 
system is the ratio of the Laplace transformation of 
the output to the Laplace transformation of the input, 
assuming zero initial conditions. For a more complete 
definition of the relating function, refer to Appendix B 
of this paper and to Chapter 17 of reference 1. 

Three guidance and control system relation functions 


are defined in Fig. 1: 


(1) The input-output relating function, 
(RF) (est ( 

d 'Wiges)in; *(ges)out 
(2) The interference-output relating function, 


(RF) (es 


3) The open-loop relating function, (RF) wesyor- 


Vint; “(ges)out 


The first two functions relate the output quantity to 
the input quantity and the interfering quantity re- 
spectively in closed-loop operation. The third function 
relates the sensing element output to the input quantity 
when the sensing element output is not fed back into 
the system. This open-loop relating function equals 
“the product of the relating functions of the signal 
modifier, the control system, and the sensing element. 
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(es)e x MU gcx) dt (t. 
Interfering Quantity Q(es)out 1d ive out | 1) 
= a ie ; 
he control system output quantities for the four 
types of guidance and control systems considered are 
uidan . : _ ~ . 
Guidonce ond Guidence ond | ong —_ oe and listed in Table 1. For example, the outputs from the 
Control Syst Control System Command t ontro ontrol system ‘ ° F 
restr Aesth Correction seen Control Syston System | Output Quantity control system to be used in the aircraft automatic | 
UgesXin) ~—“S” [(c)a]ycs | Signal 4c gcsKour) tracking system are the aircraft pitch and yaw angular j 
[Modifier ean Fe m : » ~ ¢ 
velocities. The description of a pitch rate contro] 
Guidance ° ° " ‘ 
a ae me system follows in Section II. { 
and Control System Control 
Sensing Element System 
Sensi ‘ . 
eee Element (Il) DESCRIPTION OF A PITCH RATE CONTROL System ' 
(RF). (RF), 4;- . The purpose of a pitch rate control system is to 
(RF), 7 sm (es)[ 465 4ges)out! ae Ee aie i : ! 
BCS)[Qees)in;Uges)out! = T+ Picea achieve in a satisfactorily short time a pitch angular ‘ 
(gces)OL : ° . . ° ° 
° velocity of an aircraft in proportion to an applied ] 
RF an : 
OM geaie_; « ( es)[ajnti gcs)out! input command. The pitch rate control system must 
int’ “(gcs)out! ~ x — sl ? 
aan 1+ (RF) eesOL also keep to a minimum the angular deviations in , 
where pitch which result from gust velocities. The control ‘ 
(RF)ges)OL * RP) sm'FF)(cs)[ 4.5% ges)out! (PP). system described in this section utilizes an integrating ¢ 
: . , ‘ ; ; yyro unit. ' 
Fic. 1. Functional diagram of aircraft guidance and control Byr ; : : : I 
system. The integrating gyro unit contains a gyro wheel ‘ 
that spins in synchronism with the gyro wheel supply P 
When a constant input is applied to the guidance voltage. The gyro wheel is supported by bearings ¢ 
and control system, it is desirable for the system attached to a gimbal frame pivoted with respect to 
output to change at a rate proportional to the correc- the gyro case at right angles to the gyro wheel spin f 
tion signal. In this manner, the output varies until axis. The axis about which the gimbal frame is : 
the correction is zero and thus reaches a steady-state pivoted is called the gyro output axis. In a properly 
condition without error in the system. However, for designed control system, the gimbal frame is never 
the output to change at a rate dependent upon the required to rotate more than two degrees with respect 1 
correction signal, the output must be proportional to the gyro case, and consequently the gyroscopic b 
to the integral of the correction. The required in- element produces a moment nearly equal to the product ‘ 
tegration takes place in the control system when the of the gyro wheel angular momentum, //, and the case c 
rate of change of the guidance and control system angular velocity about the input axis, which is per- s 
output, Giycs)outy IS proportional to the command signal, pendicular to the plane containing the gyro output p 
g-- In this paper, the control system output is by axis and the average position of the spin axis. The g 
definition equal or proportional to the rate of change of integrating gyro unit of Fig. 2 must be oriented with n 
the complete system output, i.e., the input axis along the J) axis of the aircraft in order fi 
g 
Gust Velocities Causing Interfering Moments and Forces on Aircraft t] 
d 
E 
Integrating Gyro Unit 
| Angular ~ | : I 
M Signal 
| loment Displacement | Modifier Ve 
Torque Applied of Torque 0 . Rate an 
Command ponents by Torque Summing | boa utput Elevator LI of T 
Voltage Torque | Current Torque | Generator | _/ rave Member Signal Signal | Voltage } Elevator [Displacement Aircraft [Pitch 
G #1 Generator u s 9g Generator >] Modifier Servo ; ircrart ° 
e& Amplifier itg qT if) Member A, | \ &sm 6 Di 
| | \ ol 
Integrating t 
| Moment Moment | Gyro ( 
Applied by| [Applied by | Output Voltage al 
| Gyroscopic] | Viscous - 
Element Damper | ’ ge 
| ge Mya | ol 
L ‘ 
1 Gyroscopic Viscous a | 
| Element —-——| Damper a] 
| : | di 
| al 
en sine te: St gin ae el 
su 
Fic. 2. Functional diagram of pitch rate control system. 
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to provide a moment proportional to the pitch angular 
velocity, 9. 

The torque summing member of Fig. 2 includes all 
the inertia elements that rotate with the gimbal frame 
and gyroscopic element. In addition to the moment 
pruduced by the gyroscopic element, two other mo- 
ments act on the torque summing member: 

(1) The torque summing member is restrained by 
the shearing action of a highly viscous fluid. The 


restraining moment applied by the fluid equals —c,Ag, 


where c, is the damping coefficient and A, is the rate 
at which the gimbal frame rotates with respect to the 
gyro case. This damping action is represented by the 
viscous damper block in the functional diagram of 
Fig. 2. 

(2) A moment is applied to the torque summing 
member by an_ electromagnetic four-pole torque 
generator. The generator stator is fixed to the gyro 
case, and the rotor is attached to the gyro gimbal 
frame. When current is supplied to the stator wind- 
ings, a moment is developed that is equal to the product 
of the torque generator sensitivity, S,,, and the torque 
generator current, 7,,. 

The moment of inertia of the torque summing 
member is kept small by design and consequently, to 


a close approximation, 
GiAg = Sigttg + HO (2-1) 


The torque generator current, 7,,, in Eq. (2-1) is supplied 
by an amplifier that maintains the current equal to 
Somp€e, the product of the amplifier sensitivity and the 
command voltage. The output from the integrating 
gyro unit is an a.c. voltage signal with an amplitude 
proportional to the angular displacement, A,. The 
gyro output voltage, e,, is provided by an electro- 
magnetic signal generator that has a four-pole stator 
fixed to the gyro case and a rotor attached to the gyro 
gimbal. The signal generator sensitivity, S,,, equals 
the ratio of the gyro output voltage to the angular 
displacement of the torque summing member, and hence 
Eq. (2-1) can be written 


I . ; Samp te ,) ? 
= (C)Sef[ (Gt )ect ala 22 


It can be seen from Eq. (2-2) that the gyro output 
voltage depends upon the integral of the pitch rate. 
This accounts for the name ‘‘integrating gyro.” 

The operation of the integrating gyro unit in the 
pitch rate control system can be analyzed by the use 
of Eq. (2-2). When a command voltage is applied 
to the control system, this voltage is first amplified 
and a proportional current is supplied to the torque 
generator of the integrating gyro unit. The gyro 
output voltage produced is then transmitted to the 
signal modifier, and the resulting voltage signal is 
applied to the elevator servo. Finally, the elevator 
displacement obtained changes the pitch rate of the 
aircraft until the moment produced by the gyroscopic 
element balances the moment applied to the torque 
summing member by the torque generator. When 
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this steady-state condition is reached, e, = 0, so that 


6 = [—(SumpS)/ He (2-2) 


Consequently, the rate of pitch is proportional to the 
“command voltage. In this steady-state condition, 
when gusts suddenly deviate the aircraft an amount 
Aé, the gyro output voltage changes an amount le, 
which can be determined from the expression 


+ 


Ae, = (£1/C,) SA? (2- 


The elevator is then displaced proportional to the 
change in gyro output voltage, thereby applying a 
moment to the aircraft which reduces the deviation 
in pitch. 

In order to obtain a large elevator displacement 
for a given pitch deviation, it is necessary—especially 
with marginally damped aircraft—to provide a signal 
proportional to pitch rate. This signal is developed 
in the signal modifier by electronic differentiation and 
summation, since from Eq. (2-2) 

6 = [—(SampSy)/Hee + [C,/(HS) lé, (2-5) 
Thus, the rate of pitch equals the sum of two terms, one 
proportional to the command voltage and the other 
proportional to the rate of change of the gyro voltage 
signal. Consequently, by differentiating the gyro 
output voltage in the signal modifier and by combining 
the differentiated signal with the command voltage, a 
signal proportional to pitch rate can be obtained. 

The relationship between the components of a pitch 
rate control system is analyzed mathematically in 
Appendix D, and the nondimensional response of a 
pitch rate control system is presented in Section (IIT), 


(III) RESPONSE OF THE PITCH RATE CONTROL SYSTEM 
TO COMMAND SIGNALS 


The selection of feedback sensitivities to obtain 
a suitable response of the pitch rate control system of 
Fig. 2 is considered next, assuming initially a rigid 
air frame and gust-free air. Later sections extend 
the investigation to include the effects of interference 
from gusts, air-frame elasticity, and location of the 
gyro unit. For the purposes of mathematical analysis, 
the system is represented as shown in Fig. 3, in which 
the parameters and paths of signal flow have been 
consolidated. The mathematical relationships in- 
dicated in Fig. 3 are derived in Appendix D. 


(3.1) Uncontrolled Aircraft Performance 


The aircraft longitudinal motion is described by 
three linearized force and moment equations defined 
and presented in Appendix A in both dimensional and 
nondimensional form, where the latter form employs 
the nondimensional time variable, Vt/uc, and the 
dimensionless aircraft stability derivatives defined by 
Eqs. (A-10)—(A-19). Numerical values of the stability 
coefficients are listed in Table 3 of Appendix A and 
are typical of values encountered with modern high- 
performance aircraft in horizontal wings-level flight. 
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Fic. 3. Mathematical representation of pitch rate control system. 
The rigid-aircraft relating functions defined and pressure. The uncontrolled aircraft is indicated as 


presented in Appendix B are functional relationships 
derived from the aerodynamic equations of Appendix 
A which relate output quantities (pitch angle, 6, and 
climb angle, y) to the input quantities (elevator de- 
flection, 6, and the ratio of vertical gust velocity to 
average air speed, w,/V). 

The rigid-aircraft characteristic equation presented 
in Appendix B takes the form 


A} + a3A° + dod" + aA + ado = 0 


and factors into the so-called ‘‘short-period”’ and “‘long- 
period” pairs of roots. Fig. 4 shows loci in the complex 
\-plane of these roots as a function of the static stability 
coefficient, m,. As defined in Fig. 4, the coefficient 
Ma iS proportional to the static margin—i.e., the 


distance between the aircraft centers of gravity and 


having static stability, neutral stability, or instability, 
according to whether m, is respectively negative, 
zero, or positive—i.e., according to whether the a 
pitching moment respectively tends to restore the 
aircraft to trim, is zero, or tends to displace the aircraft 
further from trim. It can be seer from the variation 
in the location of the roots that the influence of m, 
character of the aircraft’s performance is 
appreciable. With static stability, the aircraft has 
both short-period roots in the left-half \-plane, and, 
if m, < —0.15, the short-period pair is conjugate 
complex, characterizing damped sinusoidal motion. 
In the case of static instability (m, > 0), one short- 
period root is in the right-half plane, characterizing 
The pair of 


on the 


motion that is exponentially divergent. 
“long-period” roots (characterizing the so-called phuw- 
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goid motion) remains in the left-half plane and in the 
vicinity of X = O for all values of m,. 

The range —1 < m, < +1 is considered reasonable 
for the investigations of this paper, because changes in 
m, of unity or greater are not uncommon, especially 
for supersonic aircraft. For example, in passing from 
subsonic to supersonic flow, the center of pressure shifts 
aft one quarter of a chord length, thereby changing the 
value of m, by —1.2, assuming (OC,/Oa) + Cpo= 
48. Moreover, in the transonic region, the center of 
pressure also shifts forward in some aircraft, further 
aggravating this problem. In addition, the total 
variation in m, is increased an additional amount by 
any forward shift of the center of gravity such as occurs 
with the burning of fuel carried aft for the purpose of 
boosting the aircraft through the transonic region. 
Of course, an attempt should be made to keep in- 
habited aircraft statically stable as a matter of safety. 

Fig. 5 shows the frequency behavior—i.e., the 
sinusoidal steady-state behavior—for the uncontrolled 
aircraft with m, = —1, 0, and +1. This behavior 
is calculated from . 


((RF) 4189 ln=ja = (AR) apse? ? 4%] (3-2) 
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for various values of 2, where 2 is the nondimensional 
forcing frequency [(uc)/V]w of the assumed sinusoidal 
elevator motion. The ordinate of these curves is the 
phase angle, (PA) 415.9), of the sinusoidal pitch angle, 
6, relative to the applied sinusoidal elevator deflection, 
6. The abscissa gives the amplitude ratio (AR) 4)5.4 
and is equal to the ratio of the maximum amplitudes of 
the sinusoidal quantities, @ and 6. For comparison, 
the frequency performance is also plotted with the 
phugoid motion neglected—1.e., with 2, = x, = X%q = 
xX» = 0. Fig. 5 shows that for the aircraft considered: 

(1) The phugoid motion is not excited appreciably 
unless Q is less than about 0.5. 

(2) The value of m, in the range —1 < m, < +1 
affects only its low-frequency behavior (2 < 10), 
and this effect is primarily in (PA) 4)5.; with as much 
as 180° difference apparent. 


(3.2) Stability of the Closed-Loop System 


The response of the pitch rate control system of 
Fig. 3 to command inputs is derived in Appendix D 
and is given by Eqs. (D-20) and (D-21). Reference to 
these equations shows that the roots of the system 
characteristic equation are located in the complex 
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_Fic. 4. Loci of aircraft long- and short-period roots as a function of the static stability coefficient. In this figure the nondimen- 
sional roots are determined from the characteristic equation of Appendix B 


ri a a:\3 + 2d? + ar +a = 0 


This equation includes long- and short-period dynamics but not the effect of fuselage bending. 
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\-plane at those values of \ which satisfy 
1+ (RF) ¢s01 = 0 (3-3) 


where the open-loop relating function (RF) (so: of 
this control system is given by Eq. (D-22). For this 
system to be stable, of course, all roots of Eq. (3-3) 
must be in the left-half A-plane. The task of selecting 
suitable parameter values for this system is aided by 
(1) Nyquist’s criterion for stability? to guarantee 
that all roots of Eq. (3-3) do lie in the left-half \-plane; 
(2) the widely used methods of frequency analysis 
developed by Hall,* Bode,’ Nichols,* etc., to guarantee 
in addition a reasonable factor of safety against 
instability difficulties which is consistent with the 
accuracy to which the system’s physical behavior is 
known mathematically; and (3) analog computation 
that permits quantitative investigation of the system’s 
response to command inputs. 

For these purposes loci of frequency performance 
of the open-loop have been plotted in Fig. 6 (using 
nominal signal modifier and elevator servo charac- 
teristics) for each of four values of the sensitivity ratio 
Sies){o";51 / Ses) {0:81 The amplitude ratio has been 
normalized so that different values of the pitch angle 
feedback sensitivity, S(cs)(9.s; (with the preceding sensi- 
tivity ratio held constant) can be considered merely 
by shifting the scale of the abscissa. 

It can be seen in Fig. 6 that for certain critical values 
Sieoiesjer. Of the pitch angle feedback sensitivity, 

{(AR)csyor = 1 _ 

> with Q = OQ... 

\ (PA) cessor = —180° 
which is equivalent to the condition of Eq. (3-3) for 
A= +jQ.,. Eq. (3-4) shows therefore that with such 
critical adjustment the characteristic equation pos- 
sesses a pair of roots on the imaginary axis, indicating 
the system is capable of self-maintained oscillations 
at the frequency ®,.,... Such critical sensitivity values, 
which correspond to the intersection of a locus in 
Fig. 6 with the —1S80° ordinant, are listed in Table 2. 
Nyquist’s criterion indicates that these critical values 
are the maximum and minimum positive values of 
this sensitivity for stable operation. It is apparent, 
however, that no minimum positive value exists, if 


(3-4) 


the air frame is statically stable. 

The system’s margin of stability in terms of its 
“phase margin”’ and ‘‘gain margin’’ can be determined 
from Fig. 6 for any of the combination of stable feed- 
back sensitivities listed in Table 2. Phase margin is 
defined as 180° + (PA) ;-.0, determined for the 
frequency for which (AR)(.o. = 1. Phase margins 
of the order of 45° are often considered desirable. 
Gain margin specifies a relationship between the feed- 
back sensitivity used and the maximum or minimum 
values for which the system is stable. Eqs. (3-5) and 
(3-6) define two gain-margin relationships used in the 
discussions of this paper. 

Gain margin with respect to maximum sensitivity 


Sled [@:6]max. ~ S cs) [0:5] 9 - 


in per cent = 100~ (3-5) 


cs) [0:6] 
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Gain margin with respect to minimum sensitivity 
— § 


m (es) [4:6] cs) (6;5)min . , 
in per cent = 100 : (5-6) 

cs) [0:6] min 
Stable values of the pitch angle feedback sensitivity 
permitting operation with gain margins in excess of 
150 per cent and with various phase margins have 


also been determined from Fig. 6 and listed in Table 2. 


(3.3) Response to Step Commands 


Since more than one combination of the feedback 
sensitivities in Table 2 permits usually acceptable 
margins of stability, the final selection can be made 
from these on the basis of response to command signals. 
Fig. 7 shows the response to step command signals, 
determined by analog computer for each of the stable 
combinations of feedback sensitivity given in Table 2. 
These families of response curves were obtained with 
—1 and other system parameters as shown in 

It can be seen that there is considerable choice 
i.e., the time for @’ to reach and 
Further- 


m, = 
Fig. 6. 
in “‘response time’’ 
stay within 5 per cent of the final value. 
more, the character of the controlled aircraft's response 
to command signals is shown to be practically inde- 
pendent of the static stability coefficient m, over the 
range —1 < m, < +1 by the analog solutions in 
Fig. 8, even though a considerable variation in elevator 
response is required for differing amounts of static 
stability. 

Cases may arise, however, as in Section V, where 
uncertainties in the component's or aircraft’s charac- 
teristics warrant margins of stability somewhat larger 
than usual. Reference to Fig. 6 shows that the phase 
margin and the gain margin computed by Eq. (3-5) 
can all be increased by the reduction of S(..){¢.3; from 
the values selected in Table 2. Of course, if the 
aircraft has static instability then the minimum value, 
Sics){g.slmin. Cannot be approached too closely without 
producing a tendency to low-frequency instability. 
Fig. 9 illustrates the effect of reducing this sensitivity 
for an aircraft having static instability (m, = +1) 
and case ‘‘a’’ in this figure shows the poorly damped 
behavior that goes with a gain margin of only 115 
per cent with respect to S(.5){g.3)min. aS Computed from 
Eq. (3-6). The poorly damped behavior in the latter 
case indicates the system's tendency to a low-frequency 
instability, which instability would occur, of course, 


TABLE 2 
Stable and Critical Pitch Angle Feedback Sensitivities 
Determined from Fig. 6 


Values Selected for 
Operation with Margin 
of Stability 
(Gain Margins All 


—Critical Values Greater Than 150%) 


Values Of — $(¢5){0;3) max. -S(cs)[0;6] min Phase 

S(cs) (661 (forma = 1, (for ma = margin, 

Sics) [0:5] 0, —1) +1 only S(cs) {0;6] deg. 
0.4 9.4 0.37 3.5 45 
0.2 18 0.52 7.0 45 
0.1 36 0.70 12.0 45 
0.05 66 0 85 22.5 35 
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Fic. 6. Open-loop frequency performance of the pitch rate control system assuming a rigid aircraft. 
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if the minimum sensitivity value was actually to be 
used. An obvious remedy for this type of instability 
difficulty, if a larger sensitivity value cannot be em- 
ployed, is to decrease m, by moving the aircraft center 
of gravity forward, thereby reducing the minimum 
stable value of S(..)j9.4; and hence increasing the 
margin of stability 

It can be shown by reference to Fig. 6 that the 
magnitudes of the time lags associated with the elevator 
servo and signal modifier exert a primary influence on 
the value of S(.s)(g.5)max. and hence on the gain and phase 
margins that exist with respect to the latter sensitivity. 
It follows, therefore, that where control of aircraft 
having static instability is required, the component 
time lags may be too large to permit a sufficient amount 
of feedback sensitivity for successful control. In a 
good many applications of this type, it has been 
necessary to use hydraulic equipment in order to 
realize sufficiently small time lags. 


(IV) THe Errects or Gust INTERFERENCE ON THE 
P1ItcH RATE CONTROL SYSTEM 


In the analysis of the pitch rate control system 
discussed in Section III, it is assumed that the air 
mass is free from turbulence. Such atmospheric 
conditions seldom prevail, and, as a consequence, it 
is necessary to consider the effect of gusts upon the 
control system. The longitudinal equations of motion 
listed in Appendix A indicate the manner in which 
gusts exert forces and moments upon the air frame. 
The rigid-aircraft relating functions for a vertical gust 
velocity input are given in Appendix B for pitch and 
climb angle outputs. 

Atmospheric turbulence is a random phenomenon, 
and therefore statistical methods are usually utilized 
in the analysis of gust effects. The experimental work 
of Clementson® has been used to a great extent in 
specifying the gust characteristics discussed in this 
section. Data from reference 6 indicate that gusts 
are generated by a stationary random process and that 
the power spectral density of vertical gust velocities 
may be represented by the relationship 


(PSD) we V)} = 
\2re[(MS)we/V]}/{1 + [((Vre/uc)Q2]?} (4-1) 


where [(.17S)we/ |] is the mean square of the non- 
dimensional vertical gust velocity and rq is a time 
constant associated with the gust. Clementson found 
that the mean square of the vertical gust velocities 
is independent of the aircraft's forward speed. How- 
ever, the gust time constant, 7g, is inversely propor- 
tional to the forward speed, V. For averaged mete- 
orological conditions and for the conditions of this 


paper a nominal value of 7¢ is 2 sec. and 


[((S)we/V] = 46 XK 10-6 (4-2) 
The experimental results of reference 6 are found to 


be consistent with theoretical analyses of atmospheric 
turbulence. However, a complete investigation of the 


ce 
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effects of atmospheric turbulence would include gust 
disturbances along the longitudinal and pitch axes.” § 
The general problem is complicated by the fact that 
the various components of the gust velocity are statis- 
tically related, and cross effects must be considered. 

Once the characteristics of gusts have been estab- 
lished, their effect upon the aircraft control system 
dynamics may be determined by evaluating an integral 
of the following type. 


’ 1 f° V[(PSD)we;/V) 
((MS)gou] = | 
T Jo uC 


(AR)*(es)(wg/Vtoul@2 (4-38) 


The integrand of Eq. (4-3) is the power spectral density 
of the particular system output being considered. 


Thus 


[(PSD)qout] = [(PSD)we/ V(AR)*Ces)tw,/v (4-4) 


Mout 


The gust power spectral density is given by Eq. (4-1). 
The form of the control system amplitude ratio can 
be determined from the control system relating func- 
tions that are listed in Appendix D for various inputs 
and outputs. From Eq. (D-21), with pitch angle as 
the output quantity, 


(RF) es){w, Vv;6) = (RP) tee, Vv; 6) [1 + (RF) (es ot] (4-5) 


In Eq. (4-5), (RF )atw,/v;9 is the aircraft relating 
function for gust input and pitch angle output, and 
(RF) «so. is the open-loop relating function of the 
control system. It should be noted that the pitch 
angle output of the system is the sum of the pitch 
angle due to the elevator deflection plus the angle due 
to the vertical gust velocity input (see Appendix D). 
The system amplitude ratio is obtained from Eq. 
(4-5) by finding the magnitude of the relating function 
with A = 72. 

Fig. 10 is a plot of the normalized, nondimensional 
pitch angle power spectral density versus the non- 
dimensional frequency, 2, for three values of the 
pitch angle feedback sensitivity, S(cs)ig3. All other 
system constants have the values given in Fig. 6. 

For the case of zero feedback sensitivity—that is, 
when the aircraft is without a control system—the 
power spectral density curve peaks at a nondimensional 
frequency of 0.26, which corresponds to the damped 
natural frequency of the aircraft’s long-period mode. 
The effect of the short-period mode is seen in the slight 
distortion of the curve at higher values of ©. 

Introduction of the feedback sensitivity, Svcs)(6-3), 
greatly reduces the effect of the long-period aircraft 
mode, and for a value of S(es)(g.4; of 12 the power 
spectral density curve rises to a maximum at a non- 
dimensional frequency of 2.2. From this curve, the 
effect of either the long- or the short-period aircraft 
dynamics is not apparent. If the short-period aircraft 
dynamics influenced the power spectral density of the 
pitch angle to any extent, there would be a peak at a 
nondimensional frequency of 6.1, the damped natural 
frequency of that mode. 
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When the feedback sensitivity is further increased, 
the pitch rate control system approaches high-frequency 
instability, as discussed in Section 3. The effect of 
an increase of the feedback sensitivity on the system 
power spectral density is indicated in Fig. 10 by the 
peaking of the curve at a nondimensional frequency of 
130 for the case of S¢es){9.53 = 30. 

The integral of the power spectral density of the 
system output quantity over all frequencies yields 
the mean square of the system output quantity. 
Fig. 11 is a plot of the normalized mean square of the 
pitch angle and the normalized mean square of the 
pitch rate as functions of the control system pitch 
angle feedback sensitivity. An increase in the pitch 
angle feedback sharply reduces the value of the mean 
square of both output quantities. With S(es)i9.3) = 
12, the mean square of @ is 0.053 per cent of its value 
for the zero-feedback-sensitivity case. For larger 
values of S(.s){9.33, the mean square curve continues to 
decrease until the curve rises sharply near a sensitivity 
value of 36. This value of feedback sensitivity corre- 
sponds to the critical feedback sensitivity at which the 
pitch rate control system becomes unstable. The 
curve for 6’ behaves in a similar manner except that 
the curve rises more gradually for large values of 
Stcs)(a.4)- With Sceste.4) = 12, the mean square of 
6’ is 0.77 per cent of its value for the zero-feedback- 


sensitivity case. 


(V) EFFECTS OF AIRCRAFT ELASTICITY ON THE PITCH 
RATE CONTROL SYSTEM 


In Sections III and IV and Appendixes A and B, 
it is assumed when deriving the aircraft longitudinal 
equations of motion that the air-frame structure does 
not deform when elevator forces and moments are 
applied. Furthermore, the pitch angle feedback from 
the integrating gyro unit is assumed to be the same as 
that occurring at the aircraft center of gravity, re- 
gardless of gyro unit location. In this section, as 
supplemented by Appendix C, a method of including 
dynamic elastic effects is presented,’ and it is shown 
that elastic coupling may have detrimental effects on 
the stability of the pitch rate control system. 


(5.1) Method for Analyzing the Effect of Body Bending on 

Aircraft Performance 

The effect of body bending on the aircraft equations 
of motion can be included by employing normal modes 
of vibration as additional degrees of freedom. An 
arbitrary location on the fuselage is denoted by the 
symbol 7, and the linear displacement of the fuselage 
caused by the rth mode of vibration can be written as 
the product of two terms, ¢,‘"£,. The first term, 
¢;", is the ordinate of the rth normal mode shape at 
Station 7, and the second term, &,, introduces the time 
variation of the rth mode and constitutes the additional 
degree of freedom. 

The total linear displacement of station 7 along the Z 
axis, Z,, can be determined from the relationship 
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Z,=Zy-—l10+ Dd oe (5-1) 
r=1 

In this expression, the first term, Z,,, is the displace- 
ment of the aircraft center of gravity, where the center 
of gravity is not fixed in the aircraft as is the case 
with a rigid fuselage. The second term accounts for 
pitching of the aircraft about the center of gravity and 
consequently is proportional to the distance, /;, from 
the center of gravity to station 7. The summation 
in Eq. (5-1) includes all the normal modes of vibration 
and, in this manner, accounts for the effect of air- 
frame elasticity. 

The effect of fuselage bending on a measurement of 
pitch angle at station 7 can be determined by using 
Eq. (5-1) for 

dZ; ~  d¢," 
6,= —- =@-—) ¢, (5-2) 

dx r=1 dx 
It can be seen from Eq. (5-2) that at stations where 
dg;/dx is zero the local pitch angle, @;, equals the 
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Fic. 12. Effect of first fuselage bending on the aircraft elevator 
pitch angle relating function. 


The loci of the poles and zeros shown in this figure are deter- 
mined for an aircraft with a static stability coefficient ma = 
—1. Values of the other aircraft stability coefficients are tabu- 
lated in Tables 3 and 4 of Appendixes A and C. The loci are 
symmetrical with respect to the real axis. The loci in the lower 
half plane were omitted for convenience. The loci of the long 
period poles and zeros cannot be indicated to the scale of this 
figure. However for Qn, > 40 their location is virtually un- 
affected by first fuselage bending. By definition 


In this figure only the effect of structural elasticity on first 
fuselage bending is considered so that 


, m=4 
. méf c\* — 
(RF) ig = 4 ( ST Bnd” 
1 [050 8 \b, had 
m=0 
n=6 
* AnXr 
ws 
n=0 


It can be seen from this figure that with 2,; > 40 the high-fre- 
quency poles and zeros are nearly equal and the short-period 
zero and poles are unaffected by body bending and consequently 
(RF)Als;¢] = (RF)a4 [5:6] crigia), the relating function of a rigid air- 
craft. 





pitch angle, 6, of the complete aircraft. These stations 
are called antinodes. 

The normal mode shapes can be estimated analyti- 
cally but can be determined more accurately by vibrat- 
ing the full-scale aircraft on the ground. In addition, 
it can be shown that the time variation of each normal 
mode can be represented by a second-order differential 
equation of the form 


M,é, + Ce + Ke = YOO; & OnLy — Ly 


? 


where Q; is the force perpendicular to the aircraft 
X axis at station 7. The most significant forces 
perpendicular to the aircraft Y axis are applied aero- 
dynamically to the wing and tail as indicated in Eq. 
(5-3) where the subscript, w, refers to the wing and the 
subscript, ¢, refers to the tail. The coefficients, .1/,, 
C,, and K, of Eq. (5-3) account for the mass of the 
aircraft and the damping and elastic restraining forces 
acting on the aircraft. The damping coefficient, C,, 
is composed of viscous aerodynamic forces and the 
viscous equivalent of structural damping, and _ the 
elastic coefficient, A,, includes the effect of both aero- 
dynamic and structural stiffness. The normalized 
mass of the rth mode, J/,, can be evaluated from a 
knowledge of the normal mode shape and the mass 
distribution of the aircraft according to the relationship 


M, = > [¢: ]?Am (5-4 


The incremental mass, Am, is the mass of the fuselage 
at station 7. 

In order to determine the longitudinal response of 
an elastic aircraft, it is necessary to solve the force 
equations along the X and Z axes, the moment equation 
about the Y axis, and equations for each of the normal 
modes in the form of Eq. (5-3). These equations must 
be solved simultaneously, since the elastic deformations 
affect the aerodynamic forces and moments acting on 
the aircraft, and the aerodynamic forces introduced 
by the wing and tail influence the amount of body 
bending. This aeroelastic coupling is analyzed in 
detail in Appendix C. 

The complexity of the analysis obviously depends 
upon the number of normal modes required to describe 
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the aircraft motion. In most structures, the first 
and second normal modes are adequate, and in the 
example of this paper, only the first mode resulting from 
coupling between fuselage bending and wing torsional 
bending is considered, since it appeared to have the 
most effect on the pitch rate control system. In the 
discussion of this section the effect of this bending 
mode on the pitch angle, 6, is first considered, and then 
the effect of the bending on the pitch angle at any 
station, 9;, is investigated. 


(5.2) Effect of First Body Bending on the Pitch Rate 

Control System with Gyro Unit at the Antinode 

The longitudinal response of the aircraft with first 
bending included can be determined from Eqs. (C-6) 
(C-9) of Appendix C. The aircraft inputs in these 
equations are the elevator displacement, 6, the change 
in thrust, AZ, and the gust interference velocities, 
l; and Wg. The aircraft response is expressed in 
terms of the velocities, « and w, the pitch angle, @, 
and the displacement due to first body bending, &. 
The relationship between the various inputs and outputs 
can be evaluated by use of these equations and ex- 
pressed in terms of relating functions. For example, 
the elevator displacement—pitch angle relating func- 
tion, (RF ).4,5.9, can be written 


; -\2 4 6 
(RPlase = mm ey ~ ban/ > ad") 


Ss k, 0,1 n 0,1 


The coefficients of the polynomial in the denominator 
of Eq. (5-5) are the same for all relating functions 
derived from Eqs. (C-6)—(C-9). Values of \ for which 
this polynomial is zero are poles of the relating function 
and also are roots of the aircraft characteristic equation. 
The coefficients of the polynomial in the numerator 
of Eq. (5-5) are different for each relating function. 
Values of \ for which this polynomial is zero are zeros 
of the relating function. 

The zeros and poles of (RF).4(5.9,; are shown in Fig. 
12 for an aircraft with a static stability coefficient, 
mM, = —1. Values of other stability coefficients are 
listed in Table 3 of Appendix A and Table 4 of Appendix 
C. The loci of Fig. 12 are plotted as a function of the 
nondimensional undamped natural frequency of first 
body bending, 2,,. It can be seen that for values of 
2,, less than 14 the poles representing short-period 
motion have a positive real part, indicating instability. 
However, for values of Q,, greater than 40, the long- 
and short-period zeros and poles of (RF )a4j5.) are un- 
affected by first body bending. Moreover, the addi- 
tional two zeros and poles resulting from body bending 
are nearly equal, and consequently they do not in- 
fluence the relating function, (RF)aj5.9, since their 
effects cancel. For this reason the pitch angle relating 
function of the aircraft, (RF)a{5.91, equals the relating 
function obtained by assuming a rigid air frame, 
(RF) 415-6) (rigia)s or 


(RF) 415-6) = (RF) 45-6) rigia) (5-6) 


when ©,, is large compared to the short-period un- 
damped natural frequency. Consequently, for large 
values of 2,,, the pitch rate control system of Sections 
III and IV is not affected by body bending when the 
gyro unit is placed at the antinode. However, when 
the gyro unit cannot be placed at the antinode, the 
characteristics of the pitch rate control system may be 
influenced by body bending even for large values of 
q 


son.s 


(5.3) Effect of First Body Bending on Pitch Rate Control 
System with the Gyro Unit in an Arbitrary Location 
In order to determine the influence of body bending 

when the gyro unit is at station 7, it is necessary to 

obtain an expression relating the local pitch angle, 

6;, to the elevator displacement. This relationship 

can be determined by use of Eq. (5-2) and Eqs. (C-6)- 

(C-9) of Appendix C and for large values of {,, can be 

written in the form 


(RF) at5-0:) ad (RF), 6:6) (rigid x 


2s ?\ / 26 ? dis 
1+ A+ —J/(1+ > vA4+—. (5-7) 
Oy Qy?] | Q, | 


The poles of the relating function for the local pitch 
angle, 6;, are the same as the poles obtained from Eq. 
(5-5) and plotted in Fig. 12. For large values of 
Q,,, the poles can be readily determined as shown in 
Eq. (5-7), since they equal the sum of the poles of the 
rigid aircraft and the poles introduced by body bending. 
Eq. (5-7) also indicates the existence of a complex pair 
of zeros which only cancel the high-frequency poles 
when the damping ratio, fy, and natural frequency, 
Qy, of the numerator terms equal the damping ratio, 
6, and natural frequency, ©,,, of the first body bending. 
The numerator terms, (vy and Qy, can be determined 
from expressions contained in Appendix C, where it 
is shown that these terms depend upon the location 
of the gyro unit and only equal ¢ and ®,, at the anti- 
node. 

Fig. 13 shows the effect of body bending on the open- 
loop performance of the pitch rate control system and 
is obtained by use of Eq. (5-7). The curves are 
plotted for a pitch rate sensitivity ratio, S(es)(97.4)+ 
Sves)oa) Of 0.1, which provides the most satisfactory 
results for the rigid aircraft. It can be seen that for 
nondimensional frequencies in the region 30 to 100, 
the open-loop phase angle of the system is increased 
when the gyro unit is located aft of the antinode and 
is decreased when the gyro unit is forward of the anti- 
node. The decrease in the open-loop phase angle 
causes the pitch rate control system to become un- 
stable when the open-loop sensitivity, S(cs)¢9.3, has a 
value of 12 and the gyro unit is forward of the station 
corresponding to case C of Fig. 13. 

It can be seen that the elasticity of the fuselage 
may have an extremely detrimental effect on the 
performance of the pitch rate control system. This 
effect can be alleviated by (1) locating the gyro unit 
close to the antinode; (2) lowering the open-loop 
sensitivity, Si-s){g.4)3 (3) introducing a filter to lower 
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Fic. 13. Open-loop frequency performance of the pitch rate 
control system for various gyro locations. The aircraft bending 
characteristics are presented in Appendix C. The curves are 
computed for a first bending frequency of the fuselage Qn; = 8O0 
with ¢; = 0.045. The effects of all other bending frequencies are 
assuined negligible. The gyro locations are: 


Case Gyro Location ty Qn 
a tail 0 034 76 
b antinode 0 045 SO 
c 17, distance from the 0.060 85 

antinode to the nose 
d nose 0.099 101 


the frequency at which the phase angle, (PA). 02, 
is 180°; or (4) stiffening the structure of the air frame, 
thereby raising 2,,. 

Method (1) is satisfactory when space exists for the 
gyro unit at this location and when the second body 
bending frequency is sufficiently high. 

Method (2) is usually the most convenient method 
for insuring system stability, provided the response 
time does not become appreciably longer and provided 
the decrease in open-loop sensitivity does not induce a 
low-frequency instability for aircraft that are statically 
unstable, as presented in Fig. 9. 

Method (3) decreases the open-loop phase angle and 
amplitude ratio and causes an increase in gain margin. 
However, the phase margin is necessarily reduced for 
aircraft that are statically unstable. 

Method (4) permits the control system response of 
Sections 3 and 4 to be retained but usually involves an 
appreciable increase in weight of the aircraft structure. 


CONCLUSIONS 


The authors have reached the following conclusions, 
which are based on the material presented in this paper 
and their personal experience in the field of aircraft 
control. 

(1) An integrating gyro unit is ideally suited for 
use in a pitch rate control system. 

(2) The pitch rate control system has a response 
that is nearly independent of the static stability 


coefficient for —1 < m, < +1, provided the limits of 
the integrating gyro unit, the elevator servo, and the 
elevator are not exceeded. 

(3) The pitch angle feedback sensitivity of the 
pitch control system has both an upper and a lower 
limit when the aircraft is statically unstable. 

(4) The long-period motion of the aircraft has a 
noticeable effect on the response of the pitch rate 
control system but does not influence the system 
stability. 

5) The pitch angle deviation of the aircraft re- 
sulting from vertical gust velocities can be reduced 
appreciably by the tight-loop stabilization of the pitch 


rate control system. 

(6) The pitch angle response at the antinode is not 
affected by first-body bending when the undamped 
1atural frequency of bending is large in comparison 
with the short-period natural frequency of the aircraft. 

(7) First-body bending can have an appreciable 
effect on the response of the pitch control system when 
the integrating gyro unit is not at the antinode, and 


particularly when it is in a forward position. 


APPENDIX A 


NONDIMENSIONAL EQUATIONS OF AIRCRAFTI 
LONGITUDINAL MOTION ASSUMING RIGID FUSELAGE 


The rigid aircraft's longitudinal dynamic equations 
of motion are written in the manner and follow the 
notation of Durand (see p. 129 of reference 12). The 
aircraft axes Y, VY, and Z have their origin at the 
center of gravity of the aircraft with Y normal to the 
plane of symmetry (positive along the right wing), 
X parallel to the projection on the plane of symmetry 
of the velocity vector in trimmed flight (positive for- 
ward), and Z forming a right-hand system. The veloc- 
ity vector is a measure of the speed and direction of 
the aircraft with respect to the air mass. 

The equations have been linearized, permitting 
investigation of incremental deviations from steady 
flight conditions—i.e., from ‘‘trim’’ conditions—along 
horizontal, climbing, or diving paths. The wings are 
assumed to be level. The incremental variations in 
aircraft velocity components “ and w, from the as- 
sumed trim values Ll’) and IW along the X and Z axes 
respectively are defined by 


U = Up + u (A-1) 
V= Ww+w (A-2) 
where WW = 0 by virtue of the direction chosen for 


the X axis. Similarly, pitch angle 6 and climb angle 
lr in terms of the trim and incremental values are 
defined respectively by 


6=8®+80 (A-3) 


r=Mt+y (A-4) 
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Eq. (A-5) defines the angle of attack, a 


a = tan—'[(w — we)/(U — Ue 
(A-5) 


"oo 


a ‘ r 
a = (& We) U5 


where vq and we are the Y and Z components of the 
gust velocity measured with respect to the local air 
mass. Eq. (A-6) gives the geometric relationship 
petween the pitch and climb angle variables 6 and y. 


y = 6-—tan-'(w/L)]| 
- A-6) 
4 = 06 — (w Us) \ 


/ 


The linearized longitudinal force and moment equations 


are assumed to be as follows. 


ma = [(u — ug)Xu)] + [(w—- wea)Xu] + 
6X, + AT = ma, (A-7) 


—m(V6 — w) = [(u — ug)Zy) + [(w — we)Z,) + 


u 


0Z, + 6Z; = ma; (A-S) 


/6 = Ree ug) M, + (Ww — We) MM, | + 
(w — we) M;] + 6M, + 6M, (A-9) 


Simplification of Eqs. (A-7)—(A-9) is accomplished by 
assuming that A7, Z;, and .1/, are negligible and by 
use of the nondimensional time variable, I’t/uc, and 
the following nondimensional stability derivatives. 


vw = (V/qgS)X, = (V ¢S)(0X/ou (A-1C) 

xe = (V/gS)Xe = (V/gS)(OX/dw)  (A-11) 

zs, = (V/qS)Z, = (V/gS)(0Z/du) — (A-12) 

Sa = (V//qS)Zy = (V/gS)(0OZ/Ow) (A-13) 

% = (1/¢S)Zo = —(mg/qS) sin Oy = 

—C, tan Q@ (A-14) 

ve = (1/gS)Xe = —(mg/gS) cosO@, = —C,  (A-15) 
Me = (1/qSc)(2V/c)(OA/ Od) (A-16) 

m, = (1/qSc)(2V/c)(OM 06) (A-17) 

m, = (1/g¢Sc)(0M/da) (A-18) 

ms = (1/qSc)(OAl/05) (A-19) 


With d [d(Vt/uc)} denoted by primes, the nondi- 
mensional versions of Eqs. (A-7)-(A-9) become, 


respectively, 


2(u'/V) — [x,(u/V)] — [x Q(w/V)] — x6 = 
—[xu(ute/V)] — xe(we/V) (A-20) 
— [,(w/V)] — f2.(u/V)] + [2(w’/V)] — 20’ = 
—[s,(ue/V)] — 24(we/V) (A-21) 
= % 3 — nh. ; +2 (**): 6” — “y 0’ = 
We Mg We’ 
um; — wm, yo > 77 (A-22) 


The rigid-body parameters listed in Table 3 were used 
in the investigations of this paper and are typical of a 
high performance aircraft in level flight—i.e., @) and 
Ty are both zero. 
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> 


TABLE 3 
Dimensionless Aircraft Parameters for Horizontal Flight 


x = —0.048 z = —() SO Ma = -— 10 
ta = 0. 220 Sa = —4 SU Ma = — 36 
vt, = —0.400 Se = 0 m = —22.0 

m = — 1.54 
uw = m/pSc = 171.5 uc/V = 3.8 sec c/k, 0.67 


APPENDIX B 


RELATING FUNCTIONS FOR AIRCRAFT WITH RIGID 
FUSELAGE 


The linear differential equations of Appendix A are 
transformed into functional relationships convenient 
for the purposes of this paper by use of the Laplace 
transformation. The Laplace transform of a general 
quantity, g(/), in terms of the dimensional Laplace 
operator, p, is defined by 


O(p) = q(t)e"dt B-1) 
In terms of nondimensional time, I7¢/uc, and the 
corresponding dimensionless Laplace operator, \ = 
(uc/ V)p, the definition of the Laplace transform 


becomes 


4 lt ' Vt 
ed | MC uc 


O*(r = f V Le OP) ) 
and (B-3) 
g*(Vt/uc) = gid) \ 


where 


The relating function for any linear system relates 
any output quantity, gout, to any input quantity, din, 
all other inputs and initial conditions assumed zero, 
and is defined as the ratio of the Laplace transforms 
namely, 


c 


(RF) ¢eyststein:¢outl = . ™ 
( (s} in'%out! Qin(p) G* (r) (B-4 


Qout(P) OF lA) 


in 


Elsewhere in this paper, for convenience, asterisks are 
omitted, and both dimensional and dimensionless 
Laplace transforms are denoted by lower case symbols. 

Once transformed, the equations of motion listed in 
Appendix A may be solved by expanding the determi- 
nant of the coefficients and applying each of the inputs, 
one at a time. This procedure gives the following 
aircraft relating functions. 


c\2 
(RF) ats-01 = ‘(5 ms X 


[(XuS—~ — XeSu) — 2(Xy + Za)A + 4A? — 
: , (Bb-.)) 
(ao + QA + aod? + azd* + A‘4) 
P A /c\? 
(RF Alu, v:¢) = x 
7 S k, 
[2um.x, + (x,me — 4um,)A + (—2m,)d?*) - 
a ( -)) 


(do + aX + Qo\? + a3A3 + A*) 
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The coefficients of the characteristic equation of the 
aircraft relating functions are 


ao = —(u/8)(c/Ry) *XpeumMe (B-7) 


ay = +[(1/8)(c/k,)?]{x,[2um, — (1/2)(s,m4)| + 


((1/2)z,(mgxq — xem,)|}  (B-8) 


tum, + x,(mg + me)| + 
[(1/4)(wyZe — 2uX%e)] (B-9) 


a2 = [(1/8)(c/k,)?|[z m6 — 


—[(1/4)(c/ky)?](mp + mz) — [(1/2) (Ze + Xu)]) 


(B-10) 


APPENDIX C 


EQUATIONS OF AIRCRAFT LONGITUDINAL MOTION 
INCLUDING Bopy BENDING 


In order to determine the time variation of the 
normal modes of vibration from Eq. (5-3) of Section V, 
expressions must be developed for the lift of the wing 
and tail in terms of the dependent and independent 
variables. However, the wing lift depends upon the 
angle of attack of the wing, @,, and the tail lift upon 
the tail angle of attack, a,, and the elevator displace- 


ments, 6, since 


Lo = (dC./da) wqSa, (C-1) 


and 


L, = (dC,/da) mgS,a, + (dC. /d6) ng 6 (C-2) 


The angles of attack of the wing and tail are de- 
pendent upon the fuselage bending and wing twisting 
and bending in the form of coupled modes. In addi- 
tion, the downwash at the tail is proportional to the 
wing angle of attack (average span value) which causes 
additional coupling between body bending and the 
aerodynamic forces acting on the aircraft. By making 
use of Eqs. (5-1) and (5-2) it can be shown that 


Qy = a— >, é,(do”/dx), (C-3) 


and 
de de ) (;*) : 
Yi = [= a + = a oh - 6— 
. ( <) (= ] t J l 
(Tr) (Tr) 
> | (“ )| _ de (** ) 4 
- dx » @€ar ds, 
- be 1d dd”\ de 
ret | - ( -| (C-4) 
: ViL le dx J» da 
In Eq. (C-3) the average velocity of the wing in the Z 
direction is neglected since /,,/V and ¢,/V are small. 
In Eq. (C-4) /, is the distance from the tail to the center 
of gravity, V’, is the tail velocity, and de/da is the ratio 
of downwash at the tail to the angle of attack of the 
wing. 
The time variation, £, can be determined from Eqs. 


(5-3), (C-1)-(C-4). In cases where the body bending 
is dominated by the first mode it can be shown that 
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¢.” / ¢ 
1) 4 o¢0,, (8 2,,2(~) = 
Cc ak * Se 


ba + bga’ + bj’ + b56 (C-5) 


Eq. (C-5) is written in nondimensional form by ex. 
pressing time derivatives with respect to V,/uc, as 
indicated by the primes, and by expressing the bending 
motion as the ratio of the time variation, £,, to the wing 
chord, c. The damping ratio, ¢;, and nondimensional 
undamped natural frequency, &,,, of the first bending 
mode are nearly equal to the values of these quantities 
measured in still air. The bending coefficients, },, 
bx, bo, and 6;, are functions of the fuselage flexibility, 
the wing and tail characteristics, and the dynamic 
pressure. 

For elastic aircraft, the changes in wing and tail 
angle of attack resulting from body bending may 
influence the short-period motion of the aircraft. In 
these cases, it is necessary to solve simultaneously 
Eqs. (A-20)—(A-22) modified to include body bending 
and Eq. (C-5) in order to determine the motion of the 
aircraft. These equations can be written in terms of 
the complex operator, A, in the form* 


V) — x,(w/V) — x6 = 


(AT/qS) — xXy(ue/V) — x4(we/V) 


Ly (*) + (—2 + z,) (*) + 20+ % (*) = 
I J “"1\¢ 


(*") re =) C7 
, z (C-( 
\y a\7 

Lm; w 
“(tas n)(s)s 

2 2 | 


A k, . l 1 

2 A— ~-m,|O0+ (mir + m; ) = 
m Cc 2 i . Cc 
s) (4°) (C-8 


—(b,d + ba) (*) ~ bdO + (2 + 26,0, + 


€ WG 
2,,,”) (*) = b56 sat (b, oe bX) (¥) (C-9 
c 


The elevator displacement-pitch angle relating function, 
(RF ats}, Of Eq. (5-5) can be derived directly from 
Eqs. (C-6)—(C-9)._ The zeros and poles of this relating 
function are plotted in Fig. 12 for an aircraft with the 


(2) — x,)(u 


(C-6) 


ctr 


Ee oa ( 
—- AT — |m, + 
gSC 


lm, 


m;6 + : 
<p 


bending characteristics of Table 4. 

A relating function for the local pitch 
(RF )ats-03, can be obtained from Eqs. (C-6)—(C-9) by 
the additional use of Eq. (5-2). It can be shown in 


angle, 


*In this form the variables (u, w, 0, &, weg, ug, AT, and 6 


are functions of A, the dimensionless Laplace transform operator 


TABLE 4 
Dimensionless Bending Coefficients of Elastic Aircraft 
zt, = 2.424 ba = —154.6 
mi, = 0 593 ba = — 5.6 
mi, = 0.040 bo = — 9.4 
o& = 0.045 bs = — 41.2 
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RECENT DEVELOPMENTS IN AIRCRAFT 


the example selected for this paper that the natural 
frequency of the first bending mode is sufficiently 


great that 


(RF 1(6.0,) jelastic = (RF), 6,9) |rigid x 
r 2g ?7\ / 2%, ? : 
| i+ a+ — 1 + 1+ — (C-10) 
Qy Qy*] | Q, ie.” 
where 
Qy* j > dg; )} . ( ' ) ‘ 
= 1+? > 1 — fo (C-11) 
Ov, | i i (“4 C) ,/ \ ~ \d(x/c) 
and 
(y j ; dg; ( ] , dg; , | i 
ee ee | -- fF , 1 — fo (C-12 
P ) r ‘ | = | / \ J2 | c) ) 


It can be seen from Eqs. (C-10)—(C-12) that at the 
antinodes where d@,‘'’/dx equals zero, the pitch angle 
of the aircraft is unaffected by body bending. The 
factors, fi, fo, and f3, are positive quantities so that 
at other stations the undamped natural frequency 
and damping ratio of the numerator quadratic, Qy 
and ¢y, are increased by a positive increase in d@;'"/dx. 
The curves of Fig. 135 are determined from Eq. (C-10) 
for the case where f; = 0.031, f2 = 0.825, and f; = 1.70 
and with [d¢,'" /d(x/c)] equal to —0.127 at the tail, 0 
at the antinode, +0.146 at 25 per cent of distance from 
antinode to the nose, and +0.437 at the nose. 


APPENDIX D 
RELATING FUNCTIONS OF PITCH RATE CONTROL 
SYSTEM 
(1) Composition of Pitch Rate Control System 


The pitch rate control system of Fig. 2 is composed 
of (a) an amplifier with 


Sig = Sample (D-1) 
(b) an integrating rate gyro in which 
My = Sty (D-2) 
I,A, + C,A, = M, + HO (D-3) 
eg = 3A, (D-4) 
or 
(*) &é+ ,= (**) Siig + (=) So (D-5) 
a2 : e, C, 
(c) a signal modifier where 
Tamsm + sm = CeSism)leciegm) 1 
(@7 + OsmTsm&g)S(sm)leg:e5m) (D-6) 
and (d) an elevator servo where 
6 H+ 2€yansb + wns7d = asic (D-7) 


Combining Eqs. (D-1)-(D-6) and assuming J,/C, 
nearly zero gives in operational form 
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HSS sm) leg ;€¢m lL + QemTemP 
Csm(p) = — |} Ob) + 
l nP 


C,p + Tes 
Lo ee 
C,p ( Ll + TompP Je pr 

Ss 
e.(p (ID-S) 

Ll + T.mp 

Then by letting 

Stomiiege..) = —Senp(Seg/Co)Suqh tees lL) ron . 
(D-9 


Eqs. (D-7) and (D-S) yield in nondimensional form 


6(A) 7 | x 
(H/C) SwgSemtegseepSs 1 + [(2Es/Qns)A] + (A2/Q, 
j [—(uc/V)Samp( S/H 
\~ . é(A) 4 


l + (Va mT sm/ Me 


: a’(n)' (1D-10) 
ALL + (Vrem peor / 


where 
Qns = (uc/V)w (D-11) 
A = (uc, V)p (1-12) 
0’ = (uc/V)O (D-13) 


(2) Mathematical Representation of Pitch Rate Control 
System 


It can be seen from Fig. 3 that 
ba) 


6(A) l 
: = — " ee 
Ses) [9:8] [1 + (26;/Qns)A + (A?2/Qy5?) ] 
f  Sves)te,:61 [1 + (Scesyte-01/Stesvie-s)) .,..) 
4,-— e(A) + 6° (X) 
( N M1 + (Vrem/peya 
(D-14) 


A comparison of Eqs. (D-10) and (D-14) shows that 


Ses) [6:8] = (7 CL) Sa Gn [ps€smj 9s (D-15) 
S cs)(0';6 Ses) (0;8] VasmTsm KC (D-16) 
Sisstesn « ~Ge/VV aS) (D-17) 


However, the aircraft is influenced by both the 
elevator deflection and gust disturbances, or 


6’(A) = 4(A)(RF )ateor + [wol(A)/VICRP) atieg v0" 


(D-18) 

Solving Eqs. (D-14) and (D-18) gives 
0’(X) = e(A)(RF) ceste.01 + [Wel(A)/V (RF) ces) tug /v 01 
(D-19) 

(RF) ces) te,:6'1 | (RF) cesyon ia 
Stesyicce) L1+ (RP)coor 
l (Vr5m/pmc)r 
| > = . | (D-20) 
] + (Se [0-6 S cs) (0:6 rv 
(RF) atuc/v0’ 

Alng/¥ (D-21) 


F (es)[wg/V;e') = 
(RF) G/V 0) 1+ (RF) csyo1 
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1 + (S¢se)(0":8)/Stserfo;8)A 
1 + (Vr5m/pe)r 
(RF ) ats.) 
[1 + (26,/Qm)¥ + (A2/Qys?)] 


(RF) (esyox 
S(se) (0:8) 


(D-22 


(3) Deviation of the Aircraft Due to Gust Disturbances 


In general 


it? 
f [((PSD) qin ](A R)* es) tain :dould@ 
vis 0 


| ” | CRISD ai. 
= f q (A Ie) cs) dou JAQ 
rJ0 pec/t 


(D-23) 


( ALS) qout = 


Consequently, 


l "Tl (PSD) (we/V) |, 
(\S)0 = { = (AR) * (cs) (wg v.gd2 
0 


T pe/ | 
(D-24) 
1 {(°: (PSD) (we/V) 
(MS)0’ = [ | = (AR) * (es) tw¢ v9 d2 
rJ0 1 | 
(D-25 
where 
(PSD) (we/V) a 2[ Vre/uc|(MS)(we/ V) (9-90) 


pc/V [+ 1(Vre/pc)Q2}? 


APPENDIX E 
DEFINITION OF TERMS 
Many of the parameters used throughout the paper 
have been defined when introduced into the discussion. 


As a consequence, only the more important terms are 
listed in this appendix. 


GENERAL DEFINITIONS 
(RF) syst.) a = system relating function, defined as the 
ratio, Jout(p) to gin( p) 
(A R)oyst lain'@our! = ~‘System amplitude ratio 
(PA )gyst leintfout = system phase angle 
:. leis teout system sensitivity 
[(PSD)q] = power spectral density of g 
[(.12S)q] = mean square of g 
Vr/pc = nondimensional time constant 
Vt/ye = nondimensional time 
p = dimensional Laplace operator 
d = nondimensional Laplace operator 
d/dt = differentiation with respect to real time 
denoted by dots 
d 
d( Vt/uc) = differentiation with respect to non- 
dimensional time, denoted by primes 
w = dimensional frequency variable 
Q = nondimensional frequency variable 
OL = subscript indicating loop opened 
Ricip AIRCRAFT PARAMETERS 
0 = angle between the aircraft trim axis 


and an inertial reference axis (zero 
at trim) 


MARCH, 1955 


4 = angle between aircraft velocity vector 
and an inertial reference axis (zero 
at trim) 

= aircraft elevator angular displacement 


‘zero at trim) 


V=l = average trim air speed 

u = perturbation speed along aircraft Y axis 
% = perturbation speed along aircraft Z axis 
A = subscript referring to aircraft 


\ oe = orthogonal set of aircraft axes; also, 
with subscripts, the dimensional force 
coefficients corresponding to the in- 
dicated aircraft axis 

M = dimensional moment coefficient, mo 
ments about VY 

KX)» Be) = nondimensional force coefficient corre- 
sponding to indicated aircraft axis 


m = nondimensional moment coefficient, 1 


10- 
ments about } 


m = aircraft mass 


SYMBOLS FOR NONRIGID AIRCRAFT 


o = mode shape of rth mode at the ith sta- 
tion on the aircraft 

é, = rth harmonic mode of vibration 

do; /dx = slope of the rth mode at the /th station 

a. = the nondimensional natural frequency of 
the rth mode 

¢ = damping ratio of the rth mode 

b = normalized force coefficient, rth mode 
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Summary of Recent GALCIT Hypersonic 
Experimental Investigations 


HENRY T. NAGAMATSU?# 
California Institute of Technology 


ABSTRACT 


The GALCIT 5- by 5-in. continuous hypersonic wind tunnels 
are designed for maximum stagnation temperatures of 400°F. 
Leg No. 1) and 1,100°F. (Leg No. 2) and for a maximum reser- 
voir pressure of 1,040 Ibs. per sq.in. The hypersonic shock tube 
is being used to obtain aerodynamic information, and some pre- 
liminary results are presented. By measurements of static and 
total head pressure in the nozzle, the amount of supersaturation 
of the air has been determined for a higher range of hypersonic 
Mach Numbers. On a 26-in. long flat plate an extensive in- 
vestigation at Mach Number 5.8 was conducted to determine the 
contamination of the laminar boundary layer due to air injection 
from orifices near the leading edge and the turbulent boundary 
layer from the side walls. The transition Reynolds Number for 
a flat plate was at least 5.5 million for a Mach Number of 5.8, 
and the laminar boundary layer seemed to be extremely stable at 
this Mach Number and required strong disturbances to produce 
transition. The static pressure on a flat plate varied due to the 
interaction of the leading edge shock wave with the boundary 


layer on the plate. 


INTRODUCTION 


FT ppuseneageesian PROBLEMS at the hypersonic speed 
range have become more critical with the de- 
velopment of efficient rockets and missiles to operate 
at high Mach Numbers. Previously, one of the out- 
standing missiles was the V-2, which had a maximum 
velocity of approximately 5,000 ft. per sec. and a maxi- 
mum range of about 250 miles. In the spring of 1949 a 
two-stage rocket with V-2 (first stage) and JPL Bumper 
WAC (second stage) was launched at White Sands 
and attained an altitude of about 250 miles with a 
maximum Mach Number at the end of burning of 
about 7.5. For the long-range ballistic type of rocket 
it is necessary to have high velocity at the end of burn- 
ing in order for the rocket to travel long distances. 
For a range of about 6,000 miles, a maximum Mach 
Number of about 20 is necessary; and for orbitting 
around the earth, a Mach Number of about 25 is re- 
quired. Thus, it is apparent that fundamental knowl- 
edge regarding hypersonic flow phenomena is neces- 
sary in order to enable the designer to construct ef- 
ficient long-range rockets and missiles. 

There are various methods of obtaining hypersonic 
aerodynamic information, and with each technique 
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nual Meeting, IAS, New York, January 25-29, 1954. Revised 
and received October 25, 1954 
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there are advantages and disadvantages. One of the 
methods of hypersonic investigation is the use of hyper- 
sonic wind tunnels operated with air at elevated tem- 
peratures and pressures or operated with other gases, 
such as helium, to attain high Mach Numbers. An- 
other method that has been used extensively for the 
supersonic speed range is the ballistic range and free- 
flight firing ranges. The advantage of this system 
over the wind tunnel is that aerodynamic information 
is obtained in the actual free-flight conditions, but 
the disadvantages are the difficulties in obtaining aero- 
dynamic data. Few data have been obtained in the 
ballistic range at Mach Numbers greater than 5 be- 
cause of the problems involved in firing the models at 
higher speeds. 

More recently, the shock tube has been developed 
not only to obtain high stagnation temperatures but 
also to obtain hypersonic Mach Numbers with high 
stagnation temperatures by expanding the gas in a 
contoured nozzle. With this technique it is possible 
to obtain free-flight stagnation temperatures in a labo- 
ratory for a short duration of flow. Recent analytical 
investigations of the boundary layer at hypersonic 
speeds!~‘ have indicated that the laminar skin friction 
may increase with the Mach Number for a given Rey- 
nolds Number instead of decreasing as predicted by 
the classical theory. This increase in the skin friction 
is caused by the interaction of the boundary layer with 
the leading edge shock wave. There is an urgent 
need for experimental data at hypersonic Mach Num- 
bers to check the validity of the assumptions used in 
these analytical investigations. 

A brief summary of the recent experimental results 
for the GALCIT hypersonic research program spon- 
sored by the Army Ordnance and the Air Force will be 


presented. 


DESCRIPTION OF THE GALCIT HyPeErRSONIC FACILITIES 


The two GALCIT 5- by 5-in. hypersonic wind 
tunnels (Legs No. | and 2) are operated from one com- 
pressor plant consisting of thirteen Fuller compressors 
and three reciprocating-type compressors for the last 
These compressors can 
the 
desired compression ratio for any Mach Number set- 
For Mach Numbers greater 


two stages of compression. 
be arranged in various combinations to obtain 


ting in the test section. 
than 7, six stages of compression are necessary, and the 
maximum design pressure for this system is 1,000 Ibs. 


per sq.in. Recently, a seventh-stage compressor has 
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- he controlled within +0.05 Ib.per sq.in. by means 
. : weaver aa wet of a Minneapolis-Honeywell-Brown circular chart con- 

. ong) £. ofventt m4 troller. 

ra a ¥ The stagnation temperature is measured by a ther- 
Ft _- . ROPLE: mocouple probe located upstream of the nozzle throat. 
a sida This temperature is recorded and controlled by means 
se) z 7 ' a : F of a Minneapolis-Honeywell-Brown circular chart 
ae a nd EE ike “Gl a ee Se controller. This controller is used to control the steam 
Me gagt < t Satis P { be heat exchangers for Leg No. 1 operation, and, for Leg 
x — a ee No. 2, it is used to control the reactor that changes the 
. @v current input to the electrical heater. All the tem- 
SCHEMATIC DIAGRAM 2 aati peratures necessary for the plant operation are indi- 


OF GALCIT 5x5in HYPERSONIC WIND TUNNEL INSTALLATION 


FIG 


been installed to increase the maximum stagnation 
pressure to 150 atmospheres and to increase the operat- 
ing range of the facility to a Mach Number of about 
15. 

Leg No. 1 is designed for a Mach Number range of 
2 to 7 with a design stagnation pressure of 1,000 Ib. 
per sq.in. at S60°R.; Leg No. 2 is designed for a Mach 
Number range of 4 to 11 with a maximum design stag- 
nation pressure of 1,000 Ibs. per sq.in. at 1,560° R. 
Both tunnels are of the continuously-operating, closed- 
return type with provision for heating the air to obtain 
the desired stagnation temperatures. For Leg No. 2 
the steam heat exchanger is used as the first-stage 
heater, and a _ specially-designed 300-kw. electrical 
heater in the reservoir is used as a second-stage heater 
to heat the air to 1,560°R. 

The compressors and all valves are operated re- 
motely from a master control panel located adjacent 
to the test section. The moisture in the air is re- 
moved by a 2,200-lb. bed of silica gel in the main air 
circuit, and the silica gel is reactivated prior to each 
run by a built-in blower-heater-condenser system. 
The maximum water content of the air is kept below 
50 parts per million by weight at all times. The dew 
point is measured at the beginning and end of each run 
with the use of a standard type carbon dioxide-cooled 
indicator. 

At the present time, Leg No. 1 has a fixed nozzle 
designed for a Mach Number of 5.9, and Leg No. 2 
has a nozzle designed for a Mach Number of 8. The 
throat heights for these Mach Numbers, with test 
section dimensions of 5- by 5-in., are 0.060 in. and 0.020 
in., respectively. These nozzle blocks were designed 
by the Foelsch method with corrections made for the 
estimated boundary-layer displacement — thickness. 
Static pressure orifices at 1-in. intervals in the top and 
bottom nozzle blocks permit a check to be made with 
the calibration on each run. It is possible to vary the 
Mach Numbers in each of the nozzles by adjusting the 
throat height. 

Vacuum-referenced manometer boards using DC-200 
silicone fluid are used to measure all static pressures, 
and an atmospheric-referenced mercury micromanom- 
eter is used to measure the total head pressure. The 
tunnel stagnation pressure is measured with a Tate- 
Emery nitrogen-balanced gage, and this pressure is 


cated on a 20-point Leeds and Northrup recorder, 
The overall plant schematic indicating the heating 
system is shown in Fig. !, and the test section for Leg 
No. 2 is shown in Fig. 2. 

An optical system using a BH-6 steady source is used 
for schlieren and shadow photographs of the flow 
around the models. For heated runs the steel test 
section doors are fitted with circular plate glass observa- 
tion ports. 

In order that aerodynamic phenomena related to 
hypersonic flows might be studied at stagnation tem- 
peratures simulating free-flight conditions, a nonuni- 
form shock tube with divergent sections? was con- 
structed similar to the one discussed in reference 6. It 
is possible to undertake many types of fundamental 
investigations with the shock tube, but at the present 
time only meager data are available considering the 
possibilities of this instrument. Therefore, our in- 
vestigations have been conducted along the following 
lines of basic research: development of instrumentation 
to get useful and reliable data during uniform flow of 
extremely short duration; methods of producing 
theoretically-predictable flows at high Mach Numbers 
and stagnation temperatures; study of the deviations 
from ideal shock tube theory at high pressure ratios; 
investigation of flow over simple bodies at high Mach 
Numbers with high stagnation temperatures. 

The shock tube (Fig. 3) consists of a cylindrical 
compression chamber made out of 3-in. diameter, 
alloy-steel pipe 8 ft. long, capable of withstanding 
2,000 Ibs. per sq.in. The 3-in. compression section is 
connected by a smooth contraction section to a 1- by 
3-in. uniform-expansion chamber 11 ft. in length. For 
the diaphragm, acetate plastic sheets are used _be- 
tween the compression and the expansion chambers. 

In the uniform expansion section two observation 
stations are placed 7 ft. and 11 ft., respectively, from 
the diaphragm. At these stations it is possible to 
place either piezo electric pickups or glass windows to 
take schlieren or shadow pictures of the flow. The 
uniform expansion chamber is connected to a diver- 
gent section, which starts with a cross section of 1 by 3 
in. and expands to 40 by 3 in. at a distance of § ft. 
from the throat. In this region the flow is expanded 
after the shock wave to obtain a high Mach Number 
in the test section. The nozzle has a straight section 
15 by 3 in. for a distance of 5 ft. after the test section. 
A 7-in.-diameter test section window is located 7 it. 
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from the throat and is used for the study of the flows 
over the various models by either shadow or schlieren 
techniques and for the investigation of the pressure 
in this free stream or on the surface of the models. 
At the end of the nozzle there are valves that prevent 
the pressure from building up in the test section from 
the reflected shock waves. With the present state of 
instrumentation, the following investigations are being 
conducted: calibration of pressure pickups; speed of 
propagation of shock waves; time histories of the 
static or total pressures at various stations along the 
uniform and the divergent sections; time history of 
steady flows by means of schlieren and shadow meth- 
ods. Various Mach Numbers can be produced by 
varying the pressure ratios between the compression 
and the expansion regions and by varying the nozzle 
expansion ratios. 

The gas combinations used for the investigations 
so far have been limited to nitrogen-air and helium-air. 
In the future, other combinations of gases will be used 
to investigate the effects of high stagnation pressures at 
high Mach Numbers. At low pressure ratios across 
the diaphragm the ideal shock-tube theory’ was found 
to be in excellent agreement with the experimental 
results. Strong shock waves have been produced by 
evacuating the expansion chamber down to 3 mm. 
of mercury, and by this method pressure ratios up to 
18,000 have been obtained. 

The speed of propagation of the shock wave was 
found to be in good agreement with the ideal theory 
up to these pressure ratios. The pressure ratios 
across the shock wave were slightly lower than predicted 
by theory, and the duration of the steady flow in the 
uniform section was much shorter than expected. 

For a strong shock-wave condition, with helium as 
the driving gas and air as the working flow, a notice- 
able glow of the air was observed both in the uniform 
expansion chamber and at the test section around the 
models. The calculated stagnation temperatures when 
the glow was noticed were of the order of 6,000°R. to 
7,000°R. By utilizing the double diaphragm tech- 
nique, it is possible to increase the stagnation tem- 
perature by a factor of two to three without undue 
difficulty. Thus, it will be possible to simulate the 
stagnation temperature that is encountered in flight 
at high Mach Numbers. 

The Mach Number of the flow in the test section, 
calculated both by the shock angles and angles of 
attack for the detachment of the shock wave, have 
indicated a Mach Number of about 6 for the particular 
nozzle now installed, which is close to the theoretical 
predicted value. For the present operating range of 
the shock tube, the static pressure in the test section 
is extremely low, and the crystal pickups are not ac- 
curate enough to measure pressures of the order of 
0.15 mm. of mercury. Means of improving the ac- 
curacy of these measurements are now being investi- 
gated. The duration of a steady flow in the test sec- 
tion starts at around 1,000 microsec. at a pressure ratio 
of 2,500 and steadily increases to beyond 2,000 micro- 
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sec. at pressure ratios above 10,000. The results 
of the pressure pickup and the schlieren studies have 
indicated that the flow during this steady period js 
highly uniform. The preliminary results from the 
hypersonic shock tube have indicated that it is an ex- 
tremely useful instrument for obtaining high stagna- 
tion temperatures at high Mach Numbers. One of 
the serious disadvantages of the shock tube for hyper- 
sonic investigation is the extremely low Reynolds 
Numbers that are encountered in the test section. 
This difficulty can be overcome by using double dia- 
phragm technique with high pressures in the com- 
pression chamber. The possibilities of the hypersonic 
shock tube for obtaining fundamental aerodynamic in- 
formation with dissociation and ionization of the air 


are promising. 
CONDENSATION OF GASES IN HYPERSONIC NOZZLES 


In a hypersonic nozzle, where the air is expanded 
isentropically from the stagnation condition in the re- 
servoir to the high Mach Number in the test section, 
the decrease in the temperature of the air as it ex- 
pands is a function of the area ratio or the Mach 
Number. For high Mach Numbers, as the air cools 
in the nozzle from the throat, the free-stream tem- 
perature and pressure will cross over the condensation 
point, and the gas will eventually condense. For ex- 
ample, if the air is assumed to be perfect and isen- 
tropic, the temperatures of the air for Mach Numbers 
of 5 and 10 are 50.1°K. and 14.2°K., respectively, 
for a reservoir temperature of 300°K. Thus, it is 
apparent that to eliminate condensation of air in the test 
section for hypersonic Mach Numbers, it is necessary to 
have a high reservoir temperature that is a function of 
the Mach Number and _ stagnation pressure. To 
understand the effects of impurities on the super- 
saturation of gases, as well as the aerodynamic proper- 
ties of two-phase air (liquid droplets and vapor phase), 
intensive experimental and analytical investigations 
have been conducted over the past several years, and 
the results are presented in references S—12. 

The study of air condensation is being conducted at 
considerably higher Mach Numbers in the Leg No. 2 
facility, which is designed for a maximum reservoir 
pressure of -1,000 Ibs. per sq.in. at 1,100°F. and a Mach 
Number range of 4 to 11. The impact pressure dis- 
tribution along the centerline of the nozzle was sur- 
veved by a seven-tube rake connected to a vacuum- 
referenced mercury manometer bank, and the static 
pressure distribution was measured by a Stokes gage. 
Since the static pressure along the axis is extremely 
low at high Mach Numbers, great care was taken to 
eliminate leakage and out-gasing of the static pressure 
leads. The stagnation pressure and temperature were 
measured by the Tate-Emery gage and a radiation- 
shielded thermocouple located just upstream from the 
throat. 

The static and total pressure measurements are 
presented in Fig. 4 as a plot of the ratio of the static 
pressure to the impact pressure, p/fo’, against the 
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ratio of impact to stagnation pressure, po’/po. Thus, 
when the condensation begins in the nozzle, the experi- 
mental points depart from the isentropic curve and 
approach the values predicted for a saturated expan- 
sion. With this method it is possible to determine 
the amount of supersaturation of the air expanded in a 
hypersonic nozzle. The summary of the results ob- 
tained so far on the amount of supersaturation, in the 
j- by 5-in. hypersonic wind tunnels, is as follows: 


Ibs. per To, AT. 
sq.in. “e. M, Moa _ 2 
165 266 7.0 5.6 33 
215 600 9.7 7.2 39 
415 800 10.5 7.9 38 
AT,. = supersaturation in termsof temperature (supercooling ). 
Msc = Mach Number with no supersaturation. 


These results seem to indicate that the amount of 
supersaturation of air is increased as a function of the 
Mach Number for the present operating range. 

In order to obtain more information regarding the 
state of the air during the condensation process, the 
temperature of the air was computed from the meas- 
The method used in the calculations 


ured pressures. 
By the integration of 


is discussed in reference 13. 
the one-dimensional continuity, momentum, and energy 
equations with condensation g, amount of air condensed, 
the temperature of the air can be determined. The 
results of the calculations are presented in Fig. 5, 
which indicates the amount of supersaturation that 
exists for the conditions that were tested. 


TRANSITION STUDIES ON AN INSULATED FLAT PLATE 


An extensive experimental study of boundary-layer 
transition was made on a 5- by 26-in. insulated flat 
plate in Leg No. 1 at a Mach Number of 5.8.14 This 
work was prompted by earlier experimental surveys 
indicating extensive laminar boundary layers contrary 
to a presupposed continual decrease in critical Rey- 
nolds Number with increasing Mach Number.! !¢ 

The experimental program was comprised of the 
following: transverse contamination due to the side- 
wall boundary layer; the effect of two-dimensional and 
local disturbances in the form of air jets of variable 
intensity from a row of orifices near the leading edge, 
and single air jets and small rods normal to the surface 
at various stations along the plate centerline, respec- 
tively; and the nature of the laminar and turbulent 
boundary layer. 

Full advantage was taken of the long uniform flow 
thombus at hypersonic Mach Numbers permitting the 
extensive plate length. The model was equipped with 
22 pressure orifices evenly distributed along the center- 
line, six copper-constantan thermocouples, and a span- 
wise distribution of small orifices 11/2 in. behind the 
leading edge for air injection into the boundary layer. 
Variable reservoir pressure gave a range of plate Reyn- 
olds Numbers from 1.4 X 10° to 5.5 X 10°. The pressure 
distribution was found to be fairly flat over most of the 
plate surface, as indicated in Fig. 6, except for a uni- 
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form positive gradient over the latter third due to 
leading-edge shock-wave reflection. The temperatures 
indicated an essentially-insulated surface. 

Two techniques were used for transition detection. 
One was phosphorescent lacquer, which dries rapidly 
in turbulent areas while remaining wet for long periods 
in laminar regions but which has the property of 
luminescing when dry under ultraviolet irradiation. 
The other was seven-probe total-head rake across 
half the plate span, geared for axial transverse ') '® in. 
above the plate surface. The vertical distance of the 
probes was chosen to yield a substantial increment in 
impact pressure between laminar and turbulent bound- 
ary layers which, over a large extent of the plate, 
were of the order of 0.2- to 0.5-in. thick. 

Transition by transverse contamination due to the 
mixing of the sidewall turbulent boundary layer with 
the laminar boundary layer on a flat plate was first 
noted about 15 years ago by Charters” in low-speed 
flow. 
spread from the vicinity of the plate leading edge ata 
uniform angle of 9'/.°, independent of speed, in the 
range of 20 to 60 m.p.h. A recent theoretical treatise 
by Parker" yields an angle of 9°16’ for the spread of 
jet in a quiescent fluid, assuming a linear velocity pro- 
file. 

The present research showed a marked delay in the 
it did not occur 


These experiments showed contamination to 


start of transverse contamination; 
until stations far downstream from the leading edge 
were reached. Subsequent surveys of the sidewall 
boundary layer at several stations above the flat 
surface revealed it as partially laminar at the lower 
reservoir pressures. However, in all cases transition 
to a turbulent boundary layer on the side wall was 
found to occur substantially far upstream from the 
origin of contamination on the plate, thus precluding 
an obvious explanation for the delay. In addition, 
at the higher reservoir pressures, the sidewall boundary 
layer was verified as already being fully turbulent near 
the plate leading edge. Contamination was found to 
originate at Reynolds Numbers between 1.5 10° and 
2 X 10° in all cases. The angle of spread was found 
to be 5'/2° and uniform over the plate length; in addi- 
tion, no natural transition was noted even up to a 
Reynolds Number of 5.5 X 10° in spite of the adverse 
pressure gradient over the latter part of the flat plate. 
These results were obtained by both the phosphorescent 
lacquer technique and the total-head rake, with good 
agreement between the two. 

A study was made of the effect of air injection into 
the laminar boundary layer as a transition tripping 
device. This type of disturbance was investigated by 
Fage and Sargent!’ as a means of fixing transition by 
adjusting the rate of injection and was recently used by 
Coles” in connection with direct measurements of 
turbulent skin friction in supersonic flow. 

Air was injected through single orifices along the 
plate centerline over a range of Reynolds Numbers 
from 10° to 2.8 & 10°. The phosphorescent lacquer 
revealed that the laminar boundary layer was not con- 
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taminated for a substantial distance downstream -of 
the jet. Pitot surveys in the wake directly behind the 
jet revealed local separation and verified the boundary 
layer as turbulent in the contaminated area. That 
separation occurs in the wake is not surprising in view 
of the extremely low kinetic energy of the laminar 
boundary layer close to the plate surface and the vertical 
momentum supplied by the air jet. There was always 
a distance between the jet and the start of contamina- 
tion, but the distance was noted to decrease as the 
Reynolds Number of the point of injection was in- 
creased. The corresponding Reynolds Number incre- 
ment at which transition occurred varied from approxi- 
mately 1.2 K 10° to 0.6 X 10° as the jet was moved 
from a Reynolds Number of 10° to 2.8 & 10° Small 
rods 0.25 in. long and 0.028 in. in diameter were 
mounted normal to the plate surface at a Reynolds 
Number of 10° and 1.6 X 10°, but they proved inef- 
fective in producing transition. The wake of the rods 
was traced downstream with no sign of contamination. 
Other surface protuberances and surface roughness 
were found equally unsuccessful. This evidence leads 
to the conclusion that air injection is one of the most 
effective transition tripping devices for high-speed flow, 
with the additional advantage of control, and that the 
laminar boundary layer is extremely stable at high 
Mach Numbers. 

Air at variable rates was injected through the row of 
orifices near the leading edge (Cf. Fig. 7). Transition 
was observed to move gradually upstream as the mass 
injected was increased to 2 per cent of the boundary 
layer mass defect on the plate. A further increase 
in the rate of injection did not appreciably move transi- 
tion to a Reynolds Number much below 2 X 10°. 

Boundary-layer surveys at a distance of 10 in. from 
the leading edge with a pitot probe showed the laminar 
layer to have almost a linear velocity profile as theoreti- 
cally predicted at high speeds*'** (Figs. 8 and Sa). 
The laminar boundary layer, however, was found to be 
substantially thicker than anticipated. A careful 
analysis of these results is being made at the present 
time. The turbulent boundary layer was character- 
ized by an almost linear impact pressure profile (Figs. 
9 and 9a) and a velocity distribution of the nature of a 
power-law profile. 

This study clearly illustrates the extensive stability 
of the laminar boundary layer at a Mach Number of 
5.8, which is probably the general case in extremely 
high-speed flow, and further analytical investigations 
are being conducted to extend the boundary-layer 
stability theory to higher Mach Numbers. 


INVESTIGATION OF SHOCK-WAVE BOUNDARY-LAYER 
INTERACTION AT HYPERSONIC SPEEDS 


In recent years, a considerable amount of informa- 
tion concerning the skin friction of a two-dimensional, 
thin flat plate moving through the air at high speeds 
has been obtained through the study of compressible 
boundary-layer equations.”* The usual basic assump- 
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tion is that viscous effects are confined to a narrow 
region close to the plate surface. The flow field oyt. 
side of this narrow region is regarded as nonviscoys 
and not affected by the presence of the boundary layer. 
However, when the plate is moved at supersonic speeds, 
disturbances due to the viscous effects in the boundary 
layer can drastically influence the outside flow field. 
Thus, the retarded flow in the boundary layer causes 
the stream line to turn slightly away from the plate 
surface, and this results in a shock wave extending 
from the plate into the free stream. There is, there- 
fore, a close relationship between the rate of growth of 
the boundary-layer thickness in the viscous flow regime 
and the rate of decay of shock-wave strength in the 
outside flow field. For the hypersonic speed range 
this interaction effect is much more pronounced, as 
discussed in references 1-4. Due to the strong inter- 
action, the pressure on the flat plate surface will not be 
constant but may be much higher than the ambient 
pressure, depending on the Mach Number and Reyn- 
olds Number. Correspondingly, the skin friction at 
hypersonic speed ranges will be higher than that pre- 
dicted by the conventional boundary-layer theory. 

The analysis of this hypersonic interaction problem 
through the integration of the boundary-layer equa- 
tions of motion is possible only if some simplifying 
assumptions are made regarding the boundary condi- 
tions. The lack of adequate experimental data has 
made it somewhat difficult to choose assumptions 
for the flow picture which will correspond to the actual 
physical flow phenomena. Analysis has been carried 
forward by several authors,'~* however, and has shown 
that perhaps the important similarity parameter in- 
volved is \/*/ \/ Re, the so-called hypersonic parameter. 
Another important parameter at hypersonic speed is 
the Knudsen Number, the ratio of mean-free path to 
the leading-edge radius. 

The experimental investigation of this hypersonic 
problem has been started with some tests designed to 
give reliable data on the flow conditions between the 
shock wave and the flat plate surface. A flat plate 
model was tested in the GALCIT Leg No. | tunnel ata 
Mach Number of 5.8 and Reynolds Numbers of 6 X 
104 to 25 X 104 per inch for one-phase flow conditions. 
One of these models was 7!/» in. long by 5 in. wide and 
had a 15° bevel forming a sharp leading edge of about 
0.001-in. radius. The chordwise pressure distribution 
on this plate was measured by means of nine static 
orifices located along the plate centerline from Q.2 in. to 
3 in. The orifice at 0.2 in. corresponded to values of 
\M*/v/ Re up to 0.8. These static pressures were read 
on a multiple tube manometer utilizing silicone oil with 
vacuum reference. The static pressure was close to the 
ambient pressure, which was of the order of a few cet- 
timeters of silicone oil and was read to 0.01 em. of 
silicone. The impact pressure distribution was surveyed 
between the flat plate surface and the shock wave at 
stations close to the leading edge of the plate in order 
that some information might be obtained on the be- 
havior of the air flow in the region of the interaction 
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between the boundary layer and the shock wave. The 
pressure was measured with a flattened probe, 0.016 by 
0.050 in., connected to a silicone fluid manometer to 
measure low pressures and toa mercury micromanometer 
to measure higher pressures. The visual observation of 
the flow was made possible by use of schlieren equipment, 
which was carefully aligned with the plate leading edge. 

The shock wave at the leading edge was curved 
since the flow was deflected by the growth of the 
boundary layer along the plate surface. Therefore, 
it is apparent that the static pressure along the plate 
surface would not be constant in that region. The 
corresponding plate surface pressure distribution is 
indicated in Fig. 10, where the data is presented in the 
form of the static pressure ratio versus 1/1/ Re. This 
data seems to indicate that, for this particular Mach 
Number, there is a nearly linear variation of the static 
pressure with 1/+/Re, which is in agreement with the 
original prediction.” 

The results of a survey of the impact pressure dis- 
tribution between the shock wave and the plate sur- 
face, taken at various distances from the leading edge, 
are shown in Fig. 11. This figure indicates that near 
the leading edge of the plate there exists a strong rota- 
tional flow region between the plate surface and the 
shock wave where the viscous effects should not be 
neglected. As the boundary layer grows along the 
plate, the viscous region close to the plate and the 
shock wave are separated by a rotational flow region 
with small viscous effect. At a greater distance from 
the leading edge the shock wave and the boundary 
layer are definitely separated by a great distance, and 
for the particular Mach Number that was tested the 
static pressure was approximately constant from 4 in. 
aft. Behind this region the interaction effect was 
negligible from the pressure distribution standpoint 
but was definitely still influenced greatly by the leading- 
edge shock wave. Tests are being conducted at the 
present time to obtain the corresponding static pres- 
sure distribution between the shock wave and the 
plate surface. This information, together with the 
impact pressure survey results, will make it possible 
to determine more accurately the actual flow condition 
between the shock wave and the plate surface. These 
results should be of aid in determining the proper flow 
model to be used in order to improve the analytical 


results. 
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Tests of a Fairing to Reduce the Drag of a 
Supersonic Swept-Wing Root 


W. F. HILTON* 
Sir W.G. Armstrong Whitworth Aircraft Led. 


SUMMARY 


The shock drag of a sweptback wing originates mainly from 
the wing roots, where the ‘‘S’’ shaped wing streamlines are con- 
strained to move along the side of the fuselage. 

An accurate method of calculating the free-stream shape of 
these streamlines was developed and used to design a model of a 
fairing for an EC.1240 wing swept back 65° from a flat plate. 
This fairing reduced considerably the wing root shock drag at .V/ 
= 1.73, as determined by pressure plotting. With the fairing re- 
moved, the flat plate-wing junction gave pressures and supersonic 
shock drag corresponding to a similar but unswept wing, thus 
demonstrating the importance of fairing the wing root. A wing- 
root vortex was observed in all cases without the fairing but was 
never observed with the fairing in place. 

The fairing was equally effective in reproducing subsonic fiow 
patterns at supersonic speeds over the whole incidence range 
tested (—4° to +4°). The tests were carried out in the Ordnance 
Acrophysics Laboratory Lone Star tunnel, Daingerfield, Tex., at a 
Reynolds Number of 4.24 million. Model strength at this high 
Reynolds Number prevented tests beyond 4° incidence. It is 
shown that the theory advanced in this report would also apply to 
a fairing between a swept-wing fuselage of any cross section, at 
any speed from J = 0 to MW = Mer. sec A. 


INTRODUCTION 


dee PRESSURES OVER THE FRONT PART of any wing 
combine to produce a drag force—a pressure or 
form drag. In subsonic flow the rear part of the wing 
contributes an equal and opposite thrust, due to pres- 
sure recovery, while supersonic flow is characterized by 
the rear part of the wing contributing an additional 
drag force, which adds to the fore drag to produce a 
large pressure drag. 

If a wing is swept back behind the Mach cone, it is 
theoretically possible to obtain the low drag subsonic 
flow pattern at supersonic speeds. The practical objec- 
tion to this scheme is that a supersonic flow pattern still 
exists over both wing root and wing tip and that struc- 
tural considerations prevent the use of long wing spans 
to reduce the relative importance of these high drag end 
effects. High-aspect-ratio sweptback wings also have 
bad low-speed characteristics and the low-speed sta- 
bility may have to be improved by the use of fences, 
drooped snoots, and other devices. 

This report describes an attempt to improve the flow 
around the root of a sweptback wing at supersonic 
speeds, in an effort to improve the aerodynamic effi- 
ciency of low-aspect-ratio wings. The fairing suggested 
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should also be as effective at sonic, subsonic, and low 
speeds. 

It is well known that a wing swept back at an angle 
A in a supersonic flow of velocity Ip (as shown in Fig. 1) 
will experience pressures from only the normal velocity 
component Vy cos A; the motion Ip sin A parallel with 
the wing will produce only skin friction, leaving the 
pressures undisturbed, apart from an extremely small 
correction due to spanwise growth of boundary layer. 

It follows that the velocity vector at any point on the 
wing is composed of the constant component I’) sin A 
parallel to the leading edge, compounded with a 
variable component normal to it. Thus the direction of 
the resultant air velocity over a swept wing will be 
variable. For example, at the leading edge stagnation 
point, where the normal component is reduced to zero, 
the air will move tangentially to the leading edge. The 
sinusoidal shape of a typical streamline is shown in Fig. 
1 for the case of a ‘‘subsonic’’ leading edge, where |’) cos 
A is less than the speed of sound, although the free- 
stream velocity I’) may in fact be supersonic. 

After studying the wavy path taken by a streamline 
in two-dimensional flow in Fig. 1, it isnot surprising that 
interference drag occurs when a swept wing is joined to 
a flat wall or fuselage without taking special precau- 
tions. If the wall were shaped to the two-dimensional 
streamline, then in the absence of boundary-layer ef- 
fects, we might reasonably hope to obtain two-dimen- 
sional flow right up to wing root with zero pressure 


drag. 


THEORY FOR CALCULATION OF FAIRING 


We start by calculating the low-speed pressure dis 
tribution round the wing section normal to the leading 
edge. Although many workers have given methods for 
calculating incompressible pressure and velocity dis- 
tributions on wings of arbitrary profile, few have even 
considered the problem of calculating the whole velocity 
field surrounding a two-dimensional wing. The method 
given by Mutterperl,* was adopted for use in calculating 
the pressure field round the wing section. This is as 
sumed to be correct for 1J = 0, and requires a com- 
pressibility correction. 

It has been found by experiment (pp. 89, 90, ref- 
erence 3) that the low-speed (incompressible) pressure 
coefficient C,; on a two-dimensional wing may be 
transformed to give the pressure coefficient C, at any 
subsonic Mach Number .\/, by using the von Karman 


relation 
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Fic. 1. Plan view of streamline over two-dimensional supersonic 


swept wing. 
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Fic. 6 Agreement of pressures with subsonic swept-wing theory 
ata = 0°. 
p at Po C 
- = Pp = 
(1/2) ypoM? 
C di 


; (1) 
V1 — M? + (1/2)C,, (1 — V1 — M?) 
where the first expression is an alternative definition of 
pressure coefficient. Thus, for example, in the all-im- 
portant region at the stagnation point C,; = 1, and the 
von Karman approximation reduces to 


: 1 | ; = 
Titian sa 9 + 9 Vi M? = 


1+ M? nm M* 5Mé6 4 . 
— — + —-+ete. (2) 

4 8 64 
as compared with the correct answer for pitot stagna- 
tion pressure coefficient given on p. 165 of reference 3. 


as 


= 1 + (M?/4) + (114/40) + 
(.17°/1600) + ete. (3) 


?stagnation 


Thus, for a resolved Mach Number of 0.73, von K4ar- 
man’s equation is less than 3 per cent in error at the 
forward stagnation point, and this error is even less on 
the upper and lower surfaces of the wing. Eq. (1) was 
used in computing the theoretical curves of C, shown in 
Figs. 3, 6, 7, and 8. 

The value of 7 in Eqs. (1) to (3) is the subsonic com- 
ponent 1/) cos A (normal to the leading edge) of the 
supersonic Mach Number J), so that a further correc- 
tion factor cos” A is needed if we are to refer the pressure 
coefficients to the free-stream conditions for a super- 
sonic aircraft having subsonic sweepback. In this latter 
case, only the subsonic normal component J/) cos A is 
active in producing pressure changes. Suppose that 
this normal component produces a pressure p, and a 
local Mach Number component /, at a certain point 
in the flow field. Then the isentropic relation between 
Static pressure and the Mach Number components 
producing them is 


Pn _ E + (1/2)(y — 1)M,? cos? I (y¥-) ” 
Po 1 + (1/2)(y — 1)M,? \ 


By rearrangement and defining a false C, in terms of 
My cos A so as to be consistent with Eq. (1), we can find 
M,, the local Mach Number component normal to the 
leading edge, in terms of known quantities 
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| =. (2 4 IT! ane cos? s) x 
ly — | 2 


-— 
, (: - rs yC,M,* cos? s) oe i! (9) 


Thus we can calculate ./, from a knowledge of wing 
geometry alone. The ratio of this calculated local 
Mach Number component normal to the leading edge to 
the constant spanwise .\/) sin A component gives the 
slope of the local flow direction, relative to the leading 
edge.* 

We thus have to synthesize a streamline from a 
knowledge of its slope at all points, which is not an ele- 
gant way of drawing a curve, but nevertheless practica- 


ble. 


APPLICATION OF THEORY TO DESIGN OF FAIRING 


We have seen above how to calculate the streamlines 
past a given swept wing. We now have the problem of 
combining these streamlines to form a surface (the fair- 
ing). 

The streamline shown in Fig. | is typical of a span- 
wise family of streamlines all parallel to one another. 
Above and below the wing are other different families of 
streamlines, also parallel in the spanwise direction. Any 
solid surface defined by choosing one member of each 
family would be a possible fairing. This choice could be 
made so as to define a circular or other shape of fuselage 
at one given station. The fuselage selected would 
strictly need to extend to infinity above and below the 
wing, but in practice the effect of cutting off the wing- 
root fairing at the top and bottom of the fuselage would 
be negligible. 

The effect of this approximation may be seen from 
Fig. 10, which shows the model for wind-tunnel tests 
of the above theory. It will be seen from this photo- 
graph that the shape resulting from these calculations 
is reasonably flat, except for a small region extremely 
close to the wing root and for a ridge extending to in- 
finity downstream from the trailing edge. 

However, when this type of fairing is used at super- 
sonic speeds, the beneficial effect of this ridge can only 
extend sideways and backwards and cannot improve 
the wing root itself. Therefore in a practical design it 
could be flattened out without serious loss of efficiency. 

An important part of the fairing is the tongue stretch- 
ing forward from the leading edge junction with the 
body. This tongue sends out wavelets ahead of the 
wing which start to deviate the air to either side of the 
leading edge before the air reaches the wing, as in sub- 
sonic flow. This avoids the necessity for a pair of in- 
clined shocks to deviate the streamlines discontinu- 
ously, as when the fairing is absent. This physical ex- 
planation is discussed at greater length on pp. 303-305 
of reference 3. It is important to note that the shape of 


* This assumes the velocity of sound is constant throughout the 
flow field. This is not quite true for extremely large velocity 
variations, and, strictly speaking, velocities and not Mach Num- 
bers should be compounded. See pp. 505 and 508 of reference 3. 
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TESTS OF 
this fairing is independent of design Mach Number, 
apart from the small von Karman correction of Eq. (1). 
Thus the fairing is effective at supersonic, transonic, 
subsonic, and even low speeds, where it could at least 
avoid the loss of wing circulation in trailing vortex 
motion. 

The “S” shape of the 
in the swept-wing geometry and is only slightly modi- 
Another important 


streamlines is inherent 
fied because of compressibility. 
feature of this fairing is that it is cut away at the rear of 
the wing, and this cutaway persists in that part of the 


fairing behind the wing root. 


CHOICE OF MODEL FOR TEST 


Any Mach Number, subsonic, supersonic, or even 
M = 1 could have been chosen for testing the theory 
without loss of accuracy, but J = 1.73 was selected as 
being sufficiently high as to preclude the possibility of 
obtaining subsonic flow patterns by accident. A sweep- 
back angle of 65° was chosen so as to give a normal 
Mach Number component 1, = 1.73 cos 65° = 0.75, 
which is below the critical speed for this section at the 
design incidence of 0°. A symmetrical wing was chosen 
because it halved the computation of the shape of the 
fairing and halved the number of pressure holes re- 
quired on the wind-tunnel model. A wing cambered to 
produce the appropriate C, would be used in a practical 
aircraft design. 

The shape of the fairing was calculated by Miss B. F. 
Fuess,' but the details of the calculation have not been 
reproduced in this report since they are somewhat 


lengthy. 


DETAILS OF EXPERIMENT 


Some 3,600 pressure measurements and several flow 
photographs were made in the Ordnance Aerophysics 
Laboratory 27'/s- by 19-in. wind tunnel at a Mach 
A sketch of the model is given in Fig. 
The wing chord in the 


Number of 1.73. 
2, and a photograph in Fig. 10. 
direction of motion was 8.3 in., giving a Reynolds Num- 
ber of 4.24 million. 

Pressures were read on four rows of pressure holes at 
four different spanwise stations; wing incidence was 
obtained by rotating the whole model about a vertical 
axis, enabling experiments to be made at a = —4°, 
—2°, —1°, 0°, 1°, 2°, and 4°. The flat plate was sup- 
ported parallel with the floor of the tunnel as shown in 
Fig. 2, and, in addition, two runs were made with the 
plate tilted at 2° of yaw (Wy = 2°), so as to increase the 
Mach Number on the wing. This also represents a 
typical off-design condition. 

The model was made of steel, and the flat plate was 
held out in clean air away from the tunnel wall by a 
sweptback strut, which also carried the pressure con- 
nections. The fairing was made as a split aluminum 
collar which could be dropped over the wing or removed 
to simulate the unfaired condition. Some repeat read- 
ings are plotted in Fig. 3 from which it will be seen that 
no systematic variation occurred between repeat runs 


A 
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and that the random errors are of the order of 0.02 in 
pressure coefficient. The consistent irregularities ob- 
served in these pressure curves at certain pressure sta- 
tions presumably arise from similar irregularities in the 
static pressure distribution in the tunnel, or from ir- 
regularities in the manufacture of the model. 

A more serious form of disturbance was caused by 
the shock wave from the leading edge of the support 
underneath the flat plate. This shock spread out under- 
neath the plate until it reached the edges, and then it 
turned round the edges and propagated on the upper 
side of the plate until it hit the model. The discon- 
tinuous rise of static pressure in the shock wave upset 
all the pressures observed at holes situated behind it in 
Row III. 

It must not be overlooked that, by eliminating pres- 
sure drag on the wings, we have in fact incurred a pres- 
sure drag on the fairing. It has been shown by Dr. C. H. 
McKoen that this effect is of negligibly small amount. 
The drag coefficient for a similar wing-root fairing was 
calculated to be 3.7 by 10~*, based on unit wing area. 

The publication of this report was delayed in the 
hope of including further tests at Mach Numbers less 
than 1.73, so as to investigate experimentally the effect 
of off-design conditions of speed as well as incidence, 
but these have not been done. 

However, it must be stressed that the whole fairing 
would be similar in shape if designed for 17 = 0 (in- 
compressible case) and would make a good wing-root 
fairing to avoid the formation of a trailing vortex from 
the wing root. This vortex, seen in Fig. 9, would give 
rise to loss of lift and increase of drag even at low in- 
compressible speeds, when it would be described as 


wing-root interference drag. 


ANALYSIS OF PHOTOGRAPHS 


The flow photographs taken by parallel light, direct 
shadow, spark photography of the type shown in Fig. 9 
consistently showed two effects when the special fairing 
was removed—namely : 

(1) A separation of flow at the junction of the swept 
wing with the unfaired flat plate. This separation per- 
sisted the whole way back along the flat plate and was 
about 0.03 wing chords in thickness. 

(2) One trace and stometimes two traces are also 
visible on these photographs, and these are of a similar 
nature to the airscrew tip vortex photographed in 1935.4 
It is fairly certain that these traces represent vortex 
motion, thus indicating that the lift and the circulation 
round the wing has changed rapidly in value in a span- 
wise direction. 

Hence the junction of this swept wing with a flat 
plate caused a loss of lift. In truly two-dimensional 
flow no such loss of lift would be observed. 

While all of the photographs taken without the special 
fairing showed both vortex and separation, none of the 
photographs with the fairing showed either effect. 

This is probably the most convincing evidence that 
such a design of wing-root fairing is effective in cleaning 
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Photograph showing shape of fairing. 
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up the bad flow generated by the simple geometrical 
joining of a swept wing with a flat plate. 


PRESSURE MEASUREMENTS WITH FAIRING 


The upper part of Fig. 3 shows the pressure distribu- 
tion measured in the design condition a = 0° and y = 
0° superimposed on this theoretical curve. It will be 
seen that all pressure observations fall close to the 
theoretical curve, although a tendency exists for the 
pressure near the wing root to fall below the theory by 
about 0.05 on pressure coefficient. The pressures fur- 
ther out on the wing span are slightly closer to the 
theoretical value, as will be seen from Fig. 6. On the 
whole, it may be said that the pressures on the 65° swept 
wing plus specially shaped wing-root fairing agree well 
with two-dimensional swept-wing theory and, hence, 
that the streamline fairing proposed for the root of such 
a wing gives a good approximation to the pressures cor- 
responding to two-dimensional flow. 


PRESSURE MEASUREMENTS WITHOUT FAIRING 


The wing-root pressures measured in the absence of 
the fairing are shown in the lower part of Fig. 3. These 
are quite different from the pressures observed with the 
fairing in place, the positive slope of the curve indicat- 
ing that the air was accelerating all the way back, as in 
supersonic flow over an unyawed wing. This was un- 
doubtedly caused by the flat plate constraining the air 
to move in a two-dimensional supersonic pattern over 
the wing root and is the effect that the fairing tries to 
compensate. 

In order to demonstrate this point more effectively, 
the supersonic pressure distribution was calculated for 
an unyawed EC.1240 wing thinned down in the ratio 
cos A to allow for the reduction in apparent thickness- 
chord ratio due to sweepback. 

Busemann’s second-order approximation 


Cy _ Cid + Cop” (6) 


was used in these calculations, ¢ being the local angle 
of the unyawed wing surface to the wind. Since this 
was a subsonic wing section, the calculations were not 
taken too close to the rounded leading edge. The 
answers are thought to be accurate enough for our pur- 
pose. 

It will be seen from Fig. 3 that this approximate super- 
sonic theory agrees well with the observed pressures on 
the unfaired wing root. 
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EFFECT OF VARIATIONS FROM DESIGN LIFT COEFFICIENT 


The fairing tested was designed to work at zero lift, 
and it is of interest to study the loss of efficiency at off- 
design conditions. The comparison of theory with ex- 
periment for both faired and unfaired wings at an angle 
of incidence is shown in Figs. 7 and 8. 

No undesirable effect due to the lift is apparent from 
the results for the faired wing root. The results for the 
lower surface at a = 2° and 4° were obtained by using 
the upper surface results at a = —2° and —4° and by 
assuming that the wind direction was in line with the 
chord at a = 0°. The unfaired wing root produced the 
same undesirable flow separation, vortex motion, and 
pressure drag noted previously. 

Thus the fairing might be expected to improve the 
flow for quite a range of lift coefficients. 


REDUCTION OF PRESSURE READINGS TO GIVE PRESSURE 
DRAG 


So far the pressure readings have been plotted against 
chordwise distance in what might be termed the con- 
ventional manner. 

However, one of the characteristic differences be- 
tween wing pressure distributions at subsonic and 
supersonic speeds is that the pressures make zero con- 
tribution to the drag at subsonic speeds, while a super- 
sonic pressure distribution adds up to give a large pres- 
sure drag, which is sometimes termed ‘‘wave-drag.”’ 
Tangential (skin friction) drag is not detected by this 
method of observation. 

Pressure readings are reduced to pressure drag by 
plotting them against thickness instead of in the chord- 
wise direction, as indicated at the top of Fig. 4. Many 
of the sets of pressure readings were so plotted, in order 
to distinguish between the subsonic pressure distribu- 
tion (zero-pressure drag) and the supersonic distribution 
(high-pressure drag), which we are seeking to avoid by 
the use of this fairing. Therefore, it was felt that this 
would be the most convincing method of presenting 
the data to show the performance of this fairing. Low- 
speed theory shows that the lower central loop indicat- 
ing nose drag is balanced by the curves crossing to give 
two negative “‘ears’’ representing thrust due to pressure 
recovery on the rear parts of upper and lower surfaces, 
respectively. 

Unfortunately there were not sufficient pressure holes 
near the leading edge or the trailing edge of the wing to 
permit exact evaluation of the pressure drag coefficient 
by this means. However, several sets of pressure data 
have been plotted in this manner in Figs. 4 and 5, in 
order to demonstrate the great reduction in pressure 
drag due to adding the fairing. 

Fig. 4 shows conditions close to the intersection of 
wing and body in Row II.* In the unfaired case, every 
part of the wing having pressure holes contributed 


* Row I, which was closest to the wing root, could not be used 
for this method of plotting, because the first three holes at the 
leading edge were out of order, and these holes are essential to an 
evaluation of the pressure drag. 

(Continued on page 188) 
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The Effect of Thickness on Airfoils with 
Constant Vertical Acceleration at Supersonic 
Speeds 


JOHN C. MARTIN* ann NATHAN GERBER? 


Ballistics Research Laboratories, Aberdeen Proving Ground 


ABSTRACT 


The effects of thickness on the lift and pitching moment of 
two-dimensional airfoils at supersonic speeds with constant ver- 
tical acceleration are investigated. The airfoils considered have 
arbitrary symmetrical cross sections. and the flow is supersonic 


everywhere 

The analysis is based on a second-order theory similar to the 
second-order theory introduced by Busemann and extended by 
Van Dyke. The lifting pressure due to a constant vertical ac- 
celeration is found, and this is used to calculate the lift coefficient, 
Ca? and the moment coefficient, Cg due to a constant vertical 
acceleration. 

The airfoil’s second-order contribution to the damping of 
longitudinal oscillations in aircraft is considered, and its relation 
to the damping of oscillating airfoils is investigated 


(1) INTRODUCTION 


7 DEVELOPMENT OF THE linearized theory of super- 
sonic flow has permitted a first-order evaluation 
of a number of stability derivatives. Second-order 
theories similar to the one introduced by Busemann! 
and extended by Van Dyke? * offer possibilities of ob- 
taining second-order evaluations of certain stability 
derivatives, such as the lift, Ci.» and the moment, 
Cn» due to a constant vertical acceleration. Examples 
of the use of second-order theories to evaluate stability 
derivatives are found in references 4 and 5. 

Recently, two papers have appeared® * dealing with 
two-dimensional oscillating airfoils. The flow over 
oscillating airfoils is of great importance in flutter cal- 
culations. For stability studies, however, it is more 
convenient to know the flow over airfoils with a con- 
stant vertical acceleration and a constant rate of pitch. 

In this paper a second-order theory is developed for 
two-dimensional airfoils with constant vertical acceler- 
ation at supersonic speeds. From this theory the 
lifting pressure due to acceleration is evaluated. The 
expression for the lifting pressure enables the stability 
derivatives Ch, and Cm to be calculated. The airfoils 
considered here have arbitrary symmetrical cross 
sections; however, the analysis can easily be extended 
to include airfoils with unsymmetrical cross sections. 

The flow around an accelerating airfoil is of an un- 
Steady nature and necessitates the use of time-depend- 
ent partial differential equations. The differential 
equation studied here is not the equation commonly 
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used, since the coordinate axes employed in the an- 
alysis are attached to the airfoil. 

Recent work by Van Dyke® * indicates that second- 
order solutions of the partial differential equation of 
steady supersonic flow can be obtained by iterative 
methods. In this paper it will be assumed that iter- 
ative methods can also be used to obtain second-order 
solutions for unsteady supersonic flow. 

In addition, it will be assumed that the characteristics 
are the same for the first- and second-order solutions. 
This assumption does not appear unreasonable, since 
for steady plane flow the second-order solution’ found 
by an iterative method yields the correct pressure of 
the Busemann second-order theory. 


(2) THe PARTIAL DIFFERENTIAL EQUATION 


The partial differential equation for the second-order 


potential flow around an airfoil is* 


— BP, + BO — 


l 
Ve J vy — 1) (0.4 1) (D..7 + B,,7) + 


9 9 

26D + 2yDy + 7 P,P + 7 Py PoyP? (1) 
where x’, y’, 2’ are rectangular coordinates with x’ 
parallel to the free-stream direction, ¢t’ represents time, 
V is the free-stream velocity, c is the velocity of sound, 
M is the Mach Number, 8 = V MM? — 1, and y is the 
adiabatic exponent. The perturbation potential func- 
tion, &, has been normalized through division by the 
free-stream velocity. Also from reference 8 the pres- 
sure coefficient is given by 


pressure 





Cp 


] 
—~=—2 (+. -* 4) ~ 0, + B*4,,? + 
(1/2) pV? ] 


where p is the density. 

These two equations are associated with axes fixed 
in space. For this problem it is convenient to express 
the differential equation in terms of a new set of axes 
fixed to the airfoil. The relations between the two sets 
of axes are (see Fig. 1) 
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~ = x" 
z2=2' + (1/2)aVt? (3) 
i = jf’ 


where 4 is the acceleration. 
For the axes attached to the airfoil, Eqs. (1) and (2) 


become 
2M? 1 
— $’°,, + 2, — 7 ?,, — eo, = 
c? 
J at a (y — 1) 
2 M* 7 al®,. = = V ®., = 9 V ®, + . 9 x 


l 
(o. + <_ (b,, + P22) + 0, + P,P,2 + 


a : ne.) (4) 
py PePre F 7, P2Pery 


I 
C,= -2 E +5 + aa, | — 2 + B,? + 


2M? M? 
$,P, + 


@? (5) 
l 12 


(3) SOLUTION BY ITERATION 


On the assumption that Eq. (4) can by solved by an 
iteration procedure, the initial step will be finding the 
first-order solution. The first-order partial differential 
equation, obtained by neglecting the second-order 
terms in Eq. (4), is 

2M? l ; 
— Borer + 2 — V Ort — a ou = 0 (6) 

The solution to Eq. (6) is taken as the first approxi- 
mation to the solution of Eq. (4). We shall assume 
that the second approximation can be found by sub- 
stituting the first-order solution into the right side of 
Eq. (4) and solving the resulting nonhomogeneous 
equation, which, from Eqs. (4) and (6), is 


2M? 
= B'Wrr a W22 -> V Wet >= Ce Vin = 
ee F at a Bete. 
2M" atd,: -+- V Pet + ) yp? = 9 x 


l 
(¢. =) 1%) (Grr = a $::2) ar Prrr = PrPrz = 


l | | e 
pr P2P ee ia V PP (4) 


The solution of Eq. (7) will be referred to as the second- 


order solution. 


(4) EXAMINATION OF FORM OF POTENTIAL FUNCTION 


It is helpful to investigate the type of solution ob- 


tained from Eq. (7). The first-order solution will be 


of the form 
> = a gi(x, 2) + & gol(x, 2, t) + € g3(x, 2) 


where g;, go, g3 are functions of the variables x, z, and ¢, 


RONAUTICAL SCIENCES 


MARCH, 1955 


a is the angle of attack, € is a thickness parameter, and 
a, a, and e are small compared with unity. It follows 
from the above expression and Eq. (7) that the second- 


order solution will be of the form 


Y = ah(x, 2) + adhe(x, 2, t) + a&7h3(x, 2, t) + 
acha(x, 2) + aehs(x, 2, t) + Che(x, z 


where fy, to, .. . , 4g are functions of the variables x, z, 


and ¢. 

This paper is concerned only with the lifting pressure 
and its integrated effects. The terms a(x, 2), 
adh»(x, 2, t), and a@7h3(x, 2, t) contribute nothing to the 
lifting pressure. This can be seen from the following 
argument. The thickness parameter, ¢, is not present; 
therefore, the airfoils can be considered to have zero 
thickness when these terms alone are treated. If the 
potential of the flow on the upper surface of a flat plate 


is expressed as 
P(a, a) = agi + age + ah, + adh, + ah; 


the potential on the lower surface is given by 
@(—a, —a) = —ag; — & £2 + ah; + adh, + dehy; 
The potential difference is 


Aé = 2ag) 4 2 age 


Since for the flat plate the pressure difference between 
the upper and lower surfaces can be found directly from 
the potential difference, the terms a*/, adhe, and a*h; 
do not contribute to the lifting pressure. 

The terms on the right side of Eq. (7) which are 
multiplied by a®, ad, and & can be neglected without 
changing the expression for the lifting pressure, for 
Eq. (7) is linear, and a*h;, adh2, and &*h3 do not con- 
tribute to the lifting pressure. The terms ae/y and 
eh, are independent of @ (the acceleration parameter), 
and hence they also contribute nothing to the pressure 
due to acceleration.* 

The remaining term de/;(x, 2, ¢) will be found by use 
of Eq. (7), neglecting all expressions involving a’, ae, €, 
a’, and ad. 

A further consideration of the form of the solution 
indicates that the second-order lifting pressure is linear 
in the thickness parameter. This can be established 
by considering a first-order solution of the form 


eo) 


@ = ag,(x, 2) + adge(x, 2, ft) + erga(x, 2), €285(X, 
é d s 


It follows that the second-order lifting pressure will be 


of the form 
AC, —_ Io(x, t) + ae 1) (x, t) + deol o(x, t) + 


This equation is linear in « and e; thus the lifting 
pressure for various known thickness distributions cam 
be added to obtain the lifting pressure for new thickness 
distributions. 

This linearity has no special value for the two-dimen- 
sional airfoil, since we shall determine the solution for 

* The terms aehy, €2h¢, and a2h; are associated with steady 
supersonic flow and can be found by the Busemann second-order 


theory. 
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an arbitrary thickness distribution. It can be shown, 
similarly, that linearity of thickness distribution holds 
for three-dimensional airfoils also, and for this case it 
should prove quite useful. 


(5) BOUNDARY CONDITIONS 


Physical considerations require that the flow be 
tangent to the surface of the airfoil and that all ve- 
locity perturbations vanish upstream of the airfoil. 
These boundary conditions may be expressed mathe- 
matically as 

o(x, y, z, t) = 0 
y(x, y, 2, ¢) = O 


upstream of the airfoil and 


where v is the velocity vector and s (x, 3) = 0 is the 
equation of the surface of the airfoil. 
The equation of the surface of the airfoil may also 
be expressed as 
z= f(x) 
Thus 
Vs = —ief’ +k 
where 7 and & are unit vectors in the x and gz directions, 
respectively. Since the velocity vector, v, can be 
written as 
v=iul+o¢,+ y,) + k(@. + v. + at) 
it follows that the boundary condition on the airfoil 
surface is given by 
-(i+¢+ vy)’ +¢.+ ¥.+ at = 0 


The coordinate axes will be chosen so that the airfoil 
lies approximately in the z = 0 plane. The boundary 


conditions for the first-order solution are given by 
o(x, y, 3) = O upstream of the airfoil, and 


$2\:=0 = —a — at + ef’ (S) 


Since only the de term is being calculated, a will be set 
equal to zero. 

Similarly, the boundary conditions for the second- 
order solution are given by Y = 0 upstream of the air- 


foil, and 
Velz-0 = €hriz-0f — fo ( (9) 


For the airfoils considered in this paper the first- 
order velocity components are discontinuous across 
the Mach sheet from the leading edge. The effect of 
these discontinuities on the second-order solution will 
be evaluated by assuming that a small transition sur- 
face is attached to the leading edge. This small surface 
is so shaped that the discontinuities in the velocity com- 
ponents are removed. The effect of the discontinuities 
will be evaluated by a limiting process in which the 
width of the attached surface approaches zero. A sur- 
face that will remove the discontinuities in the first- 
order velocity components in the flow over the upper 


surface is 


z= (1/2) (—oxt + xa) 
where the surface extends from x = 0 to x = / and 
ol = &, al =e 


Evaluating the second-order potential function asso- 
ciated with the flow over the small transition surface 
introduced above will lead to the determination of the 
second-order leading-edge discontinuities. The first- 
order solution for the flow over the transition surface is 

(x — Bz)? o M(x + 262 : ; | 
6 = - - + g¢lt — Va 
261 3B" 
(10) 


The second-order solution can be found by substituting Eq. (10) into Eq. (7) and solving the resulting nonhomo- 


geneous equation; however, in this case it is easier to use a different approach. 


In Appendix (A) it is shown that 


the solution for planar problems to the three-dimensional nonhomogeneous equation 


— B’y,, -- Woy . V2 = (217? V) Wee — (1 Wit = Fi X,Y, 3, t) (11 
can be written as 
((f —— if 1 Pov(é, 7, 0,4) . Olé, n, 0, 2 | 
(x ¥, 3,4) = — F(t;) + F(ts)] dv — + la (12) 
v > fare JdJit R | bal | ; Dr Jd SS R or O& a 


where F is an arbitrary function; é, n, ¢ are auxiliary variables for x, y, and 2; 


S 


(x — &)M R (x — &)M R 
h=t— — —-—, h=t— = + — 
B°c BC B°c B-c 

R = \ zs — &)* — B?(y — n)* — Bs _ ¢)? 


v% is the volume inclosed by the forward Mach cone from the point (x, y, 2); and the surface S, is the area of the 


transition surface in the forward Mach cone from the point (x, y, 2). 


Z) 


When Eq. (12) is applied to an arbitrary point located on the Mach surface from the trailing edge of the transi- 


tion surface, the result can be expressed as 
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2 


1 2 Bo+l /a)Va—?-B@=H? | 
¥(x, 2, ¢) = f ac f ae f [F(t,) + F(te) |dn — 
4m Jo Br -(1/82) Va - 8? —we- 9 R 


)2 
| l (1/8) V(x — 8? — Bs? | (E, 0, 0, ¢ (E, 0, O, to 
i) ae | |e ae 4 ee 
2r Jo —(1/8) V(x — 8? — Ba R or o¢ 


Since F, (Ow /O¢) (&, 0, 0, t:), and (Ow /O¢) (é, 0, 0, f2) are independent of n, the preceding equation can be reduced to 


1 4 Be +l : : 1 1 ra) Pa 
¥(x, z,t) = f ac f F(é, ¢, th) + FE, §, t2) dé — - Y (é, 0, O, t:) + v (£, 0, 0, 0, f2) | dé (13 
4J0 Bt 270 LO¢ rer 


The integrand of the first integral in Eq. (13) can be written as 


F(é, cs ty) + F(é, c, ty) ial Jo(&, &; t) + (€ tF BO) SLE, - t) + (E ra, BC)? So(E, e. t) + (€ = BC) FT3(&, é. t) 


where Jo, J;, J2, and J; are continuous functions and are not zero at & = BF. 

Since the é in the integrands of Eq. (13) varies over a small range, the functions Jp, J;, J2, and J3 can be con- 
sidered as constants with respect to integration in the & direction. Therefore & can be replaced by 6¢ in these 
integrands, thus permitting integration to be performed with respect to —& Eq. (13) becomes 


» 


] V4 eZ 
W(x, z, t) = { Jy(pe, ¢, de + { Ji(B¢, ¢, de + 


13 PY 0 SB « 


a ne l on EPeMt (x  1\ éVt , a 
128 Jo J (Be, c t)dé + 168 Jo J3(B¢, ¢, =) di + > 3? > — 3 — 9 + 9 (14) 


The value of the second-order potential function on the downstream side of the leading-edge Mach sheet above 
the airfoil can be found by taking the limit of Eq. (14) as / approaches zero and as ¢ and a approach infinity so 
that 

lim o/ = & 
(1 > 0) 


(o¢ —~> ~) 


and 


lim a/ = e« 
(1 > 0) 


(a—> ~) 


Thus the discontinuity in the second-order potential function is given by 


¥(x = Bz) = lim a J)(Bé, ¢, Ede + lim 5 f Ji (Be, &, tdt (15) 
(1 —> 0) (1 > 0) 
(¢ > ~) (¢ —> ~) 
(a—> @) (a—> @) 


The integrals containing Jz and J; are zero in the limit, since J. and /3 are linear functions of ¢?, ca, a’, a, and 


ad. 
From Eqs. (10) and (7), the functions J) and J, are (neglecting the o’, ac, a*, and &? terms) 


j= - falas —t + [M(x — Be) ‘B?c]} 
J, = —4A(y + 1)Mtca[Vt — (M*x/B?) ]/( VB?) 


Substituting the preceding expression into Eq. (15) and performing the indicated operations yield 


M? Mx +1 MM Mx 
| : 7 ioe " >) 
Vix =p: = de 3? | - + a | Bz + de 5 3! [ — a | Bz (16) 


Eq. (16) is the value of the discontinuity in the second-order potential function across the leading-edge Mach sheet 


above the airfoil. 
(6) SOLUTION OF THE PARTIAL DIFFERENTIAL EQUATION 


The part of the second-order potential function which contributes to the lifting pressure because of a constant 
acceleration will now be determined. Since the method of obtaining the potential is essentially the same for both 
the upper and lower surfaces, only the flow over the upper surface will be considered in detail. 
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The first-order solution is 


A i) Ec) f(x — Bz 7 
9 = ae . 22+ Vt(x — Bs) | — €° (14 
BI a" 2 8 


’ 
It follows from Eqs. (7) and (17) that the second-order potential function must satisfy the nonhomogeneous 
equation 
f’ — 


= B'Yrr = ie V22 _- > Vr seiealiie™ Vir - 34 


2M? l Mae k + y1I/7(1 — B 
| c* V 


I : - 
go (VM? — MP + 2) Vif + 2. (y+ 1) wt BI (18) 


where the a’, ad, ae, €’, and &* terms have been neglected. 

The potential function y will be divided into two parts 6; and #2 such that 6; satisfies the nonhomogeneous equa 
tion and 6, satisfies the homogeneous equation. 

By inspection a solution of the nonhomogeneous equation is found to be 


ae? ae Al? 


i= — (M?y — M? + 2)ztf’ — y(81 — 2) + M — 3)ef — 
l 233 Y if 48° [y 3 |; 
ae? adel? 
Bt — y(M+ — 1) ]2?f’ + — (y + l)wsf’ (19) 
iL | Y |s* 27° Y . 
The boundary conditions are given by Eqs. (9) and (16). From Eq. (19) and the boundary conditions it fol 
lows that 6.|, -g: = 0 and 
06» ve Ma? ig ae . 
: eee | -ary — M*+ 2) Vif’ — (—2M*y + 464 — 2)f + — (M*y — AM? + 2)xf | 


2B? . 5” 


Oz |, =0 2VB* 
By inspection 


ae jf 2 M? 
§. = <(Mty — M+ + 2)Vtf+ 


5V6 ) — [— MP*(y — 1)6?2 — (M*y — M? + 2x) f + 


Mw y ithe \ 

(— 2Mty + 6M* — 10M? + 2) f(ujdu (20) 
28? Jo 

It follows from Eqs. (19) and (20) that the part of the second-order potential function which contributes to the 
lifting pressure is 

a en a —_ ae = i4 oy 
y= —— ~(M*(M*y — M? + 2)62Vif’ — yM'*p?2*f’ — (y + 1)xBef’ + (M*y — M* + 2)Vb¢f - 
2VB4 \ 8? 

M 


2 — x Bz 
[((Mty — M* + M? + 1)B2 + (M*y — M? + 2)x/f a (Mty — 3M* + 5M? — vf f(u)dut (21) 
B? / 0 
(7) LIFTING PRESSURE 
The lifting pressure distribution can be expressed as 


AC "ih Sala C p upper = Cy lower 


surface surface 


From Eqs. (2) and (21), it follows that for symmetrical airfoils the lifting pressure distribution at ¢ = 0 is 


: aAx 2aeM? ° " outa F ‘ — si . oe 
AC, = 793 + 79 (Mey — M? + 2)xf’(x) + (2yM? — 3M? + 5) f(x)] (22) 


(S) THE Force AND MOMENT 


The effect of thickness on the lift, C, ., and the moment, C,, ., due to a constant vertical acceleration can be found 
a a 


by use of Eq. (22). The stability derivatives CL. and Cm, can be expressed as (c, = chord length) 
: normal force 
CL, =0 = a 
is (1/2) pV°c, 21 “ ‘ 
: = -" AC, dx 
O(ac,/2V) —" ac,” Jo 


and 
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moment 
. (1/2) pV°c,? Ay 
C., = — —= = — (x — d)AC, dx 
- O0(ac,/2V) a~o &¢,? Jo , 


From the preceding relations and Eq. (22) the C;, ‘ and the C,, ; of an airfoil are found to be 
. oF | eM? 
c B3 jo 283 


. S J l d el? ? ’ d fy E 
c m = ’ — M? — M? y l = ‘ 
«p83 26, 7 2B [ _ ( 7 c, 


l ioe "er sef( sx) 
. (M*y — 2M* + 3) | = dx — 3° | ” ax |' (24) 


C; 0 ¢," Ye 4 


Cr | (ry — M? + 2) Ire + (Mey — 2M? + 3) is) ax | (23) 


Cr 0 ¢,* 


The sum of the stability derivatives Cm; and the damping in pitch, C,,, provides a large contribution to the 

qd 
damping of the longitudinal oscillations in aircraft. The second-order damping in pitch for airfoils with sym- 
metrical thickness distributions is given by Eq. (23) of reference 5. The sum of this equation and Eq. (24) can 


be written as* 


: . S| (1 — B*) (28° — 1) fd a\? 
“4 m i Cc m = — 
= . | 3p" ” 26? i [ | 


; iF (—M* + 4M* — 6M? + 4+ yM*) + : (M?y + M* — 4M? + »|(a _ a) ize -- 
i Zz c » 


: G4 
d l ‘ ; - "er F(x) 
) = [yW4(2\77 — 1) + 247% — 11474 + ISA? — 8] | x 
Cr §* / 0 C,* 
Pore ee “ vf(x) |) 7 
dx +> 2(M v + Mt — 44/° + 4) = dx ¢ (25) 
J 0 vee 
* The integral in the C,,, expression [Eq. (23), reference 5] was changed by letting & = x — d, where d indicates the position 
of the axis of pitch. This put it into the same form as the integrals in the expression for C,, 

Fig. 2 presents the variation of the C,, a Cy of a the upper and lower surfaces of the airfoils treated are 


independent of each other, the aerodynamic properties 


wedge airfoil with Mach Number for various thick- 
due to constant vertical accelerations of airfoils with 


nesses for the case where the axis of pitch is located at 

the mid-chord point. Fig. 3 presents the variation of 

the C,. + Cm, of an airfoil with a parabolic cross sec- 
a 4G 


unsymmetrical thickness distributions can easily be 
determined from the results obtained here. 
The limitations of the Busemann second-order theory 


tion with Mach Number for various thicknesses for the 7 
have been investigated (see reference 9). Since the 


case where the axis of pitch is located at the mid-chord 
theory contained in the present paper is closely asso- 


ciated with the Busemann second-order theory, it 
seems likely that the results presented here have similar 


point. These two figures indicate that thickness has a 
destabilizing effect for a wedge airfoil and has a sta- 
bilizing effect for an airfoil with a parabolic cross section. nahin gh 
Fig. 4 presents the regions of possible instability for limitations. 
a 10 per cent wedge and a 5 per cent parabolic arc 
airfoil. The curves are lines of zero damping found by tanta 
placing Cm + Cm, equal to zero. This figure indicates 1 Busemann, A., Aerodynamischer Auftrieb bei Uberschall 
geschwindigkeit, Atti die convegni 5, R. Accad. d'Italia, pp. 328- 


that the effect of thickness increases the region of in- 
360 1936, (also printed in Luftfahrtforschung, Bd. 12, Nr. 6, 


stability for the wedge airfoil and shifts the region to a je : : ae 
: -* : see be ™ pp. 210-229, October 3, 1935; available in translation as British 

lower value of d/c, for the parabolic airfoil. Air Ministry Translation No. 2844.) 
The relation between Cn, + Cm, and the damping 2 Van Dyke, M. D., A Study of Second Order Theory of Super 


sonic Flow Past Bodies of Revolution, NACA TN No. 2200, 1951 
3 Van Dyke, M. D., First- and Second-Order Theory of Super 

sonic Flow Past Bodies of Revolution, Journal of the Aeronautical 

Sciences, Vol. 18, No. 3, pp. 161-179, March, 1951. 

: + 4 Martin, John C., and Gerber, Nathan, The Second-Order 

ortional to the sum C,, On: ike ; Meee “es 

I meg ™q Lifting Pressure and Damping in Roll of Sweptback Rolling 

Airfoils at Supersonic Speeds, Journal of the Aeronautical 
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rical thickness distributions. But since the flow over matics and Physics, Vol. XX XIII, No. 1, pp. 46-56, April, 1954 


of a slowly oscillating airfoil is investigated in Appendix 
(B). It is shown that, to the second order in ampli- 
tude and thickness and to the first order in frequency, 
the damping of a slowly oscillating airfoil is directly pro- 
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APPENDIX (A) 


In reference 10 several expressions are given for the scalar potential function of the three-dimensional time 
dependent linearized partial differential equation of supersonic flow. In this appendix the results of reference 10 
are extended to include Eq. (11). The notation of reference 10 will be used in this appendix. 

Eq. (17) of reference 10 can be written as 


’ I R dd, R oer) | M [0¢: . Of 
V¢—(¢ — + “)¥ -) | +i, da = 
sf E , (+ B% dt, | Bc db) R ’ BXR & . | — 
y 2V O° 1 0° 2V Ode 1 O°d\ 1 
WJ (ore = = ua = . _ ; oe —- *) dv (Al) 
c* O£0f; c* Ot,” c? O£0ls c* Of.” R 


where the closed surface of the integral of the left side encloses the volume of the integral on the right side. 


The potential function, ¢, is required to satisfy Eq. (11). In this case Eq. (A1) reduces to 


oo lic ( EO | Boh yi + 
R ? ? B’c Ot B°c Obs R 
_ M dd, oer) 4 /- : ‘ 
- ft, da = (F, + Fs) dv (A2) 
’ BCR (* Le) Rehedee Y) ) oat e 


Eq. (A2) will be applied to a volume (denoted by w%) enclosed in the forward Mach cone from the point (x, ¥y, 3) 
(see Fig. 5). This volume is bounded by the forward Mach cone, an arbitrary surface S; enclosed in the forward 
Mach cone, and a surface given by = x — 6, where 6 is small. 

The finite part of the surface integral is zero, since on the Mach cone f; = fs. 
& = x — 6 reduces to a time-independent problem in the limit as 6 approaches zero. 


lim ff |; a —¢ x () | da = — 27r(x, y, 2, t) 
70 R On, On, \R 


The details of the above integration may be found in reference 11 or in reference 12. 
It follows from the preceding paragraph and Eq. (A2) that the potential function at the point (x, y, z) can be 


The surface integral over the area 
Thus this surface integral in 


the limit becomes 


expressed as 


l _— ; 
d(x, y, 2,¢) = — re Wy R (F, + F.) dv + 


| ff F ( R 09; R 2s ) (;) oM = ag! 
— = ‘fi, da (A3 
on IIs, LR’? \9~ Be on 7 atc on)” \R) Tar \on * on) 


For planar problems in which the point (x, y, 2) lies above the z = 0 plane and the disturbing surface lies in the 


z = O plane, Eq. (A3) can be reduced to 


| I | : 
G.%5,43) == F F,) dv — 
P(x, ¥ ) ‘a » R's it fe) ¢ 
| l R 0¢, R 0¢,' R 0O¢» R =a] 
es ae ae — 1é d (A4) 
Qn If | )RT (¢ o — 5% at, | Be dh | Boh Be w/ | 


where the surface S, is the z = 0 plane and the potential function, ¢’, is a potential below the z = 0 plane. 
The potential function ¢’, as yet undefined, will be defined so that 


o' (x, ¥, —do, t) = G(x, ¥, ao, t) 
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where do 1s an arbitrary constant. In this case Eq. (A4) reduces to 


l fff I - l if (* oe) l * 
(x, y, 2,0) = — (Fi + Fe)dv — lE d (A5) 
PRAY edad. 2 TO ~ Soda ee * BE 


APPENDIX (B)—-THE RELATION BETWEEN THE C,,. PLUS THE C,, AND THE DAMPING OF OSCILLATING AIRFOILS 
a q 


The damping of oscillating airfoils calculated on the basis of the linearized theory can easily be shown to be di 
rectly proportional to the sum Cn, + Cm, calculated by the linearized theory. The question arises: Does this 
hold true for second-order calculations? This question will be answered by investigating the partial differential 
equations and the boundary conditions for the three types of motions involved. 

Eq. (1) can be used in the analysis of oscillating airfoils; however, for comparison purposes it is more convenient 
to express the equation in terms of axes fixed to the airfoil. The relations between the two sets of axes are (see 


Fig. 6) 
/ wl 5 wil 
x = x’ cosa,e” — 2’ sin ane’ 
° wl i wl 
z = x’ sin a,e” + 2’ cos ane’ (Bl) 
t = ?’ \ 


where a, is the maximum angle of attack of the oscillating airfoil and w is the frequency of oscillation. 
For the axes attached to the airfoil, Eq. (1) becomes (to the second order) 
' 2M? l iat = j | 1W 1w La) 
“PO + 95 — On — 5 te = eel | 8 +o tn — 5 te + Ge tS Oy 
2 


2tw w” o l ; 
— (2@,, — x®.,) + — (26, — x@,) | + M2] (y — 1) (4+ &,) (,, + #,,) + 


Cc Cc 
») ») 
26,9,, + 26,8,, + — O,%,, + —®,9,,| (B2) 
. | | 
For slowly oscillating airfoils the time variation can be expanded in the form 


ee = 1 + swt — (w¥?/2) +... 


The result of substituting this relation into Eq. (B2) and retaining only the terms up to the first power of w yields 


| 2M? . l 
— B°®,, + ®,. “ys V ®P,, ee ®P,, = Am2M? ®,, + pr Pet = 
Cc 
; t Z < l 2 x 
a,tw2M* | = pr Pee radi V ®:, + 1’ P,, T lV ®, i y2 ®,, + V2 | is 


2 2 
M4 () _ (4. aa *,) (Bee + Pez) + 2b,Pre + 26, b,2 + 7, OO + <P (B3) 


The boundary condition on the airfoil surface can be expressed as v-Vs = 0. In this case 


xt 


v= 21(1 + ¢,) (—ef,) +7 (a. = + @o.+ ®.+ an — anit 


Therefore, on the airfoil’s upper surface, 


¢: | 2 0 = —an — [amiw(x + Vt)/V] + ef, (B4) 


y.|. =) = r|: ofr —_ efo-: z=0 (Bd) 


The discontinuity in y across the leading-edge Mach sheet need not be evaluated, since the differential equation 
and the boundary conditions on the airfoil surface are sufficient to determine its value. 

The first-order boundary condition on the airfoil’s upper surface for a constant vertical acceleration and a steady 
pitching is [from Eq. (8), and Eq. (8) of reference 5] 


$.|2<0 = —a—[dat + (qx/V)] + é&f; (B6) 
Comparing Eqs. (B4) and (B6) it can be seen that, if we let 
a& = 7 = anla, (B7) 


to the first order, the oscillating motion can be considered as the sum of a constant vertical acceleration plus a 
steady pitching motion. 
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The second-order boundary condition on the airfoil’s upper surface for a constant vertical acceleration and a 
steady pitching is [from Eq. (9) and Eq. (9) of reference 5] 
¥.|. = e(," z=0 ar or" z=0) — ef(o-24 2 
where ¢’ denotes the first-order solution associated with pitching and ¢* denotes the first order solution associated 
with a constant vertical acceleration. Comparing Eqs. (B5) and (BS) it can be seen that the second-order bound- 
ary condition on the airfoil surface for a slowly oscillating airfoil can be considered as the sum of a constant ver- 


o + d2z'|2 =0) BS 


tical acceleration plus a steady pitching motion. 
The second-order partial differential equation for the combined motion of accelerating and pitching can be 





written as 


2M? l 2 
— B°P,, + ®.,. —_ V ®P,, — 
Cc 


- Oo, = 7 MP anlo(x®,, + &. — 2b,, + Vtb,. + tb.,) + 


— | l l l 
2M? J 2 (*. + Vy »,) (?,, = ®..) +- ?,P,, = P.P,. - \" ?,P,, + V P.P 4 (B9) 


where g and @ have been replaced by a,iw. Comparing Eqs. (B3) and (B9) it can be seen that they are not the 
same. The two equations agree except for the terms a,,2.M/°®,,, an(2/V)M°®,,, aniw(2/V)M°2®,,, and ap,iw X 
(2/V*)M*x®,,. These terms do not contribute to the part of the second-order solution, which is of the form 
anlweh(x, y, 2,t). This can be established by analyzing the form of the second-order solution in the same manner 
that the form of the solution for the accelerating case was analyzed previously in this paper. An analysis of the 
form of the second-order solution for the slowly oscillating airfoil will also show that only that part of the second- 
order solution which is of the form a,iweh(x, y, 2, t) will contribute to the lifting pressure which is out of phase 


with the instantaneous angle of attack and, hence, to the aerodynamic damping of the motion. 


It follows that, 


to the second order, the damping of slowly oscillating airfoils is directly proportional to C,,. + Cy». 
” 7 a q 





Tests of a Fairing to Reduce the Drag of a Supersonic 


Swept-Wing Root 


(Continued from page 178) 


positive drag, but, when the fairing was added, the net 
pressure drag of this part of the wing was approximately 
zero. 

This result shows the use of the fairing to produce a 
subsonic type of pressure distribution, giving good pres- 
sure recovery on the back part of the wing. Without 
the fairing the pressure drag is obviously much larger, 
corresponding to a supersonic flow pattern. 

Similar pressure drag results for a = 4° are shown in 
Fig. 5. Again the reduction of pressure drag due to the 
addition of the fairing is obvious. 

It should be noted that although the fairing eliminates 
pressure drag from the wing, the wing produces an ex- 
tremely small pressure drag on the fairing. Probably 
the chief advantage of the fairing is to eliminate the 
bad flow over the wing root. 


CONCLUSION 


The fairing is considered to have fulfilled its purpose 
of maintaining a subsonic two-dimensional flow pattern 
up to the junction of wing and plate for the following 
reasons. 

(1) Separated flow and a trailing vortex were always 
visible at the junction of the wing with an ordinary flat 
plate, but both were removed when the fairing was 


added. 


(2) The chordwise pressure distribution near the wing 
root was of a supersonic (accelerating) nature without 
the fairing and was of a subsonic nature (decelerating 
over rear half) with the fairing in place. 

(3) The pressure drag of the inner part of the wing 
was large without the fairing and approximately zero 
with the fairing in place. 

(4) This fairing will presumably reduce the induced 
wing drag, because it prevents the formation of a trail- 
ing vortex from the wing root, which must decrease the 
apparent wing aspect ratio of unfaired swept wings. 

(5) The fairing worked well at off-design angles of in- 
cidence and yaw. 

(6) Finally, the experiments show how these effects 
approach two-dimensional flow as we proceed in a 
spanwise direction away from the wing root. 
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Theoretical and Experimental Pressure 
Recovery of Sweptback Normal Shock Inlets 


GEORGE D. DICKIE, JR.* 
United Aircraft Corporation 


ABSTRACT 


A theoretical analysis of the performance of normal shock in- 
lets having sweptback faces is presented. The analysis, which is 
based on a continuity and stream-thrust balance without regard 
to the internal shock patterns, permits the prediction of the 
critical total pressure recovery of a sweptback normal shock inlet 
as a function of face angle, angle of attack, and free-stream Mach 
Number. The theory shows that the performance of a normal 
shock inlet at positive angles of attack can be significantly im- 
proved without penalizing the performance at zero angle of attack 
by sweeping back the inlet face. The theory also shows that 
the supercritical mass-flow ratio for the normal shock inlets is 
independent of Mach Number and depends only upon the geom- 
etry of the inlet and the angle of attack 

A series of normal shock inlets having face angles of 0°, 10°, 
20°, 30°, and 40° was tested at Mach Numbers of 1.62 and 1.99 
The experimental results are shown to be in good agreement with 


the theory 


NOMENCLATURE 


A = area 

D = diameter 

f = stream thrust per unit area 

F = stream thrust 

g = acceleration due to gravity 

m = mass flow 

m = Mach Number parameter defined in Eq. (3) 
\J = Mach Number 

P, = static pressure 

P, = total pressure 

R = gas constant 

7, = total temperature 

V = velocity 

W = airflow 

a = angle of attack 

8 = inlet face angle 

y = ratio of specific heats 

p = density 

@ = parameter defined in Eq. (14) and tabulated in 


reference | 


Subscripts 
() = free stream 
1 = inlet station 
2 = exit station at end of constant-area passage 
3 = end of subsonic diffuser 


INTRODUCTION 


6 bes VARIATION OF TOTAL PRESSURE RECOVERY with 
mass flow of normal shock inlets for jet engines 
operating at zero angle of attack at supersonic speeds 
is generally well known. In particular, the critical 
pressure recovery, defined as the maximum pressure re- 
covery that can be obtained at maximum airflow, can 
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be predicted from normal shock relationships by apply- 
ing a correction for the subsonic diffusion losses that 
occur downstream of the shock. Tests of normal 
shock inlets have shown that significant reductions in 
the critical pressure recovery are obtained at angles of 
attack other than zero. This angle-of-attack effect, 
which is generally obtained with all types of nose inlets, 
is particularly undesirable, since high engine thrust is 
required to overcome the induced drag at angle of 
attack. Since high pressure recovery is required for 
high engine thrust, some method of improving the 
angle-of-attack performance of normal shock inlets is 
desired. 

Visual observations made during previous investi- 
gations at the Research Department of United Aircraft 
Corporation have shown that the detached shock 
wave at the entrance of normal shock inlets during 
subcritical (reduced mass-flow) operation at angles of 
attack other than zero was neither parallel to the face 
of the inlet nor perpendicular to the free-stream flow 
direction. Measurements of the shock-wave angles 
during subcritical angle-of-attack operation showed 
that these waves were inclined at angles that approxi- 
mate the angles of the strong, or second, family of two- 
dimensional oblique-shock waves. These observations 
indicate that partial turning of the flow was accom- 
plished externally by an oblique shock. Consequently, 
it was concluded that significant improvements in the 
angle-of-attack performance of normal shock inlets 
might be obtained by inclining the face of the inlet 
such that the expelled shock wave would be parallel to 
the face of the inlet and all of the turning of the flow 
would be accomplished externally. 

The investigation reported herein was conducted 
at the Research Department of the United Aircraft 
Corporation to provide a theoretical basis for evalu- 
ating the performance of normal shock inlets having 


swept faces. 
THEORETICAL ANALYSIS 


Determination of Critical Pressure Recovery 


The theoretical analysis that was used to determine 
the critical total pressure recovery of a normal shock 
inlet having a sweptback face is based on a continuity 
and stream-thrust balance between two stations of the 
inlet. No assumptions were made regarding the type 
of shock pattern between these two stations. A sketch 
of the inlet is presented in Fig. 1 to show the various 
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flow parameters used in the analysis. Station 2, which 


is designated the exit station, is located at an unspeci- 
fied distance from the inlet, Station 1, at the end of 
the constant-area passage. Wall friction forces in this 
passage are assumed to be zero. The equations that 
satisfy the conditions of a continuity and stream- 
thrust balance are derived in the following manner. 
The entrance stream tube area of the inlet, Ay, changes 
with face angle and angle of attack whereas the exit 
area, A», is constant 
Ay = (rD,.?/4) [cos (8 — a)/cos B] 
A» = rD»? } 
Therefore, the inlet-to-exit area ratio is 
Ay/A» = cos (8 — a)/cos B (1) 
From continuity considerations, it is known that 
Wo _ W, = W. 
It can be shown that 
, = ; Fa . 
W=PAm/vVT, (2) 
where 
» Ps a ; 9) ‘ 
m = (P,/P,)gMV (y/R) 31 + [(y — 1)/2]M?} (3) 
Eq. (3) is a function of Mach Number only and is tabu- 


. ° Th . . . bad 
lated in reference 1. The continuity equation then 


becomes 
P,Aotin/WVT, = P,,Aste/+/T, 
but, since 
T, = 71, 
the total pressure recovery of the inlet is then 
P,./P, = (Aomto)/(Aome) 
or, by substituting Eq. (1), 


P,,/P, = |cos (8 — a)/cos B| (i10/ me) (4) 


The stream thrust in a given direction at any station 


is defined as 
F = PA + mV 
or 
F = P,A(1 + yM?) (5) 
The stream thrust at Station 2 is, therefore, 
F, = P,,AxX(1 + yM,") (6) 


Since the flow is not parallel to the inlet axis at Station 
1, the following substitutions are made. 
V1 = Vo COS @ 


a 
Pr, = fe 
A, = A;/cos B = Ao/cos (8 — a) 


The stream thrust parallel to the inlet axis at Station | 
is then 


«J 


F, = P,,A; cos B + m,Vo cos a ( 
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In order to satisfy the continuity equation, 
Mm, = M = poAoVo 
Therefore, Eq. (7) becomes 
F, = P,,A2 + (piAoVo) Vo cos @ 

This reduces to 

F, = P,,Ao[1 + y(Ao/A2)Mo? cos a] 8) 
Substituting Eq. (1) in Eq. (8), 
F, = P,,A2}1 + yM,2[cos (8 — a)/cos B] cosa! (9 


Since the stream thrust parallel to the inlet axis at 
Station 1 must equal the stream thrust at Station 2, 
Eqs. (6) and (9) are equal. 


P,Ax(1 + yMs*) = P,As 
\1 + yM,?[cos (8 — a)/cos 8] cos a} 


By rearranging and introducing total pressure, the 
equality reduces to 


Pe _ yee 4 op ee ~ * | 
P,, IP, y-Ho a COS @ { 


ae 
a oe Me | (10) 


To facilitate computations, two flow parameters that are 
tabulated in the tables of reference 1 were used in this 
One of these is the ratio of stream thrust 
The stream thrust 


analysis. 
per unit area to the total pressure. 
per unit area is by definition 

f = F/A = Pl + yM*) 
then 

S/P, = (P;/P,) 1 + yM?) 11) 
Substituting in Eq. (10), the total pressure recovery 


then becomes 


F.. _ jP., ' v2 £2 (B — a) a \ / f 
P,, ie VP + YMo aan COS @ (/ P,); 


lo 
(12) 


Since Eqs. (4) and (12) are both expressions for the 
total pressure recovery, equating them and rearranging 
provides the following expression. 


k P,) ro | cos 8 4 My | 13 
= stig~ COS (15) 
m 9 mP, Leos (8 — a) ’ wee 


Since both f/P, and m are functions dependent only 
upon Mach Number, their ratio is a function of Mach 
Number and is proportional to a parameter, ¢, listed 
in reference 1, where 


@ = 0.41942 [(f/P,)/m] (14) 


Substituting Eq. (13) in Eq. (14), 


Fr. ‘os B 
gd. = 0.41942 ( ° YI = + y\," cos «| 
m)P7,/ Leos (8 — a) 


(15) 
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Fic. 6. Comparison of theoretical and experimental values of 
critical pressure recovery versus angle of attack for sweptback 
normal shock inlets at Mo = 1.62. 


The calculation procedure consists of solving Eq. (15) 
for values of ¢» at different angles of attack for a given 
face angle and free-stream Mach Number. Then, from 
the tables of reference 1, values of m7» or (f/P,)2 can be 
obtained which permit solution of Eq. (4) or Eq. (12) 
for the total pressure recovery. 

The results of this analysis are presented in Fig. 2 
to show the critical pressure recovery as a function of 
angle of attack for various face angles at free-stream 
Mach Numbers of 1.40, 1.62, and 1.99. Since the an- 
alysis is based on the assumption of continuity, the re- 
sulting curves apply when the oblique shock either 
lies along the face of the inlet or is inside the inlet at 
the critical shock condition. The maximum value of 
inlet face angle shown in Fig. 2 for each of the three 
Mach Numbers represents the limiting condition where 
the shock wave lies along the face for all positive angles 
of attack. For cases where the oblique shock lies 
along the face of the inlet, calculations based on con- 
ventional two-dimensional oblique-shock theory will 
yield the same results as those calculated by the above 
analysis. For face angles larger than the maximum 
values shown in Fig. 2, the oblique shock is always ex- 
ternal to the inlet, and the total pressure recovery is 
independent of face angle and is dependent only upon 
the angle of attack and free-stream Mach Number. 
The pressure recoveries calculated by two-dimensional 
oblique-shock theory for both the first and second 
families of shock waves are shown by the dashed lines 
in Fig. 2. Since the two-dimensional oblique-shock 
theory can be used to determine the critical pressure 
recoveries only when the inlet shock lies parallel to the 


by sweeping back the inlet face. 


Determination of Supercritical Mass-Flow Ratio 


The supercritical mass flow of a normal shock inlet is 
defined as the mass flow through a stream tube in the 
undisturbed flow whose area, Ao, is equal to the pro- 
jected area of the inlet in a plane normal to the free- 
stream flow direction. The supercritical mass-flow 
ratio is defined as the ratio of the supercritical mass 
flow at a given angle of attack to the supercritical 
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Fic. 7. Comparison of theoretical and experimental values of 
critical pressure recovery versus angle of attack for sweptback 
normal shock inlets a Mo = 1.99. 
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Fic. 8. Comparison of theoretical and experimental values of 
supercritical mass-flow ratio versus angle of attack for sweptback 
normal shock inlets. 


by definition, the supercritical mass-flow ratio is 
Mo/ Mog=0 = pvAo | 0 (pA J e= = Ay Ava “a. % 16) 


From Fig. 1, 


Ava “@ = A» 


Also, from Eq. (1) 


Ay/A» = cos (8B — a)/cos B 


Substituting in Eq. (16), the supercritical mass-flow 


ratio becomes 
Mo/Mog=) = cos (B — a)/cos B (17) 


Therefore, it can be seen that the supercritical mass- 
flow ratio of sweptback normal shock inlets is inde- 
pendent of free-stream Mach Number and depends 
only on the geometry and angle of attack of the inlet. 
Theoretical curves of supercritical mass-flow ratio as 
a function of angle of attack and face angle are shown 


in Fig. 3. 


EXPERIMENTAL RESULTS 


To check the validity of the theoretical analysis de- 
scribed above, a series of sweptback normal shock inlets 
was tested in a 4.5- by 4.5-in. supersonic wind tunnel. 
The tests were conducted with the models mounted on 
the special airflow measuring system shown schemati- 
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cally in Fig. 4. With this system, the airflow through 
the model is measured by means of standard bellmouths 
that are attached to the end of an insulated plenum 
chamber. The total pressure recovery of the model is 
computed from calibrated free-stream conditions and 
from the measured airflow, static pressure, and cross- 
sectional area at the diffuser exit. The models are 
supported by a strut from a circular plate on the side- 
wall which permits rapid changes in angle of attack. 

The normal shock inlets used during this investiga- 
tion had circular cross sections and sharp lips. A con- 
stant-area passage approximately one model diameter 
in length was provided downstream of the lip. The 
faces of the inlets were swept back at angles of 0°, 10 
20°, 30°, and 40°. Each model was tested at free- 
stream Mach Numbers of 1.62 and 1.99 over a range of 
angles of attack. 

Schlieren pictures of the 20° sweptback inlet with the 
shock fully detached at angles of attack of 0°, 6°, and 
12° at a Mach Number of 1.99 are presented in Fig. 5. 
At an angle of attack of 12°, the expelled shock wave is 
parallel to the face of the inlet as predicted by the theo- 
retical curves of Fig. 2. 

A comparison of the theoretical and experimental 
values of critical pressure recovery as a function of 
angle of attack for each of the models tested is shown 
in Figs.6 and 7. The theoretical curves were corrected 
for a subsonic diffuser efficiency of 83 per cent. This 
value was determined from calibration tests of the un- 
swept inlet at zero angle of attack. It can be seen 
that the agreement between experiment and theory is 
good, especially at angles of attack near those for maxt- 
mum critical pressure recovery. Although no precise 
explanation of the discrepancy between theory and ex- 
periment at angles of attack other than those near the 
peak of the curves can be given, it is believed that the 
discrepancies result from changes in the subsonic 
diffuser efficiency of the models over the range of angles 
of attack. Such changes in subsonic diffuser efficiency 
are to be expected, since the velocity profile of the flow 
entering the subsonic diffuser is probably changed with 
angle of attack as a result of shock-wave boundary- 
layer interactions in the inlet. 

Experimental values of the supercritical mass-flow 
ratios obtained during these tests are compared with 
theory in Fig. 8. It can be seen that the agreement is 
good at all angles of attack and that the supercritical 
mass-flow ratio is independent of Mach Number as 
predicted by the theory. 
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Buckling Efficiencies of Plate Materials at 
Elevated Temperatures 


‘ 


GEORGE GERARD* 
New York Unwersity 


SUMMARY 


The structural efficiencies of plates that fail as a result of elastic 
or plastic buckling are considered. Several materials are evalu- 
ated on the basis of short-time elevated temperature properties 
using structural efficiency criteria that depend upon whether 
failure is due to elastic or plastic buckling. Efficient temper- 
ature ranges of application for the materials considered are de- 
termined. Relative weight penalties resulting from decreased 
physical properties at elevated temperatures are indicated. 


NOMENCLATURE 


a length of plate, in. 

6 = loaded edge of plate, in. 

E = elastic modulus, lbs. per sq.in. 
E plastic modulus, Ibs. per sq.in. 
E, = secant modulus, Ibs. per sq.in. 
K = buckling coefficient 

N = applied loading, lbs. per in. 

t = plate thickness, in. 

W = weight, lbs. 

a = coefficient 

p = density, lbs. per cu.in 

» = structural efficiency parameter 
@ = failing stress, lbs. per sq.in. 

oy = compressive yield stress, lbs. per sq.in. 
ocr = Critical stress, lbs. per sq.in. 


INTRODUCTION 


hes STRUCTURAL EFFICIENCIES of various materials 
at elevated temperatures have become of increas- 
ing interest in anticipation of thermal flight of super- 
sonic aircraft. Several publications have appeared 
recently which have considered the temperature ranges 
for which various aluminum, titanium, and steel alloys 
would have efficient application. 

Hibbard,' among others, has compared various ma- 
terials based on short-time ultimate tensile strengths. 
Heimerl and Barrett? evaluated the efficiency of 20 per 
cent cold-worked titanium as compared to aluminum, 
magnesium, and steel alloys. The criterion used to 
evaluate the structural efficiencies of column and plate 
elements subject to buckling was the stress-density 
ratio as a function of a structural index appropriate to 
the type of buckling. Evaluations were made on the 
basis of short-time stress-strain characteristics of vari- 
ous materials for a series of selected temperatures. 
More recently, Heimerl and Hughes? have investigated 
the compressive, stiffness, and plate and column buck- 
ling efficiencies of the more promising aluminum, ti- 
tanium, and steel alloys. An analysis of buckling 
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strengths of plates and columns for large and small 
loadings similar to that of reference 2 was made over a 
temperature range up to 1,200°F. using the stress- 
density ratio as a criterion of efficiency. 

The buckling of plates under load for short-time 
elevated temperature conditions is similar to room 
temperature conditions. Under elevated temperatures, 
Heimerl* has shown that for plates that buckle elasti- 
cally, the failure stress is related to the buckling stress 
and compressive yield stress of the material by the same 
relationship that holds for room temperature. Simi- 
larly, failure of plates that buckle plastically is ade- 
quately predicted by substitution of the appropriate 
value of the secant modulus in the buckling stress 
equation. 

Since the short-time buckling behavior of plates is 
consistent over a wide range of temperature, it is the 
objective of this paper to consider the structural effi- 
ciencies of plates that buckle elastically and those which 
buckle plastically. It is found that different weight- 
strength criteria are obtained for the two cases. 


BUCKLING EFFICIENCY ANALYSIS 
The weight of a plate element is 
W = ptab | 
The buckling stress of the plate 
On = KE (t/b)? (2) 
For a structure of minimum weight it is assumed that 
the applied loading, N, and failure stress o are related 


by 


By combining Eqs. (2) and (3) 


t = (Nb?/KE) (4 


Elastic Buckling 


For plates that buckle elastically, the value of the 
modulus & = E appropriate to the temperature under 
consideration. Since the failing stress of plates that 
buckle elastically is related by the same equation at 
room and elevated temperatures,‘ the relation proposed 


by Schuette® can be used 


- V/, 3/5 - 
= 2 Ce’ Cee" (5) 
By substituting Eqs. (2) and (3) into Eq. (5) 
rrie/ e/a, 2/5 pp al/e 1/o >. 
b= CNR a) ee ae) (6) 
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FIGURE / 
WEIGHT- STRENGTH EFFICIENCIES OF PLATES WHICH 
FAIL BECAUSE OF ELASTIC BUCKLING 
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Introducing Eq. (6) into Eq. (1), the weight of the plate 
that will fail at the applied loading is 


W = (N/aK')*“ab‘"(p/E) “oe,” (7) 


In comparing the weights of two structures designed 
of different materials for the same loading and geomet- 
rical parameters, the values of NV, a, 6, and K are fixed. 
Test data for a wide variety of aluminum and mag- 
nesium alloys at room temperature’ and for 75S-T6 at 
elevated temperatures‘ indicate that the coefficient a is 
constant to a good approximation for comparable struc- 
tural elements of different materials. Therefore, the 
remaining variables that are a function of the physical 
properties of a material are p, EF, and o,,. Thus, the 
proper efficiency criterion for evaluating the failure of 
plates that buckle elastically is, from Eq. (7) 

te = (0/E) "ey" (8) 


Plastic Buckling 
For plates that buckle plastically, the failure stress at 
both room and elevated temperatures can be obtained‘ 


TABLE | 
Room Temperature Properties of Alloys 


Cc, 
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by substituting E, = £ in Eq. (4). Thus, the weight 
of a plate that fails by plastic buckling can be obtained 


by combining Eqs. (1) and (4). 
: W = (N/K)*ab’*(p/E,”") (9) 
For this case, the efficiency criterion is simply 


a» = /E,” (10) 


RELATIVE EFFICIENCIES 


The short-time elevated temperature compressive 
properties usually stated for various materials are the 
compressive yield strength, o,,, and elastic modulus, £. 
Based on such data, given in references 3 and 6, the 
structural efficiencies of the aluminum alloys 245-T4, 
75S-T6, and XA78S-T6, titanium alloy RC-130-A, the 
steel Stainless W, and Inconel X were computed ac- 
cording to the criteria given by Eqs. (8) and (10). 
The room temperature properties of these alloys are 
given in Table 1. 

Fig. 1 represents a structural efficiency comparison 
of plates that fail by elastic buckling. It can be ob- 
served that according to the criterion used, aluminum 
alloy XA78S-T6 is most efficient up to 400°F. Between 
400°F. and approximately 970°F., RC-130-A is superior, 
and beyond 970°F. Inconel X is the most efficient ma- 
terial of those considered. It is to be noted that Stain- 
less W is competitive with RC-130-A in the temperature 
range of 400°F. to S00°F. 

To use the efficiency criterion for plates that fail by 
plastic buckling, it is necessary to specify the values of 





FIGURE 2 


WEIGHT-STRENGTH EFFICIENCIES OF PLATES WHICH 
FAIL AT THE COMPRESSIVE YIELD STRESS 
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FIGURE 3 
RELATIVE WEIGHT OF EFFICIENT MATERIALS AT 
ELEVATED TEMPERATURES 
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the secant modulus in Eq. (10). In elevated tem- 
perature applications, it may be efficient to design the 
heavily loaded portions of the structure to the com- 
pressive yield. Therefore, the value of the secant 
modulus is taken as that corresponding to the com- 
pressive yield stress in this analysis. 

Since the yield is defined as the 0.002 in. per in. offset 
stress, the secant modulus at yield is 


Exey = Fey/ (0.002 + ocy/E) (11) 


Scy 


By use of Eq. (11), only the values of o,, and F are re- 
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quired to compute the secant modulus at compres. 
sive yield. It is not necessary to use the stress-strain 
curve in this case. 

In Fig. 2, a structural efficiency comparison for 
plates that fail at the compressive yield stress is pre- 
sented based on Eqs. (10) and (11). In contrast to 
the previous case, aluminum alloy XA7S8S-T6 is effi- 
cient to 250°F., after which 24S-T4 is superior up to 
700°F. Titanium alloy is efficient beyond 700°F 

Fig. 3 was prepared to indicate more clearly the dif- 
ferent ranges of structural efficiency for the two cases 
considered. The relative weights at elevated tem- 
peratures were obtained by comparing the values of the 
weight-strength criterion of the efficient materials at 
elevated temperatures to XA7S8S-T6, which is the most 
efficient material at room temperature. It can be ob- 
served that the temperature ranges of efficient appli- 
cation of a given material are appreciably different for 
the two cases considered. Fig. 3 also indicates some 
of the weight penalties associated with the design of 
aircraft for thermal flight. 
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Extension of the Statistical Approach 
to Buffeting and Gust Response 
of Wings of Finite Span 


H. W. LIEPMANN* 
Douglas Aircraft Company, Inc., and California Institute of Technology 


SUMMAR\ 


The response of a wing of finite span moving in a turbulent 
field is discussed. On the basis of thin-wing theory and using 
concepts from the theory of homogeneous turbulence, the mean 
square lift can be expressed in terms of a generalized aerodynamic 
admittance and one component of the power spectrum tensor of 
turbulence. 

Explicit results are given for several limiting cases. 
that the mean square lift coefficient of a wing with a span large 
compared to the scale of turbulence is simply proportional to the 
The analytical procedure 


It is found 


ratio of scale of turbulence to span. 
for the general case of a finite wing in a turbulent field is outlined 


(1) INTRODUCTION 


¥ A PREVIOUS PAPER’! statistical concepts were used 
to formulate and study the buffeting of an airfoil 
of infinite span due to turbulent fluctuations. The 
spatial structure of the turbulent field was neglected. 
In the present paper the analysis will be extended to 
wings of finite span, and the structure of the turbulence 
will be taken more fully into account. It is intended 
to formulate the problem in a way similar to that in 
reference 1 and to obtain explicit solutions for some 
simple idealized cases. 

The main purpose here is the study of the influence of 
finite aspect ratio and of the three-dimensional char- 
acter of the fluctuations on buffeting and gust response. 
As in reference 1, the emphasis will be mainly on the 
concepts and on the general features of the problem 
rather than on detailed computations of specific cases. 
Since the concepts used are closely those of reference 1, 
it is unnecessary to repeat much of the detailed dis- 
cussion. References to the general statistical concepts 
are also given in reference 1 and will not be repeated 
here. 


(Il) GENERAL APPROACH TO BUFFETING OF WINGS OF 
FINITE AND INFINITE SPAN 


To formulate the problem of the present study, a 
condensed version of the approach to the buffeting of 
wings of infinite span is given first and then generalized 
to wings of finite span. 

Buffeting is considered the response of a linear sys- 
tem to a random force field. The force field is due to 
the random fluctuations present in turbulent flow pro- 
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f(w), ete. 


duced by the wake of a body or by the gustiness of the 
air. It is then assumed that the relation between, say, 
the fluctuating angle of attack of the wing and the re- 
sulting lift is given by a linearized theory—e.g., thin- 
wing theory. 

The problem of the response of two-dimensional air- 
foil as formulated in reference | can then be expressed 
as follows: The airfoil is located at the origin of a co- 
ordinate system x, y, 2; the mean velocity of flow is L’ 
and is directed along the positive x axis. The airfoil 
is infinitely extended in the y direction. Within the 
frame of nonstationary thin-wing theory the airfoil is 
replaced by a /ifting point at the origin, sensitive to 
angle of attack fluctuations in the x, z plane. The 
effect of the finite chord, c, is absorbed in an admittance 
function that relates the frequency and magnitude of 
the lift to the frequency and magnitude of the angle of 
attack fluctuations at the origin. The problem is thus 
formulated in terms of a one-dimensional, stationary, 
stochastic process. That is to say, the lift, LZ, of the 
airfoil is related to the fluctuating angle of attack, 
a(t), by a linear operator. a(t) is assumed to be a sta- 
tionary, stochastic function defined by certain mean 
values like the mean square a’, the power spectrum 
The response of the airfoil is represented by 
an “admittance” or “impedance”’ function, ['(w), giv- 
ing the lift of the airfoil corresponding to a sinusoidal 
fluctuation of the angle of attack. Often it is more 
convenient to use the indicial admittance, h(t)—that 
is, the lift response to a delta function variation of the 
angle of attack. The latter, /(¢), is the Fourier trans- 
form of I'(w). The lift L(¢) at a time ¢ is then related 
to the angle of attack by a Duhamel integral: 


Lit) = { drh(r)a(t — r) (1) 


from which the mean square lift is obtained as 


L? = ff dr, dro h(7)h(72) a(t — ry)a(t — 72) (2) 


a(t — m)a(t — Te) = y( T1 — Te!) 


is the auto-correlation function of a and thus related to 
the power spectrum of a by 


¥(r — o) -{ dw f(wee — ? (3) 
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Introducing Eq. (3) into Eq. (2) and inverting the order 
of integration yields then 


“= ae = >| I (w)| 2 f(w) ) 
Z(w)|? | os tw) | Se 


where the “impedance” Z or the “‘admittance’’ I’ is 


given by 


l + « 
Tiw) = Z(w) = i- dr e'*" h(r) (5) 


Thus to compute L?, the spectrum of the fluctuation has 
to be obtained from turbulence measurements and the 
‘“admittance’”’ from nonstationary wing theory. 

To generalize these results to wings of finite span 
means to take into account both the fact that the fluc- 
tuations at any instant will vary along the span and 
that the admittance function / of the wing depends on 
the spanwise coordinate y. 

The problem can thus be formulated as follows: The 
wing is replaced by a /ine extending from —db to +0) 
along the y axis. The response of this line is described 
by an indicial admittance h(t, y). The fluctuating 
angle of attack a will now be a function of both / and y. 
Thus generalizing Eq. (1), we write for the contribution 
to the lift of an element dy of the span 


+ = 
dL = ay f dr h(r, y)a(t — 7, y) (6) 


L and L? thus become 


+b + « 
L= f f dy dr h(r, y)a(t — 7, y) (7) 
_" _— 
+6 7+ « 
L? = ff ff dy, dy2 dr, drz h(11, vidh (72, V2) X 
= =f 


a(t — 7, yi)a(t — 7, Vo) (3) 


To proceed from Eq. (8) we assume that in addition to 
being stationary—that is, homogeneous in time—the 
fluctuations are homogeneous in space. Then 
a(t — nn, W)alt — T2, ve) = W(r1 — T2 V1 — Ve) (9) 
that is, the correlation function Y depends upon the 
difference of the variables only. It is now a simple 
matter to generalize Eq. (3). In order to keep within 
the standard nomenclature of turbulence (e.g., references 
2 and 3) one introduces a correlation tensor R,; defined 
by 
Ri = uj(r)uj(r + s) (10) 
= ai 
and the corresponding spectrum tensor ¢;; (k), where k 
denotes the wave number vector 


Ri;(s) - fff (Re am dkdkodks => 
[ | dV(k)bi(k)e** (11) 


Furthermore, in relating y to R;;, use is made of the fact 
that within certain limitations (discussed in references 
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1 and 2, for example) the difference of time 7; — 7» in 
Eq. (9) can be replaced by a difference in the x-coordi- 
nate (7; — 72)U = x, — x2. Now the angle of attack 
a is related to the z component of the fluctuating ve- 
locity by 


a= U3 U 


Hence y corresponds to the correlation R3;—1.e., the 
correlation between the two z-components of the fluc- 


tuating velocity taken at two points a distance s apart, 
ee 


where s has the components 


$= S(% — Xe, Vi — Ve, V) 


The turbulent velocity fluctuations 1, u#2, “3 are func- 
tions of x, y, z,¢. But because in thin airfoil theory the 
boundary conditions are applied in the plane of the 
wing, z = 0, a = u;/U enters only as a function of x, 
y, t. The previous transformation has related ¢ to x, 
and, hence, the correlation function for a@ depends 
finally on two variables only. 

Using Eq. (11) ¥ can thus be expressed by 


V(71 — T2 Vi — Vo) = R(x — Xe, V1 — Ye, O) = 


{ dV(k) a3(k)e'! , got y2) ho) (12) 


Inserting Eq. (12) into Eq. (8) yields the proper general- 


ization of Eq. (4) 
L? = | dV (Rk) b33(R)| T (Ri, Re) |? (13) 


with 


l e+h x , on 
rie) = = | dx dy h we" | 6h 
l = ‘il 


(III) Lrmitinc Cases; IDENTIFICATION WITH PREVIOUS 
RESULTS 


(A) Consider first the case where the airfoil admit- 
tance / is independent of y—that is, two-dimensional 
airfoil theory is used—but the airfoil itself moves in a 
homogeneous field of three-dimensional fluctuations. 
The results that follow are more general than those 
given previously! but still do not take into account the 
effect of the finite span in the response to the fluctu- 
In this case Eq. (14) becomes 


1 te ev. 
U dx h U eo ae 


+b 1 a ‘ 
f dy eth?) on 1 (R1) U7 | dy elke (14a) 
—b - 


I'(ky, Ro) = (2/U)T,(Rki) (sin kob/ke) 


ations. 


T'(ky, Re) 


and thus Eq. (13) becomes 


4 — , + — sin? kob 
i? = = | dk; |T,|?2 {f dk dks 33(k) Bo? 


(13a) 
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From Eq. (13a) one obtains the two interesting limit- 
ing cases for b large and 6 small compared to a scale of 
turbulence. The former case corresponds to a wing of 
extremely large span, so large that the main contri- 
bution to the lift comes from turbulence components 
with wave lengths much smaller than 6. The latter 
case corresponds to a wing with span small compared 
to a scale of turbulence but for which two-dimensional 
theory is used to obtain the response of the wing. 
Hence in this case the fluctuating angle of attack at any 
instant is constant across the span. This is the case 
previously treated. ! 

For the first case introduce the variable 86 = Dk. 


Then Eq. (13a) becomes 


4 7. + o 
L? = l? If dk; dk3\ T\|\*b X 


{ ws (; Bb, sin? B 
a ap P33 Iv py Ka 32 


which, for large values of 6, tends toward 


arb Pr 
lL? = U2 ff dk, dks Tr; . 33(R1, 0, ks) (l3aa) 


In the second case, where bk. < 1, the sine in Eq. 
(13a) may be expanded in series and one obtains 
1h? 


[? = [72 } 7 dk; Vy : i] | dko dks $33(ki, ko, ks) 


Now from Eq. (11) it follows that 


U3? = R33(0) = iff dk, dk» dk; 33(R1, ko, kz) 


while on the other hand y%;? is expressed in reference 1 


by 


us = U? | f(o) da = F(k,) dk, 
/7 0 /7 0 


Comparing these two equations 


Fu) = 2 {f dk» dks $33(R1, ko, k3) 


so that for small 6 the mean square lift can be written 


Dp2 4 

Ll? = = f dk; 1(k1)| ?F(Ri) (13ab) 

This agrees exactly with the result given in reference | 
except for a slight change in notation. 

It should be noted that the root mean square lift is 
proportional to +/ for large 6 and proportional to 
forsmall 6. This is a general feature of the behavior of 
linear systems under random loading. 

(B) Retaining the assumption that h/ is independent 
of y, one may simplify the problem further by assuming 


~~» 
the fluctuations to be isotropic. In this case $33(R) can 
be obtained from the general form of the spectral tensor 


_— 


¢:;(k) for isotropic, homogeneous turbulence :* 


— E(k) .. j=l 4=7 

b. (Rk = [(R76;;) a (Rik ’ bay" ; 

ee yy N20 ij 
(15) 


where E(k) is the so-called ‘“‘energy’’ or ‘‘Heisenberg”’ 
spectrum and is related to the total kinetic energy per 
unit mass by 


») 
0” 
9 


u2 = f E(k) dk, k? = ky? + ko? + ks? 
0 


— 


From Eq. (15) it follows that for the component ¢33(k) 


33(k) = (E(k) {rk (k;" . 2 ko") (16) 
With ¢33; as given by Eq. (16) the inner part of the in- 
tegral in Eq. (13a) becomes 


+a — sin? kob 
if dks dks 33(k) ° = 


ky 


| If thy dy (by? + by?) C7 
dR» ak: (rR 2" \ 


By introducing polar coordinates in the kz, ks plane the 
double integral in Eq. (17) can be written as the product 


l f : E(x? + ky”) 
K dk 
Ar 0 (x? + k,?)? 


- 9 a — sin? (bx cos @) 
dd(k,? + x® cos? #) = 
0 (17a) 


Kk? cos? 


ko = xcos vd 


kg = xsind 


which can be transformed easily into another form often 


more convenient 


l f dx E(x? + k,?) - 1 sin? bkt ky? + xt? 
re 4 7 
2r Jo a (x* + &,?)? J; t° V1-£ 


(17aa) 


Both the forms (17a) and (l17aa) are useful because the 
second integral in these is independent of the spectrum 
of turbulence /(k) and is thus a function that can be 
tabulated once and for all.* Thus the mean square 
lift for a wing in isotropic, homogeneous turbulence, 
using the impedance corresponding to an infinite aspect 


ratio, becomes 


L? = a [ dk T,(ky) 2 { kK dk x 
rl Y a anal JV 


E(x? + ki?) 2" sin? (bk cos #) 

- _ = I dd(k\? + x? cos? #) = 

(x2 + ky)? Jo x” cos” # 
(18) 


(C) Finally, one may choose a simple spectrum of 
turbulence to obtain explicit results. A reasonable 


* The measurements of the spectrum of turbulence obtained 
with a hot wire of finite length 2) can be analyzed by an identical 
approach leading to formulas similar to (17a) or (17aa). Using 
gu instead of ¢33 in Eq. (13a), one obtains the equation found by 
Uberoi and Kovasznay* who employed a somewhat different ap- 


proach. 
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simple spectrum which applies approximately to iso- 
tropic turbulence is that used in reference 1, which, 
in the notation of this paper is 


$33(R) = (213" m”) A®[ (Ry? + Rk?) (1 + A*k?)3] 
(19) 


where A denotes the scale of turbulence. Using Eq. 
(19), the previous limiting cases for 6 large and 6 small 
can be given explicitly. Eq. (13aa) becomes 


U3? Bic Il. k,? 
(2 = 3 — DA4 1k,| 1 \(R,)| ? ; 
aiid J eV 1 At) 


and using for | l’,| ? the approximation 
(Ti(ki)|? = (9/4) ((pU%e)?/(1 + meki)}] (20) 


one obtains 


or ., U3" s? ds 
L? = — put pact 
7" wT’ J. i+ @a/sian+ 


(21) 


C 6 A f s’ds aes 
2 = Gr%q? — (22) 
7 eb So [1 + (re/A)s] (1 + 52)” 


Eq. (22) is valid for \ <b. If in addition A < cc, Eq. 
(22) can be simplified further and yields the simple re- 


- 70 


sult 

Cry = 2a?r(A2/bc) (23) 
If on the other hand the aspect ratio is extremely large 
and the scale of turbulence such that c << A <b), Eq. 


(22) can be evaluated also and yields 


Ci? = 29a2(A/b) [1 + 2e(c/A) +...] (24) 
These results are essentially independent of the form 
of the spectrum chosen provided the scale A is defined 
in a suitable way. 
If b is small, one uses Eq. (13ab) and inserts for 
U3°A ] a 3A°R)? 
wr (1 + A?k,*)? 


ll 


pi ff dks dks 33 ( k ) (25) 


F(ki) = 


which is the one-dimensional spectrum corresponding 
to Eq. (19). Thus Eq. (13ab) becomes 


I, (1) > dk, YH 


2. 2? 1+ 3A%,? 
L? = 45? sa \ | us a 
, ( 


Uz (1 + A%%,?)? 


and hence 


or with n = me A 


, us" a3 - f 
e . a Ag? i 
7 U? Fee + 1) - 


n+ 3 a i 
= > (n log n° + 7) 2S) 
2a" + 1)? 


which is identical with Eq. (16) of reference 1. 


(IV) GENERAL DISCUSSION 


In general it is necessary to obtain the two functions 


I (ky, kz) and ¢33(k) in order to compute the mean square 
lift from Eq. (13). The admittance I has to be ob- 
tained from thin-wing theory and represents essen- 
tially a generalization of Sear’'s results to wings of finite 


ab 
span. 3;(k) has to be obtained from a knowledge of 
the turbulent field. 

In limiting cases the problem can be much simpli- 


fied as illustrated by Eq. (22). 
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Prediction of Shock-Induced Turbulent 
Boundary-Layer Separation 


Artur Mager 

Engineering Specialist, Marquardt Aircraft Company, Van Nuys, 
Calif. 

October 1, 1954 


SYMBOLS 


H form factor (6*/6) 

VM Mach Number outside boundary layer 

p pressure 

Ap pressure difference across shock (ps pi 


dynamic pressure 

velocity outside incompressible boundary layer 
r velocity inside incompressible boundary layer 

distance along the surface 

distance normal to the surface 

oblique shock angle 

ratio of specific heats 

displacement thickness 
‘ V;2 sin? 8 -— 1 

Gruschwitz form factor 

momentum thickness 


Subscripts 


initial 


eparation 


INTRODUCTION 


pine RECENT EXPERIMENTAL information (reference 1) indi 
cates that the pressure coefficient required to separate a 
wholly turbulent boundary layer through shock is essentially 
independent of the Reynolds Number, but varies inversely with 
the Mach Number. 


this finding is in close agreement with the semiempirical rules 


The purpose of this note is to show that 


used to describe the incompressible turbulent boundary-layer 
behavior when these rules are modified by the use of the well- 
transformation. 


known Dorodnietzyn-Stewartson 


DISCUSSION 


In practical computations, the behavior of turbulent incom- 
pressible boundary layer is usually described in terms of momen- 
tum-integral equation, a skin friction law, and an additional semi- 
empirical relation describing the variation of the form factor. 
Although there are a number of these semiempirical relations for 
the form factor, almost all of them can be rearranged, so that 
the change in form factor with velocity outside the boundary layer 
is separated from the change with distance along the surface. 
For example, the well-known Gruschwitz relation can be put in 
the form: 

dx : 

dn = (0.005 — 0.0097) = * oe d(U?) (1) 

Now, in cases of very rapid diffusion, one may neglect the changes 

in the form factor due to the distance, so that the form factor re- 

lation may usually be integrated in a closed form. For instance, 
Eq. (1) becomes: 


(U/Us)? = n/n (2) 
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Since the initial and separation values of the form factor are 
known, it is thus possible to evaluate the velocity ratio l’.// 

necessary for separation of the incompressible, turbulent bound- 
ary laver, when the diffusion occurs with a very steep velocity 


To proceed with the above chosen example, since 


; (‘“.) a H-1 |"! 
aaa U > H(H +1 7 


one may obtain the appropriate values of n from the knowledge 


gradient 


of the more common //. As is well known, the initial value of 


H depends very slightly on the Reynolds Number (it varies in- 
versely with the Reynolds Number), so that for application in 


the range of Reynolds Numbers which is of interest, one chooses 


H; =~ 1.26 which gives 7 0.465. The separation is generally 
agreed to occur at H 2.6 or 7 0.941 Using these values, 
the Gruschwitz relation gives (U,/U;)? 0.494. This value 
is in agreement with the value given in reference 2 (L’,?/[ 


0.5) which was also derived from the Gruschwitz relation, but on 
the basis of somewhat different reasoning. It can be shown that 
almost all of the semiempirical relations give, in this fashion, a 
value of l’, which falls into this general order of magnitude 


It has been shown in reference 3 that use of the modified 


Dorodnietzyn-Stewartson transformation converts the com- 
pressible turbulent boundary-layer equations into the usual in- 
compressible form. This means that all the rules for incom- 
pressible boundary-layer flow will be applicable to the compres- 
sible flow, providing the significant incompressible quantities are 
interpreted in terms of their transformed compressible equiva- 
lents. Among others, this transformation requires that the in- 
compressible free-stream velocity L’ be thought of as the product 
of the stagnation speed of sound and the local free-stream Mach 
Number. Since the stagnation speed of sound through the shock 
is a constant, therefore the new compressible flow rule for this 


example is simply: 
V2 = 0.494./,? (4) 


From the knowledge of Mach Numbers in front and behind the 
shock it is possible to determine the pressure rise coefficient re- 
quired for separation. This is particularly easy if it is realized 
that since the boundary layer is incapable of sustaining very 


large pressure gradients, the usual linear approximation: 
M,? sin? 8B = 1 + ¢; e< 1.0; M; > V2 5) 


should give accurate results. It is fairly elementary to show by 
use of this approximation that: 

Ap Appi . 24 2 

qi Pi i ri var 
The small quantity « may be evaluated in terms of 17; and .\/, 
by use of the relation given in reference 4: 


1 + [(y — 1)/2] M.? 
1+ [(y — 1)/2] M2 
(y + 1)? 


(27M? sin? B — y+ 1){y-1+4 [2A Mesintay} ‘? 


so that finally, after some manipulation, using Eq. (4): 
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FIGURE |. SEPARATION PRESSURE COEFFICIENTS 


Ap 0.506 
Pa (8) 


gag 1+ ((y¥-—1)/2|M2 

This last relation is plotted in Fig. 1 together with the data 
reported in reference 1. It may be seen in this figure that not 
only relation (8) is in fair agreement with the experimental data, 
which after all may be due to the particular semiempirical rela- 
tion chosen here, but that also the trend of the data is properly 
represented. Note that due to the restrictions on Eq. (5), one 
may expect Eq. (8) to give poor results for MW; < +/2, but this 
should not be interpreted as reflecting on the validity of Eq. (4). 

Finally, it is important to comment that Eq. (4) indicates 
that for normal shock and Mach Numbers smaller than approxi- 
mately 1.2, no shock induced separation will occur. 
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Comment on the Note ‘Electrical Resistance 
Networks for Beam and Column Problems’’ 


R. H. Scanlan 

Centre National de la Recherche Scientifque, Ch€atillon-sous- 
Bagneux, Seine, France 

October 12, 1954 


ie A RECENTLY PUBLISHED NOTE,! we neglected to acknowledge 
the existence of a somewhat similar work by Swenson.’ It 
is of interest to note the following points of difference: 

(1) Swenson obtained his beam and/or column effects through 
the direct application of moments (i.e., their electrical equiva- 
lents) while we work simultaneously with both loads and mo- 
ments, since our second, superposed network, absent in Swenson’s 
system, converts loads into moments for-us automatically. 

(2) The existence of the second line in our network permits us 
to face the question (not brought out in Swenson’s formulation ) 
of the interaction between support moments engendered by 
buckling and those moments due to lateral load (i.e., the ques- 
tion of distinguishing between reactive moments and active mo- 


ments). This is of particular importance in the study of the 
buckling of columns having other than ‘‘classical”’ end fixities, 

(3) Swenson employs negative resistances in solving column 
problems by successive potentiometer adjustment. We employ 
positive resistances at all times while using iteration procedures 
that resemble closely the standard matrix iteration technique, 
the critical load dropping out as the characteristic value. 

Apart from the above, we have found our iterative procedures 
suitable for the solution of thin beam lateral buckling problems of 
the type described by Timoshenko’ and more recently by Martin‘ 
and Salvadori. In the last three references beam rigidity is 
limited to uniform or to linear variation. Our system, like 
Swenson’s, is completely general in this respect, being composed 
of finite differences. 
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On the Intermittency of a Two-Dimensional 
Jet* 


William H. Hannum and Wayland Griffith 
Palmer Physical Laboratory, Princeton University, Princeton, N.J 
October 13, 1954 


TT INTERMITTENCY OF A two-dimensional jet has been in- 

vestigated with a hot-wire anemometer and the statistical 
probability for the location of the edge of the jet compared with 
measurements using a pitot tube. Intermittency refers to the 
variation in time and space of the boundary of a turbulent flow 
The location of the frent has been found to have a nearly Gaussian 
distribution for turbulent boundary layers,! wakes,? slipstreams,’ 
and axially symmetric jets.‘ This investigation shows that the 
location of the boundary of a two-dimensional jet is also very 
nearly Gaussian. 

The hot-wire signal obtained in a turbulent region is char- 
acterized by high-frequency oscillations that are easily dis- 
tinguishable from the relatively smooth signals arising in regions 
outside the jct. Intermittency may be measure in terms of the 
intermittency. factor y, introduced by Townsend? and defined as 
the fractional time spent by a probe in turbulent fluid. A direct 
determination of the intermittency can therefore be made from 
the percentage of time a turbulent signal is recorded, measured 
over a reasonable length of time, for any given position. When 
the velocity of the turbulent fluctuations in the stream is small 
compared with the mean motion, these irregularities will be 
swept over a stationary probe without deformation of the form 
of the front so that time and space fluctuations are interchange- 
able. 

A hot-wire unit was constructed along the lines of Kovasznay’s 
transformer compensated “poor man’s’ amplifier.5 We are 
grateful to Dr. Kovasznay for the loan of a suitable transformer 
The hot-wire used was Wollaston platinum wire 0.25 mil in diam- 

* This work was submitted by W. H. H. as a Senior Thesis to the Depart- 
ment of Physics, Princeton University. Part of the equipment used was 
made available through ONR Contract N6oril05, Task II We are grateful 
to Dr. G. B. Schubauer for the suggestion of this problem and for advice on 


procedure during the course of the investigation 
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Typical intermittency films ranging from fully turbu 
lent near the center of the jet (top) to fully nonturbulent outside 


Fic. 1 


the jet (bottom The markers below the trace are 60 cycle per 
sec. flashes. 

eter and (0.2 in. long. Since all of the experiments were carried 
out at low velocities and the hot-wire was used primarily to dis- 
tinguish between turbulent and nonturbulent flow, many of the 
refinements necessary under more severe conditions could be dis- 
pensed with. Before each run a sine (or square) wave signal was 
superimposed on the heating voltage and the circuit adjusted 
to give the best response as seen with an oscilloscope. 

Permanent records of the signal were made with a drum camera 
using 35-mm. Linagraph Ortho film. The film speed was from 
5to7 ft. per sec., the record at each position being about 2 sec. in 
duration. 

The nozzle used is 1 in. wide and 24 in. high 
trifugal blower passed into a 24- by 24-in. settling chamber and 
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Fic. 2. Intermittency factor y measured with the hot-wire 
and relative velocity distribution obtained with the pitot tube 
plotted as a function of the angle a defined in the inset. Data 
Were taken at three axial distances from the exit plane of the 
nozzle. The solid line is a Gaussian distribution. 
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thence through a smooth entrance duct to a 12-in. straight sec- 
tion before emerging as a horizontal jet. Traverses with a pitot 
tube showed that the flow is symmetrical and two-dimensional 
over the central 14 in. of the jet out to a distance of 18 in. from the 
exit of the nozzle. All of the data reported here were taken in 
the plane midway between the top and bottom of the nozzl 
Preliminary surveys showed that the jet was fully turbulent 
at the exit plane. 

Fig. 1 shows four typical intermittency films ranging from 
fully turbulent (top) to fully nonturbulent (bottom In the 
intermediate cases the distinction between the two types of signal 
is clearly discernible. Such records were made with a jet ve 
locity of 68 ft. per sec. at three distances from the exit plane 
(x = 6, 12, and 18 in 


angle a between a line joining the point in question to the edge 


for several lateral positions defined by the 
of the nozzle exit and the plane of symmetry. For each point 
the intermittency factor y7 was determined by measuring the 
fraction of the trace which was turbulent Fig. 2 shows the 
variation in y with a at the three axial stations rhe solid curve 
is a Gaussian distribution that fits the data rather well except 
near the inner edge. No attempt was made in this investigation 
to explain the deviation from a Gaussian curve 

The velocity distribution obtained from a pitot traverse is 
also included in Fig. 2. Here Vo is the velocity at the center of 
the jet and V 
measurements the pitot was relied on as in any intermittent re 


the velccity at a point given by a For mean 
gion, or even outside of the jet, changing velocity can be much 
more easily averaged than with a hot-wire, provided velocity 
fluctuations are not too violent. Rough observation with the 
hot-wire seemed to indicate the same general mean profile of the 
jet 

Some data were taken at other velocities than that used for 
the foregoing analysis and revealed no qualitative differences in 
the conclusions already drawn. The angular spread of the jet, 
however, was found to increase approximately linearly with ve 


locity in the range 35 to 70 ft. per sec., indicating that the in- 
tensity of the turbulence was probably strongly dependent on 


operating velocity. 
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On the Boundary Condition for the Flow Along 
a Perforated Wall 


P. F. Maeder and James B. Carroll 
Brown University, Providence, R.! 
October 13, 1954 


[Ts TRANSONIC WIND TUNNEL TEST SECTIONS, perforated walls 

have been used to eliminate choking. The boundary condi- 
tion at a perforated wall has been investigated, both theoretically 
and experimentally, in reference 1. This note presents another 
theoretical treatment plus some experimental results. In pre- 
vious works, the walls were assumed to be porous and the flow 
through the wall to be of a viscous nature. A linear pressure- 


velocity relationship was thus assumed 
d= —xAp 


In experimental work a perforated wall was used; for this wall 
the hole size is no longer small compared with the wall thickness, 








204 JOURNAL OF THE AERONAU 


TICAEL SCLENCES 




















Fic. 1. Velocities through a two-dimensional slot. 
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Fic. 2. Characteristics of 1-in. slot and 4-in.? hole. 
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and the treatment of the individual holes as orifices was con- 
sidered. This resulted in the quadratic relationship 
d=ea /(2/p)Ap 

This relation is achieved by assuming that the Ap is created by 
the perturbation velocities only. This assumption is incorrect 
if there is an appreciable flow in the x direction, as is found in 
wind tunnels. 

If at first a wall containing a single hole is considered, then, 
assuming incompressible flow, the Bernoulli equation, applied 
between the hole cross section and a point where undisturbed 


flow exists, is 


Ap = (p/2)c? — (p/2)U.,2 (1 
where 
c = velocity at hole 
U., = free-stream velocity 
Ap =?) — Po 


The perturbation velocities in the x, y, z, directions are denoted 
by u’, v’, w’, respectively, and upon substitution into Eq. (1) 
give 
Ap = (p/2) [(U, + u’)? + 0% +w”? — U3] (la) 
Since U,, is much larger than wu’, v’, and w’, second-order terms 
can be neglected. Dividing Eq. (la) by the dynamic pressure, 
(p/2)U.*, gives the dimensionless coefficient. 
Ap/(p/2)U,,2 = 2u'/U, (2) 
Since u’ ~ v’, a linear pressure-velocity relationship is obtained. 
To find the constant of proportionality between u’ and v’, two 
assumptions are made. 
(1) The hole is considered two-dimensional—i.e., a transverse 
slot in the wall. 
(2) The pressure Ap is constant over the slot cross section, and 
only linear terms of the Bernoulli equation are retained. 
Referring now to Fig. 1, at y = 0a slot extends from x = —1 
to x = 1 ina solid wall lying along the x-axis. If the flow along 
the wall is said to consist of two parts—namely, a basic parallel 
flow of velocity U.,, which is perturbed by a flow possessing a 
complex conjugate perturbation velocity c*(z), then the flow 
velocity far away from the slot will be u = U, + a’ and 
aty = 0, 
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=| <2< i, ai’ = Rec* 
i ge —v’' = Imc* = 0 
yo, c* = ii’ 


From an examination of the physical problem, it is seen that two 
conditions must be satisfied: 

(1) The complex potential must have suction points or poles 
sts = 1. 

(2) At x = —1 the Joukowski condition holds, since the edge 
of the slot can be regarded as the trailing edge of a wing; thus 
the complex velocity must possess a zero or stagnation point 
atx = —1. 

The only complex function that satisfies the above two condi- 
tions, transforms the physical plane into one Riemann sheet in 
the hodograph plane, and remains finite at infinity is from refer- 


ence 2. 


ct = Ale + 1/Vz + 1)(s-—1)] = AV (2+1)/2— 1) (3) 


If the vectors z + 1 and zs — 1 are considered in polar form and if 
y— o andc* = a’ = A, then, 
ci=U, 4+ a'’V(24+ 1)/(s — 1) (4) 


and at y = O and |x| < 1 
/ , / 
v= —w' V(x + :1)/(1 — x) 


The average velocity i’ across the hole is obtained from integra- 


e 4 le +1 P ious - 
v= -dx = — ii 5 
2J-1 Vi-«x 2 


‘ 


tion 


Thus from Eq. (2), if u’ = a 
Ap/q = —2u'/U, = 40’/r-U, (6) 


Fig. 2 shows some experimental results of tests made with a 4-sq. 
in. circular hole and a 1- by 4-in. rectangular slot in the upper 
wall. 

If the wall is now considered as containing an infinite number 
of transverse slots, uniformly spaced along the x-axis, as shown 
in Fig. 3, the method described in reference 2 (pages 77 to 93) can 
be applied to obtain the velocities due to the slots. 

The flow must have suction points at z = nl, + h, and the 
stagnation points must be located in such a way that the velocities 
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Fic. 3. Wall containing infinite number of slots. 




















Fic. 4, Wall characteristics of several types of perforated sheets. 
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at the trailing edges of the solid portions remain finite. For 
this condition to be realized, terms like z — (nl, — /;) must ap- 
pear in the numerator. If the sine is used to replace infinite prod- 
ucts, the conjugate velocity c* will be 

sin (2/lz) (z + hh) ae 


1 
N sin (m/le) (2 — hh 


ila 
i |sin (x/ls) (2 + },) 7 mh, gw Oh . 
¥ = \ sin r/lo(z — 1) a le ee lz wi 
and 
vy = = —asin r(h,/l2) (9) 
u’ = —acos r(1,/Is) (9a) 
v’/u’ = tan ra/2 (10) 
where 
2),/le = o = open area/total area 


Extending the analysis to the compressible case by applying the 
Prandtl-Glauert rule to the linearized subsonic flow, one obtains 


v’ To / : 
— = 6 tan ‘ B= V1— MM 11) 
u 2 


Reference 3 gives the beundary condition required fer no block- 
age interference for a two-dimensional doublet in a transonic 
tunnel with two perforated walls as 

o’/a’ = —R = —B cot(r/2VW/3) = 0.781 8 (12) 
From Eqs. (11) and (12), ¢ can be evaluated independently of 
Mach Number. 

A comparison between the experimental and the theoretical 
work shows that a linear pressure velocity relationship does exist 
at a perforated wall. The constant of proportionality between 
u’ and v’ depends on B in such a way as to cancel the 8 obtained 
in the boundary condition for interference free walls, within the 
range of linearized subsonic theory. Results of tests made with 
different types of walls are shown in Fig. 4. For v’/U,, < 10 
per cent (small perturbations), the theory shows good agreement 
with the experimental data, and it shows that the boundary con- 


dition is indeed linear. 
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The Static Pressure Variation in Compressible 
Free Jets* 


Walter R. Warren, Jr. 

Research Assistant, Department of Aeronautical Engineering, 
Princeton University, Princeton, N.J. 

October 15, 1954 


M“™s STATIC PRESSURE PROFILES taken during investiga- 
tions!~* of the mixing of compressible jets with quiescent 
air were found to give insight into the physical structure of a 
spreading jet and information as to the validity of the widely 
accepted assumption that the static pressure of a jet does not vary. 

Results were obtained for initial Mach Numbers of 0.97 and 
2.60. The nozzles used to obtain these conditions are referred 


* This study was conducted in connection with research being performed 
under Contract PR-N6ori, Task Order III, for Project SQUID 
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to as the 7 = 1.0 and \f = 2.6 nozzle, respectively, and both 
have the same external shape and exit diameter of 2.55in. The 
M = 2.6 nozzle was designed to give shock-free parallel flow. 


The weak shocks that were found to exist in the nozzle did not 
affect the overall jet and mixing characteristics. The ten-probe 
static pressure rake used was calibrated in order to test its sensi- 
tivity to direction, and it was found that a rotation of 5 deg 
to stream direction at \J = 0.97 did not affect appreciably the 
pressure reading. Total pressure profiles were taken with a 22- 
probe rake of standard design. For all tests, both jets were oper- 
ated at the isoenergetic condition 

Figs. 1 and 2 are representative of the static and total pressure 
data obtained from the J = 1.0 nozzle experiments, which were 
carried out at test positions that ranged from 0 to 30 nozzle exit 
diameters. Each profile was determined from six separate tests 
and indicates the repeatability of data and symmetry of jet. 
The calculated velocity profiles are shown in Fig. 3. The radial 
location of the negative static pressure peak and the velocity at 
which this peak occurs are presented in Fig. 4 as functions of 
downstream distance. The total and static pressure profiles 
at 15 nozzle diameters for the 1J = 2.6 nozzle are presented in 
Fig. 5. These profiles are shown to illustrate the similarity to 
the 17 = 1.0 nozzle tests, which suggests that the phenomena 
occurring are the same in both cases. 

The turbulence effect on the static pressure readings when 
considered qualitatively tend to accentuate the static pressure 
profile peaks but do not appreciably alter the overall shape. 
This can be deduced from Figs. 1 and 8. It is observed that 
the negative peak pressure occurs always at about the position 
of the maximum radial velocity gradient. Since this is the posi- 
tion of maximum momentum transfer (and maximum shear 
stress) at any particular downstream station, the turbulent ve- 
locity components will be larger there than at other radial points 
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Fic. 1. Static pressure profiles a, b, and c. (14, = 0.97, 
pi/pa = 1.19, T3:1/Ta = 1.00) where p = density, 7 = tempera- 
ture, and subscripts S, 1, and a refer to stagnation conditions, 
nozzle exit conditions, and atmospheric conditions, respectively. 
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Fic. 2. 


of that station. This would imply that the largest correction to 
an observed static pressure profile (which would be a negative 
correction since the effect of turbulence on a static pressure probe 
is to increase the reading from the actual mean pressure) would 
occur at the negative peak. The correction on the axis should 
be small because the velocity gradient is zero there, and at the 
outer edge of the jet the correction obviously goes to zero as the 
velocity goes to zero. 

It seems apparent from these tests that it is not completely 
accurate to assume the static pressure remains constant through- 
out a high Mach Number free jet. At the four-diameter station 
for the 17 = 1.0 nozzle (Fig. la), the measured pressure variation 
from centerline to negative peak is more than 4 per cent of at- 
mospheric pressure. The actual variation in mean static pres- 
sure is undoubtedly larger in view of the previous statement 
concerning turbulent effects. At 15 diameters, this maximum 
variation drops to less than 1.5 per cent of atmospheric pressure 
The maximum static pressure change measured at 15 diameters 
for the 7 = 2.6 nozzle is about 10 per cent of the atmospheric 
value, which indicates that the assumption of constant pressure 
across a free jet becomes less valid with increasing jet Mach 
Number. 

It is interesting to notice that this last conclusion is predicted 
in the development of Prandtl boundary-layer equations, which 
also govern jet flows. In the Prandtl analysis, the radial pres- 
sure gradient is shown to be of the order of magnitude of 17°6 
(where 6 in this case can be considered to be the width of the 
mixing region nondimensionalized by the axial distance). For 
low-speed jets, this value will, in general, be much less than 1, and 
the static pressure may be assumed constant without appreciable 
error. However, the Mach Number effect should be considered 
as one goes to high-velocity jets. 

Study of the various profiles makes it possible to draw inter- 
esting conclusions concerning the physical structure of free jets. 
Consider the fluid that is ejected from the nozzle. This airflow 
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must spread outward from the axis; therefore, the pressure must 
decrease as one goes away from the axis, since a flow from one 
region to another must be accompanied by a pressure drop in the 
This drop exists in the data until the negative 
The air entrained by 


flow direction. 
static pressure peak is reached (Fig. 1). 
the jet starts from zero velocity and atmospheric pressure and 
must experience a pressure drop as its velocity is increased—i.e., 
the effect is noticed as one comes in from large radii toward the 
et axis. This drop is also apparent in the measured profiles 
Fig. 1 

It might, then, seem that the jet can be divided into two re- 
gions: the central region is composed almost entirely of the 
fluid ejected from the nozzle, and the annular region includes only 
the entrained fluid. However, this simple picture is contra- 
dicted experimentally by the consideration that the mass flow at 
zero diameter (that mass injected into the system by the nozzle) 
cannot be fully accounted for at downstream stations by the 
mass flow inside of the position of the static pressure peak 
These mass flows are obtained readily by integrating experi- 
mental density and velocity profiles. Also, one could assume 
that the outer flow was isentropic and calculate the pressure drop 
needed to accelerate the flow to the measured velocities. It is 
found that the experimental deviation from the calculated drop 
is much too large to be explained by turbulence, indicating that 
more than the induction effect is influencing the outer flow. 
Therefore, the static pressure profile that is observed is a com- 
promise of two opposing effects that are equal at the negative 
static pressure peak. The radius at which the peak occurs can 
be considered, not as a boundary between two fluid streams, but 
as a boundary between the regions of dominance of two oppos- 
ing effects: the spreading outward of the jet stream and the 
induction of surrounding fluid. This boundary is well defined 
for some distance downstream of the nozzle. 


Further downstream (x/r; > 22) the sharp distinction between 
the influence regions, as manifested by the clarity of the static 
pressure peak in the data, disappears. It is believed that the 
jet then becomes a highly turbulent core surrounded by a non- 
turbulent annular region as has been suggested by Corrsin and 
Kistler.4 Such a structure is indicated by the roughness of the 
central data in Fig. le. As would be expected, the outer portion 
of the jet still exhibits the induction effect as its strongest 
influencing factor. The spreading effect, although present 
throughout the jet, is no longer dominant in any region. 


Consider the velocity profiles shown in Fig. 3. It is seen that, 
until the 12'/.-diameter station, the negative static pressure peak 
occurs at the inflection point of the profile. The change in pro- 
file shape at this point indicates a change in the mechanism 
governing the flow. Furthermore, the concave down shape in 
the central portion indicates a spreading mechanism and the con- 
cave up shape in the annular region indicates an induction 
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Fic. 4. Correlations of static pressure peak position: (14, = 
0.97, pi/pa = 1.19, T5i1/Ta = 1.00). 
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Fic.5. Static pressure and total pressure profiles at 15 diameters: 
(M, = 2.60, Pi/ Pa = 2.30, Ts le = 1.00), 


effect. This substantiates the conclusion drawn previously by 
investigation of the static pressure profiles only. 

Further conclusions concerning characteristics of the jet can be 
reached by considering Fig. 4. It is noticed that the radius at 
which the negative static pressure peak occurs increases slowly 
with axial distance. It is believed that, if the induction effect 
were the stronger, the peaks would tend to go inward as one goes 
downstream, and, conversely, if the spreading effect were the 
stronger, the peaks would spread further outward than is ob 
The spreading affect in the 17 = 2.6 nozzle tests would 
be larger than in the M = 1.0 nozzle tests, since the momentum 
injected into the system is considerably larger. It is found that 
= 2.6 nozzle occur near the 


served. 


the static pressure peaks for the 17 
value of 7/r; = 1.3 to about 20 nozzle diameters downstream 
This is an appreciable increase over the J = 1.0 nozzle value 
of about 1.05. It is also found that the value of u/u; at which 
these peaks occur is about 0.65 for both series of tests. On the 
basis of the previous conclusion that the inflection point in the 
velocity profile occurs at the position of the static pressure peak, 
it is now reasonable to conclude that the velocity profiles at least 
for the upstream portion of free jets will become fuller (have in 
creasing radii of curvature) near the axis as the jet Mach Number 


increases. 
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Erratum—‘‘Wing-Body Interference at 
Subsonic and Supersonic Speeds—Survey 
and New Developments’’! 


H. S. Bluston and John B. Large 
79 Onslow Road, Richmond, Surrey, England 
November 27, 1954 


= (157), PAGE 310, should read: 
dLz pul b R 


dx x?/R? 1 ~ ‘x2 + hb? Vx? + R? 
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Gaussian Outputs from Non-Gaussian Inputs 


Leo J. Tick 
Research Division, College of Engineering, New York University 


October 18, 1954 


I A RECENT PAPER by Mazelsky,! in the JOURNAL, a proof was 
seemingly given for the Gaussian output of a linear device, 
for a non-Gaussian input. The theorem is incorrect. 

In order to prove his result, the author invokes a central limit 
theorem (there are many). Every central limit theorem in ex- 
istence requires some form cf independence in the sequence of 
random variables that are summed, if not for every set of vari- 
ables, at least for variables that are more than a certain distance 
apart.2 If we examine the assumptions made in the mathematics 
(rather than in the word statement), we see that the author as- 
sumes that the spectrum is flat out to a frequency that is allowed 
to grow indefinitely large. This is an assumption of uncorrelated- 
ness, not independence. Therefore the central limit theorem 
cannot be invoked. 

This distinction is not a fine point. There is a result in random 
process theory? which states that any stationary random process* 
can be represented by the infinite sum of a sequence of uncor- 
related random variables. 

If results such as Mazelsky’s were true, it would imply that 
most processes were Gaussian. Uncorrelatedness implies a sort 
of “independence” in the linear structure of the process. In 
processes with strong nonlinear aspects, there can be uncorre- 
latedness with strong dependence. Consider the sequence‘ 
}X,{ of random variables, where each x,, is distributed uniformly 
This is an 


on the interval (—1, 1), and x,4; = 1 — 2\x,. 
uncorrelated process, but one for which there is an exact (non- 
linear) predictor. In the atmosphere there are strong nonlinear 
components. 

In addition to this basic error, there is a certain looseness in the 
statements made in the paper which do not agree with the mathe- 
matical operations performed. 

Assumption (a) on page 150 states: ‘‘The input power spec- 
trum (two-dimensional), G,(f), is constant with respect to the 
system transfer function H(f)”. I do not understand what the 
“with respect to’’ part means. Mazelsky uses the assumption 
that the power spectrum is constant. From Sketch (a) Eq. (40) 
does not follow. G,(f) must be constant to zero. 

Finally, in summing up, the author states his result as follows: 
“For the special case of a non-Gaussian forcing function that has 
a wide bandwidth compared with the bandwidth of the system 
transfer function, it is shown that for most cases the output re- 


sponse will be Gaussian.” 


* With continuous spectrum 


In the proof, wideness alone of the bandwidth of the forcing 
function is not used, but wideness plus flatness is used. Atmos- 
pheric turbulence spectra may be wide, but they are not flat. As 
for the applicability of the result for ‘‘most cases,”’ there is no 
evidence presented for this. 

These remarks refer only to the second portion of Mazelsky’s 
paper. The first part on higher order moments contains correct 


and worth-while results. 
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Calculation of Span Loading for Arbitrary Plan 
Forms 


John DeYoung 

Aeronautical Research Scientist, Ames Aeronautical Laboratory 
NACA, Moffett Freld Calif: ~ \ 

October 22, 1954 


_ WEISSINGER METHOD (for example, see NACA Rep. 921) 
is limited to use on wings with a quarter chord line of con- 
stant sweep. An attempt to follow the same procedure for a 
quarter chord line having a variation in sweep across the span 
cannot result in a general solution. However, a different treat- 
ment of the integral equation expressing the necessary total 
downwash variation can lead toa simple solution. It is assumed, 
as in the Weissinger method, that a bound vortex at the quarter 
chord line duplicates the wing loading and that the total down- 
wash induced at the three-quarter chord line equals the section 
angle of attack. The integral equation is derived for a quarter 
chord line of varying sweep. 
The downwash due to the load vortex is (see Fig. 1) 


Thds dy : dx 
dw, = = (x — ¥) — (y — 9) 
4 ar 4nr3 "" dy 


/ 
where r = VV (x — £)? + (y — §)? 
trailing sheet is 


The downwash due to the 
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dv dv ( x — # 
dw. = — 1 - 
V Tr 


2x(y — 7) 4al(y — #) 


in which the second term is finite for all J values. The total down- 
wash at the point (+, y) is the sum of the integration across the 
wing span of dw, and dw. In the past, a typical procedure for 
simplifying this integral equation was to specify * and integrate 
«, by parts in terms of I'’(¥). However, if the second term of 
w» is integrated by parts in terms of I'(¥) and combined with zw, 


the equation becomes very simple and dx¥/dy disappears. The 


result is 
b/2 ee b/2 o — 
1 ras) .. . & x— ry) 
y= - d§ + 1— —— Oo 
2rd —6/29 ~— 9 4nJ —b/2 r (y — 9)? 
This equation gives the downwash due to arbitrary loading, I, 
distributed along an arbitrary curve + = f(%) (%, ¥ pertains to 


quarter chord line and x, y to three-quarter chord line). 
The numerical procedure for determining I is similar to that 


of Multhopp and Weissinger.* In trigonometric coordinates, 
5/(b/2) = cos ¢, 
y/(b/2) = cos gy, ¥/(b/2) = He), x/(b/2) = Her) + 
((cy/2)/(b/2)], and (1'/bV) = G(¢) is assumed given by G(¢,) 
m m 
 ® (2hny/sin gy)G,, where hy, = [sin g,/(m + 1)] > 
n=1 wi = 1 


sin “ign SIN wigy = 


COS Na Sin ¢g, Sin gy sin (m + ley 


a Ou ~ ¢¥ 
20m + 1) (cos gu — COS gy) 
SIN gy ra 
St Oe = Em OF Sy = on =F 
where ¢, = [mx/(m + 1)], gy = |uw/(.M + 13], 2 and yu are in- 


tegers, m (an odd integer) is the number of span positions at 
which loading is determined, and .V/ (an odd integer) is the num 
ber of terms used in the numerical integration of the loading at 
Fu 

At the spanwise control point, 7 = cos ¢y, ¢» = [vr/(m + 1)], 
the induced downwash is equal to (w,/V) = a,. The integral 
equation for arbitrary loading can be represented by 


m 


a = > A va Gus y= 1,2,...™ 
n=1 

The simultaneous solution of the m equations gives the value of 
loading, Gr, at m points. The 4,, coefficients are equal to the 
sum of those obtained from the well-known solution of the first 
integral in the w equation and those obtained from a numerical 
integration of the second integral. The coefficients of the first 
integral are given by 


m 


9 
A yy = . > My sin Mig sin Mi¢gy = 
(m + 1) sin @ , | 
—sin ¢ b= (— 3)? v? 
. _n# Fv 
(COS ¢g, — COS ¢,)* m +1 
m+ 1 
: »n=rv 


2 sin ¢y 


Let the multiplicative factor of T in the integrand of the w equa- 
tion be designated N,, (in trigonometric coordinates ) 


l E(¢v) — eu) + bc,/d) 


(cos ¢g»y — cos ¢,y)? r 


With the G(¢,) series, the numerical integration gives the coeffi 


cients of G,, 


VV 
i 

Noyhnpde = a = 
7 J0 M+1,y,=1 


Then the total coefficient is 

* DeYoung, John, and Harper, Charles W., Theoretical Symmetric Span 
Loading at Subsonic Speeds for Wings Having Arbitrary Plan Form, NACA 
Rep. 921, 1948 
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, oe = Aovn 
M+1, l 
The analysis is simplified if the loading is symmetric or anti- 
symmetric. For symmetric loading (G, = Gm +1 —- » 


(m + 1)/2 


a, = 5 tide. y= 12, .(m + 1)/2 


, 1 
where 
. V } 
Ayn = Aoyn t+ Ny, (M + 1)/2 An, (M + 1)/2 + 
M+1 
(M — 1)/2 
> (Nop t+ No, +1 0) hing 
yu l 
~ m+ 1 
lorn, M = 
2 
Adve = » nw 
m+ i 
~~ - Ay m+ 1 ~~ 5 
m+ 1 
hip 2 = - 
m+ 1 
hipt hy 1 — n, we e 
For antisymmetric loading (G, = —Gm 4+ 1 n 
m — 1)/2 
a,= >) PG v=1,2,...(m—1)/2 
n l 
where 
1 VU 1)/2 
Pon = Poon + : } (Now Ny, M+1—p) hn 
M ws 1 , = 1 
Pov» = Ap» —_ Avy, m+il1 Ns Nip = Ray —h m+l1l—n, pu 


In performing the computations, V,, can be evaluated from 
trigonometric tables by use of the relation 


y 1 — sin y ? 1 
Nog ” ; ait N = o 
(cOS gy — COS ¢y)* 2(c,/b)e2 
where 
Egy) — Hep) + le, b) 
tan yy = 
| COS Gy — COS Py 


where £ is a function of the shape of the wing quarter chord line 
(e.g., gy) = €(¢,) = 0 for a straight wing, &¢,) = | cos ul tan 
A, &(¢») = |cos g,| tan A for a swept wing), and c, is the wing 
chord at cos ¢>». 

Once G,, is determined, the lift coefficient, induced drag, and 
spanwise center of pressure for symmetric loading and rolling- 
moment coefficient for antisymmetric loading are determined from 


the following numerical integrations: 


rA ’ — 
C, = Gim + 1)/2 + 2 - G, sin ¢ ) 
m-+ 1 os 1 
. m + 1)/2 
rA Gim + 1)/2 ; 
Cp; = > ¥ 40(m + 1)/2,nG T 
m+ i 2 cat 
m 1)/2 m+ 1)/2 
p ® G, y & ag G sin Fy 
y= il n l 


$4 | 
nep = Gm +72 - 
(m + 1)C; s l,odd \4 — #* 
m — 1)/2 ” (—1)™ + 1)/2 
. tT « FF. tm 
n l “ 1, odd = Mir 
= m— 1)/2 
Cc, = ¥% G, sin 2¢ 
2(m + 1 - 1 


The equations are general for arbitrary odd values of m and 
M. The m gives the number of terms of a loading series and has 
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primary control of the accuracy of the loading distribution. 
However, the accuracy of the a,, coefficients is controlled pri- 
marily by M. It is necessary to find optimum values of m and 
M that give the least computational effort for a sufficiently accu- 
rate prediction. In the present method, it has been found (based 
on a few computations) that m + 1 should be either one-half or 
one-quarter as large as WM + 1. With this ratio, the number of 
simultaneous equations are kept at a minimum. If digital com- 
puters are available and the number of simultaneous equations 
(fixed by the choice of mz) is of little importance, then letting m 
equal M considerably simplifies the computations of the a,, co- 
€.8., Ayn = Aovna + [Non sin ¢n/2(m + 1)]—which 
must be done by hand. 

For swept wings with straight one-quarter chord line a m = 
7, M = 31 solution results in loading that is essentially the same 
as that given by the method of NACA Rep. 921. The compu- 
tations in the present method are easier (since N,, is simple); 
however, the time required is of the same order (since J is 


efficients 


larger ). 

To illustrate the effect of various m’s and M’s for a particular 
variable swept wing, the span loading has been computed by 
m — M procedures of 7-31, 11-11, and 11-47, and presented in 
Fig. 2. The results show that for m = 7, difficulty would arise 
in properly defining the span loading, but this difficulty would be 
largely eliminated by making m = 11 (at the expense of more 
simultaneous equations to solve). It can also be seen that some 
increase in accuracy is obtained when m = 11 by increasing / 
from 11 to 47 (at the expense of computational effort in computing 
Nyy) 


Curved Shocks in Pseudostationary Flows 


Ryuma Kawamura and Hakuro Oguchi 

Cornell University, Ithaca, N.Y. (on leave from the University of 
Tokyo), and Institute of Science and Technology, University of 
Tokyo, 856 Komaba Meguro-ku, Tokyo, Japan, Respectively. 


November 1, 1954 


groan A. H. Taub! made a theoretical study on curved 
shocks in pseudo-stationary flows and derived the relations 
between the nth derivative of the curvature of a curved shock 
and that of the pseudo-stream line behind it under the assump- 
tion that the flow is uniform in front of the shock. We have 
independently treated the same problem in a somewhat different 
way, and naturally arrive at essentially the same conclusion as 
that of Taub. However, in applying the theory to shock-wedge 
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interaction phenomena, we have found some interesting results, 
which seem to be useful clues in obtaining the complete solution 
of the phenomena. 

A two-dimensional nonstationary motion of a gas in the x, 
y plane is called pseudostationary when all flow parameters are 
functions only of the two independent variables, £ and 7 defined as 


f= x/t and n = y/t 


where ¢ is the time. The flow produced by interaction of a plan 
shock with a wedge becomes pseudostationary in case the effects 
of viscosity and heat conductivity are negligible. 

In a mathematical treatment of a pseudostationary flow, it is 
convenient to use the pseudovelocity (U’, V) in place of the par- 
ticle velocity (1, v), which is defined as 

U=u-é and V=v-7n 
A line given in the &, 7 plane by the relation 
dn/dé = V/U 


is called a pseudostream line which corresponds to a particle path 
in the x, y plane. Furthermore, a line element, do, in the &, » 


plane is defined as 
(da)? = (dé)? + (dn)? 


Therefore, it has the dimension not of length but length divided 
by time. This has to be noticed in interpreting the definitions of 
curvature and pressure gradient in a pseudostationary flow, which 
appear in this paper. 

We assume a pseudostationary flow where a shock moves 
outward into a uniform region. At a certain point on the shock, 
let the curvatures of the shock and the pseudostream line behind 
it be K and Ks, respectively, and the pressure gradient along the 
pseudostream line behind the shock be (p2);. Here, we will pre 
sent the relations of these quantities without showing any de- 
tails of the calculations, because our method of obtaining them is 
essentially the same as Taub’s in the sense that they are derived 
from the equations of motion, the equation of continuity and the 


shock relations. They are written as 


B( M2? — 1) cos yo ; 
A sin Yo = » K 
p2Ws 


W.(1 — My? cos? yo) Ks (1) 


[A peW22 cos yz — BW, sin wiJK = 


Wo(1 — Ms? cos? Yo) (pe)s 
where 
A = —W2 + Wi cos yi cos Yo + 
Wy sn ce th = sin yy, sin ye 
(y + 1)\;? cos? yy \ (2) 
P= + piM;? cos x1 sin yr 


and where W is the pseudovelocity, 7 the pseudo-Mach Num- 
ber or W divided by the local velocity of sound, a, p the pressure, 
p the density, y the ratio of the specific heats and y the angle 
measured from the normal of the shock to the pseudostream line 
at the point considered (Fig. 1). The subscripts 1 and 2 are used 
to denote the states in front of and behind the shock, respec- 
tively. 

By virtue of the shock conditions, such quantities as W2, p2, 
and J, in the region behind the shock can be given in terms of 
the quantities in front of the shock together with y as follows: 


sin y 
W. = W,—., 
sin ye 
(y + 1)p1.-W? cos? yy 
ioe 1 9 3 


(y — 1)M,2 cos? yy; + 2 
M, = 


} 1/2 


ne. + 1) .,2 sin yy cos Y, cosec Yo a 
[2y1,2 cos? % — (y — 1))"? [Cy — 1) M2 cos? Yi + 2 
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It must be noted here that these relations are obtained under 
the assumption of the uniformity of the flow in front of the shock. 
Eq. (1) can be reduced to Taub’s zero-th order relation and, as 
he pointed out in his paper, it coincides exactly with the results 
obtained in a stationary case. 

Now, let us apply the above results to some special cases of 
shock-wedge interaction phenomena. 

(a) Curvature of Reflected and Mach Shocks at the Triple Point 
in Mach Reflection: Because of the existence of a common pseu- 
dostream line behind the triple point, called the slipstream, and 
also because of the uniformity of the flow in front of the reflected 
shock, Eq. (1) can be applied to this case and written down in the 
following forms: 


QO Ky = OKs and DiKy = Dp»); 


for Mach shock and 


O3Kr = GiKs and D;Kr = Di po)s 

for reflected shock, where Ky, and Ap are the curvatures of the 
Mach shock and the reflected shock at the triple point, respec- 
tively, and the C,S and D,S are finite valued coefficients in gen- 
eral. The solution satisfying the above equations simultaneously 


1s 
Ky = Kz = Ks = (p:), = 0 (4) 


Hence, it can be concluded that the shocks have zero curvature 
at the triple point. This is correct insofar as the flow variables 
remain regular to, at least, their first derivatives behind the 
shock. Once a singularity appears at the triple point, the present 
theory loses its validity. 

(b) The Case of a Weak Shock: As the strength of the shock 
becomes very weak, .\/2 cos ¥2 approaches unity. Hence, in the 
limiting case the right sides of Eqs. (1) become zero, while the 
coefficients in the left sides do not vanish except for the case of a 
normal shock where ¥; = 7. Accordingly, in a general case, both 
K/Ks and K/(p»), tend to zero for a vanishingly weak shock 
This means that, for a finite value of Ks or (p2)., the curvature of 
such a weak shock becomes infinitesimal 

Applying this result to the triple point in the trivial solution 
of a Mach reflection, where the reflected shock is a sound dis- 
turbance and the Mach shock is the continuation of the incident 
one at the triple point, we can conclude that the curvature of the 
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reflected wave must be zero while that of the Mach shock can 
take a finite value at the triple point. In this case, the Mach 
shock can be connected to the incident one with a finite jump of 
curvature. The pseudovelocity behind the triple point of a 
trivial solution is always supersonic, and hence, the singularity 
mentioned in (a) will never appear in this case 

(c) The Case of a Normal Shock: When the shock is normal to 


the pseudovelocity, 


Then, in Eq. (2), 
A= W, — W, and B=0 
Hence, from Eqs. (1), we have 
Ks/K = 0 (5) 
and 


(p2)s = —[p2WAW, — W2)/1 — M,22)|K, 
or by the use of the shock conditions, 


2a,;? peWs 1 Vv? , 
(po). = : : R (6) 
y+1W,1- MM? 

where a is the local velocity of sound 

These results can be applied to the point of intersection of the 
Mach shock with the wedge surface. By virtue of the first rela- 
tion, the curvature of the Mach shock can be finite at this point 
and the pressure gradient on the wali is calculable provided the 
curvature of the Mach shock is obtained by experiment 

(d) The Case of a Sound Disturbance: When a plane shock 
passes over a wedge of vanishingly small vertex angle, the reflected 
shock becomes almost a sound disturbance, the front of which is a 
circle with its center floating with the same velocity as that of 
flow behind the incident shock. The wave front is everywhere 
normal to the pseudovelocity in front of it and accordingly, the 
results in (c) can be applied. Moreover, in this case, it is appar- 
ent that 

K= ay 1 = q.~! 


Hence, from Eq. (5), Ks = 0 

In Eq. (6), on the other hand, both the numerator and the 
denominator approach zero simultaneously when the shock 
strength becomes very weak. On carrying out a limiting process, 


we have the following result 


(po), = — - pid} (7) 
7 +s 


It is well known that the linearized theory? gives us an infinite 
value of (ps2), behind such a sound disturbance. Eq. (7) gives its 
exact value 

The authors wish to express their thanks to Professors W. R. 
Sears and Y. H. Kuo and other staff members of the Graduate 
School of Aeronautical Engineering, Cornell University, for their 
stimulating discussions and also to Professors C. C. Lin of M.I.T. 
and J. D. Cole at C.I.T. for their helpful advice. 
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An Additional Note on Rolling Pull-Out 
Calculations 


Ralph W. Stone, Jr. 
Aeronautical Research Scientist, NACA, Langley Field, Va. 
November 22, 1954 


T° REFERENCE 1 an analysis of a rolling pull-out calculation 
presented in reference 2 has been made by what might be 
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termed an inverse procedure, the resulting motion having been - 


used to establish the yawing moment which must have existed 
to produce the motion. In this inverse analysis, however, the 
author of reference 1 has failed to note that the motions in refer- 
ence 2 as defined by standard NACA symbols were calculated 
relative to the principal axis such that the values ¢, 6, and y 
and ¢, 6, and y presented in reference 2 are the angular velocities 
about the principal axes, and their integrals with time. The 
author of reference 1, however, erroneously has assumed that the 
values of ¢, 6, and y and ¢, 6, and yp presented in reference 2 are 
the Euler angles and their rates of change with time, whereas 
actually the following relations exist in reference 2 

¢=p=¢,— ¥, sin 8, } 

7] qg = 0, cos ¢, + y, cos 8, sin ¢, > (1) 
y =r= Y, cos 6, cos @, — 6, sin ¢,/ 


where ¢,, 9,, and y, are Euler’s angles. 

Thus, the yawing moment (Fig. 1, reference 1) evaluated by 
the use of Eq. (2) of reference 1, by using values of ¢, 6, ¥, etc., 
from Fig. 4 of reference 2 is incorrect for the reasons stated 
above, and the yawing moment expressed by 


N =I/,y — (I, — 1,)¢6 (2) 


{[Eq. (3) of reference 1] is correct and actually is the moment 
which produced the motion in Fig. 4 of reference 2. The ap- 
proximation of the value of the yawing moment expressed by 

N=I,y (3) 
[Eq. (5) of reference 1] to the actual yawing moment implies 
only that the nonlinear term (J, — /,,)66 was not relatively large 
for the case presented (Fig. 4, reference 2). 

If it is desirable to determine what the Euler angles and the 
rates of change of Euler angles are in such a motion as shown in 
Fig. 4 of reference 2, Eqs. (1) above must be solved at any in- 
stant of time for the desired values. The rates of change of 
Euler angles may be obtained at any instant of time from a 
rearrangement of Eqs. (1) as follows 


6, = cos $,q — sin ¢,r 
o, = tan 8, (sin ¢,q + cos gr) + P 
sin.d@eq + cos or 


cos 6, 


and the Euler angles 0,, ¢,, and ¥, may be obtained by a step-by- 
step integration of these equations. If values of the Euler angles 
and their rates of change are available, then their use in Eq. (2) 
of reference 1 would correctly re-evaluate the yawing moment 
acting and, in fact for the case in point, would have agreed with 
the results given by Eq. (2) above (Eq. (3) of reference 1) and 
shown in the figure of reference 1. 

It should be noted that the equations of motion of a rigid 
body can very readily be expressed in terms of the Euler angles 
such that the solution of the equations result in motions relative 
to the Eulerian coordinates. Eq. (2) of reference 1 would be 
one of the 6 equations of motion so expressed. It would appear, 
however, that such equations would be more cumbersome to 
handle than those expressed in terms of angular velocities about 
the principal axis as were used in reference 2. The implication 
of reference 1 that certain omissions and linearizations were made 
in reference 2 is correct only to the extent that certain degrees of 
freedom and higher order stability derivatives (as, for example, 
B derivatives) were neglected. The three moment equations 
expressed about the principal axes and in terms of angular ve- 
locities about those axes, used in reference 2, are complete and 
correct except for the omission of higher order stability deriva- 
tives just mentioned. No linearizing approximations exist in 


them. 
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Summary of Existing Cone Temperature 
Recovery Factor Data Obtained in 
Continuous-Flow Wind Tunnels Plus 
Additional Data 


Leslie M. Mack 

Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, Calif. 

November 4, 1954 


ges TEMPERATURE recovery factor is the ratio of the actual 
temperature rise from the edge of the boundary layer to the 
adjacent surface, for the case of a surface with no heat transfer, 
to the corresponding adiabatic temperature rise. Since a state 
with no heat transfer seldom exists in many practical cases, the 
main importance of knowing the recovery factor is that the re- 
covery temperature is the zero point of the temperature differ- 
ence which is multiplied by the heat-transfer coefficient to de- 
termine the heat transfer rate. Thus knowledge of the recovery 
factor is basic to the whole heat-transfer problem 

The existing published temperature recovery factor data for 
cones obtained in continuous-flow supersonic wind tunnels are 
summarized in Fig. 1. These data were obtained at four dif- 
ferent laboratories, in ten wind tunnels, and with five cone 
models. The Mach Number range is 1.21 to 6.0, and the Rey- 
nolds Number range is 0.1 & 10° to over 9 X 108!~* In addi- 
tion to these previously published results, recovery factor data 
obtained at the Jet Propulsion Laboratory are presented. These 
data are taken from experiments performed in the two JPL wind 
tunnels using a 13° included angle lucite cone, and two 5° in- 
cluded angle Fiberglas cones. The lucite cone results were ob- 
tained in the course of a recovery factor investigation which is 
reported in fuJl in a JPL report. The Fiberglas cone results are 
taken from a current wind-tunnel transition investigation con- 
ducted by Laufer.’ 

The laminar recovery factors are seen to be in the range 0.825- 
0.865. The result predicted by the laminar boundary-layer 
theory (e.g., reference 8) is Pr /*, where Pr is the Prandtl Number. 
The value of Pr’’*, based on the recovery temperature for a stag- 
nation temperature of 100°F., is near the middle of the experi- 
mental range. The turbulent recovery factors are in the nar- 
rower range from 0.875 to 0.890. The solid points in Fig. 1, 
which give turbulent recovery factors for Reynolds Numbers 
less than 2.5 X 10°, were obtained with artificially tripped 
boundary layers. The value of Pr'/3, which is the simplest re- 
sult given by an analysis of the turbulent boundary layer,’ based 
on the recovery temperature, is higher than the experimental 
values of the turbulent recovery factor. 

The maximum measured value of the recovery factor is 0.91, 
and tt occurs as a peak in the transition region. It was obtained 
at ./ = 4.50 in the JPL 20-in. wind tunnel. The magnitude of 
this transition-region peak varies with Mach Number,® 7 and it 
appears to be different in different wind tunnnels. For instance, 
in the Ames 10- by 14-in. tunnel, the peak value at WM = 4.48 is 
only 0.883.! The wide difference in the transition Reynolds 
Numbers on models in different wind tunnels! is quite apparent 
from Fig. 1. 

The extensive use of cones for the measurement of the re- 
covery factor is due in large part to their simplicity as test in- 
struments compared with another common constant surface- 
pressure shape, the flat plate. In addition, when the working 
surface of the flat plate is aligned with the flow direction, it does 
not provide a boundary layer with zero heat transfer everywhere 
because of the nose-heating effect." 1! For these reasons only 
a limited amount of flat plate data are available for comparison 
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with the cone data." '2 The turbulent recovery factor of refer- VACA Supersonic Wind Tunnels, NACA TN No. 3020, October 
ence 12 agrees with the cone value, but the laminar recovery 1953. 
é : ~ 2 Evve . C., Tucker, M., < sess, W. C atistical Study 0 
factor is almost equal to the turbulent value. The turbulent : Byverd, 5. € Fuster 5 aud Bangees, ™ anise ~ — . * hey . 
Z 7 Transition-Point Fluctuations in Supersonic Flow, NACA TN No. 3100 
recovery factors of reference 10 are also in agreement with the March, 1954 
3 Wimbrow, W. R., Experimental Investigation of Temperature Recovery 


cone results. The distribution of the recovery factor along 
the plate in the laminar region is suggestive of the nose-heating 
effect, which is not the case in reference 12, and when corrected 
for this, a value of 0.85 is obtained in agreement with the cone 
laminar recovery factor. 

The remaining cone and flat plate recovery factor data were 
obtained in intermittent wind tunnels.'*~ The cone results of 
Eber,'*® which are all for laminar flow, are in agreement with 
continuous-flow tunnel results with a scatter about the same as 
in Fig. 1. The flat plate results of Fallis't do not include any 
laminar results, and the turbulent recovery factors are much 
higher than those of Fig. 1, and are in agreement with the cone 
cylinder results of Eber.'* However, an investigation by Gruene- 
wald'® in the NOL wind tunnels with a series of cone-cylinder 
models has shown that the turbulent recovery factor measured 
in the intermittent tunnel is always higher than that measured 
with the same model in the continuous-flow tunnel. Since 
the wind tunnel used by Fallis has a shorter running time than 
the NOL intermittent tunnel, it is not surprising that his results 
agree with Eber rather than with the continuous tunnel results. 
The well-known cone results of Stine and Scherrer were ob- 
tained in the Ames 1- by 3-ft. No. 2 tunnel which has a running 
time of 5-18 minutes, depending on Mach Number, compared 
Eber’s 40-60 seconds, and Fallis’ 22-24 The 
thin-walled cone later used in reference 


with seconds. 
model used was the 10 
1, and the results are in agreement with Fig. 1. 

In summary, the laminar recovery factor, as measured on cones 
in continuous-flow wind tunnels from 14 = 1.2 to 6.0 and Rey- 
nolds Numbers from 0.1 * 10® to over 9 X 108, is 0.845 + 0.020, 
and the turbulent recovery factor is 0.884 + 0.006. 
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Normal Force and Center of Pressure on Right 
Circular Cylinders 


J. Leith Potter, William D. Murphree, and Norman M. Shapiro 
Research Division, Ordnance Missile Laboratories, Redstone 
Arsenal, Huntsville, Ala. 


November 2, 1954 


— NOTE CONCERNS right circular cylinders with major 
axes at small angles of attack in either subsonic or 
supersonic streams. Much of that which is written has applica- 
tion to more general flat-faced bodies of revolution. Although 
the quantity of data for specific conditions is small, the material 
presented herein is thought to enable reasonably accurate esti- 
mates of the static longitudinal stability of right circular cylinders 
having fineness ratios in the usual range of practical interest. 


The right circular cylinder represents a type of aerodynamic 
body which, because of its high drag, is not frequently used. 


’ 


However, it is sometimes desirable to find a “‘finless” and ‘‘spin- 
less’? body of revolution possessing a strong restoring moment, 
drag being a secondary consideration. If it be further required 
that this configuration have a simple and compact geometry, 
there is good reason to consider the right circular cylinder. Vari- 
ous structurally simple modifications which reduce the drag of 
these shapes are well known, but it should be realized that there 
are usually accompanying decreases in stability. Where such 
compromises are feasible, some very satisfactory designs may be 
achieved. 

The writers have completed an investigation which included 
in its objective the collection of data to use for predicting the 
static stability of flat-faced bodies. Both published and unpub- 
lished information available from other sources has been com- 
bined with the results of our own experimental work. Aside from 
the possible usefulness of the data, certain unusual-and inter- 
esting fluid mechanical aspects of the flow over these bodies are 
expected to attract the attention of aerodynamicists. In the 
absence of any theoretical method for predicting the behavior 
of flat-faced profiles, it is necessary to rely on experimental data, 
and the remaining discussion is an attempt to present the more 
essential material in a suitably compressed form. 
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Fic. 1. Local normal force distribution on a right circular 
cylinder. (Experimental points omitted for clarity.) 
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Fic. 2. Center of pressure on right circular cylinders as a fune- 
tion of Mach Number. 


A typical example of the local normal force distribution on right 
circular cylinders is shown by Fig. 1 where the following notation 
applies: 


Cr = (local normal force per unit length) /(qga D) 

D = diameter of cylinder (caliber) 

L = total axial length 

Vo. = free stream Mach Number 

qa = free stream dynamic pressure 

X any axial distance measured from front face of cylinder, positive 


toward base 


a angle between major axis and free stream 


Similar data are available from other measurements by _ the 
writers using a cylinder having an L/D = 2, also in low-speed 
flow. Pressure distributions on cylinders in supersonic flows 
= 2.86.! Data 
for the same cylinder at /,, = 1.88 and only one angle of attack 


have been obtained for an L/D = 5.76 and M 


were measured by the writers in order to verify a phenomenon 
to be discussed later. Rather than present a large number of 
curves, these other pressure distribution data will only be de- 
scribed qualitatively by comparison with Fig. 1. A summary of 
normal force and center of pressure data for right circular cylin- 
ders is given on Figs. 2 and 3. On these figures, Yep is the dis- 
tance from the front face of the cylinder to the center of pressure, 
and Cyq represents the initial slope of the normal force coefficient 
curve. Body cross-section area, 7D?/4, is the reference area in all 
cases. 

Load distributions measured ii incompressible flow on the 
cylinder with L/D = 2 were quite similar to those on the for- 
ward or first two diameters of the longer cylinder which had 
an L/D = 5.15. By inspection of Fig. 1, it is obvious that the 
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Fic. 3. Initial slope of normal force coefficient curve for right 
circular cylinders as a function of Mach Number. 
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short cylinder should possess remarkable longitudinal stability 
and this is verified by Fig. 2 where the integrated results of the 
low-speed tests are presented near the left side of the figure 

The negative loads which were found near the front of both 
cylinders tested in low-speed flow appear to be caused by un 
symmetrical separation of the flow over the nose. Experiments 
with tufts located in the plane of symmetry showed a much more 
extensive separation burble extending back from the nose on 
the leeward side of the body than that on the windward side. 
The fluid, being able to flow around the rim of the face adjacent 
to the windward side of the cylinder in more nearly potential 
flow fashion, produces lower local pressures than those on the 
Based 
on the somewhat incomplete supersonic pressure measurements 
at Mach Numbers of 1.88 and 2.86, it seems that the negative 
local forces at the nose are first smaller, more concentrated, and 


opposite area of the cylinder and negative loads result. 


later disappear as Mach Number increases 

The peaks in positive load which appear roughly one diameter 
from the face in Fig. 1 are to be found at supersonic speeds as 
well, although they are located slightly nearer the face. Unfor- 
tunately, there are no pressure surveys available for these bodies 
in the high subsonic speed range, but the data from force meas 
urements’ shown on Figs. 2 and 3 for Mach Numbers between 
(.5 and 0.9 indicate increasing positive normal force and rear- 
ward movement of center of pressure until the critical Mach 
Number is reached 

Loading on the rearward portion of these bodies is greatly 
dependent on Mach Number. Referring to Fig. 3, it is seen that 
the longer body has less lift than the short one in the low super- 
sonic speed range. In seeking the explanation for this, it was 
noted that reference 2 contains the suggestion that this may be 
due to negative loads on the rear of longer cylinders at these 
Mach Numbers. To verify this, a brief investigation of the 
pressure distribution on a cylinder of L/D = 5.76 was conducted 
with ./,, = 1.88 and four-degrees angle of attack. Substantial 
negative local normal forces were found on the aft 65 per cent of 
the cylinder. However, the opposite was true of measurements 
at \J,, = 2.86! where relatively large positive loads exist on the 
rearward portion of the same body. 

A further point of interest is the diminishing normal force 
curve slope, Cyq@, With increasing angle of attack at subsonic 
speeds. Again, the opposite is true of the data taken at 7, = 
2.86. 

More detailed presentation and discussion of the experimental 
data alluded to herein will be included in a forthcoming report 


from this Laboratory 
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Similar Solutions of Turbulent Boundary 
Layers 


J. C. Rotta 

Aerodynamics Department, The Glenn L. Martin Company, 
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November 14, 1954 


—psetienee IN THIS JOURNAL, Clauser! reported on some very 
interesting results in turbulent boundary-layer research 


FORUM 215 


The purpose of this note is to point out that Clauser’s experi- 
ments are in agreement with the theory of references 2 and 3, 
which had already been published in 1950. The results of these 
papers are the basis for the methods for numerical boundary 
layer calculations which are described in the subsequently pub 
lished paper of reference 4. 

In the earlier papers, references 2 and 3, it is shown that, for 
certain pressure distributions, the equations of turbulent bound 
ary layers have so-called “similar solutions’’ which are character 
ized by the fact that the velocity profiles preserve their shape 
along the wall. Theoretically, these solutions exist only if the 
local skin friction coefficient is constant. These velocity profiles 
depend on two parameters, namely, the dimensionless pressure 
gradient, [x(dp/dx)]/q, and the local skin friction coefficient 
Because the skin friction coefficient changes only slightly as a 
consequence of Reynolds Number variation, it was supposed 
that the similar solutions exist approximately in real boundary- 
layer flow, and that the influence of skin friction coefficient on 
the velocity profiles is unimportant if an appropriate plot is 
chosen. A systematic evaluation of available measurements of 
turbulent boundary layers with arbitrary pressure gradients 
formed the basis for an estimation of the similar solutions. Now 
Clauser has accomplished an experimental verification of these 
solutions which he calls ‘‘equilibrium velocity profiles.’’ Since 
he uses the same procedure which has been suggested in reference 
2 for plotting and characterizing the velocity profiles, a quan- 
titative comparison can easily be made 

The following notations of the two works are identical: 


Reference — ] 2 and 3 
Characterizing parameter G I, 
Universal thickness Asé (5, U,)/(év*) 


The shift of semilogarithmic profiles caused by pressure gra- 
dient, which Clauser denotes by Aus/ux, is expressed in reference 
2 by a change of the magnitude of B. The relations between 
I,, (6, U1) /(6v*), and change of B can be calculated by Eqs. (6, 15), 
(6,16), and (6,19) (with « = 0.4), reference 2. These are based 
on boundary layers with arbitrary pressure distributions. The 
following table gives a summary of the comparison: 


A/é = (6,U,)/(év*) Ate/ux = B— Bex 

Ge-If, Ref. 1 Ref. 2 Ref. 1 Ref. 2 
6.1 3.6 3.5 0 0 

10.1 6.4 6.7 3.2 3.65 

19.3 12.0 13.7 13.6 13.3 


The agreement is reasonably good. This statement is important 
in view of the difference in profile shapes, shown in Fig. 24, refer- 
ence 1. Unfortunately, this figure is not made in the universal 
plot. A direct comparison of the outer portion is, therefore, not 
possible. 

Results of Clauser’s measurements show that equilibrium 
boundary layers without separation may exist up to a nondimen- 
sional pressure gradient (x/q) (dp/dx) = 0.58, whereas the esti- 
mation of reference 2 indicates that similar solutions without 
separation are possible up to (x/q) (dp/dx) = 0.52 for a Reynolds 
Number U6*/vy = 10°. Thus, with regard to this point also, the 
discrepancy is not serious 
Clauser—that two different shape-preserving boundary layers 
already resulted 


Moreover, the case detailed by 


may exist for the same pressure distribution 
from this calculation, as is seen from Fig. 10, reference 2. The 
author did not, however, attempt to discuss this phenomenon be- 
cause he was not sure at that time whether it had a real physical 
meaning or whether it was merely caused by the use of the extrap- 
olating formulas. 

The pressure distributions used for evaluation in reference 2 
are rather different from those of the ‘“‘similar solutions."” There 
fore, the fact that, in spite of this difference, the results agree 
satisfactorily with Clauser’s measurements, is at least some 
support for the hypothesis which treats the outer portion of the 
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velocity profiles as a one-parametric family. Despite this, furthér 
investigation should be undertaken to show whether or not serious 
deviations of the velocity profiles from the one parametric family 
are demonstrable. The next step then would be to explore these 


deviations. 
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Temperature Field in a Laminar 
Incompressible Round Jet and in Convergent 
and Divergent Round Nozzle 


Dragutin Stojanovic 
Asistent Masinskog fakulteta, Beograd, Yugoslavia 
November 9, 1954 


— EXACT SOLUTION of the Navier-Stokes equation for all 
plane flows in which the stream lines of the viscous flow are 
the same as in potential flow has been given by G. Hamel.! H. 
Squire? has used the same transformation in axisymmetrical flow 
and by interpretation in spherical coordinates has given the ve- 
locity distribution in a laminar incompressible round jet. At 
the same time he has given the solution of the energy equation but 
neglecting the viscous dissipation term. For a case of plane flow 
through a divergent and convergent channel as a special case of 
flow in logarithmic spirals it is possible to find the exact solution 
of the energy equation® by using the transformation: 


l 
R? F\(0) + Ty (1) 


In the following it will be shown that this transformation is valu- 
able also in an axisymmetrical two-dimensional Hamel flow. 


The energy equation for this case in spherical coordinates 
(R, 6, @) with the identical vanishing of the velocity component 
V4 and the differential quotient for coordinate ¢ is: 


; (v - Vo OT = 1 O (oF 
wade ihe! el 2S eee ar) | 
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where there are the strain tensor components: 


OVe Vr Vo 
oR €oo = R + R cotg @ 


Ver 1 OVe 1 OVe ' re) Vo 
+ 7R0 = a 
R R 060 R 00 OR\R 
The transformation (1) and the transformation given by G, 
Hamel and H. Squire for velocity: 
V F.(6 V F;(6 
= — F,(6) = — F,(6) 
R R 6 R 
change the partial differential Eq. (2) to the ordinary differentia] 
equation: 


pep( F\’F; — 2F\ Fo) = k(F,” + F,’ cot g @ — 2F,) + 


2u E + 2F.F;’ + F;’? — 2F.F; cotg 6 + 


2F,'F; + F;2(2 + cotg? 0 | (3) 


By this it is shown that the transformation (1) makes possible 
the exact solution of the energy equation for axisymmetrical 
Hamel flow. For a divergent and convergent round nozzle the 
velocity component Vg vanishes identically and with that the 
function F; and Eq. (3) gives: 

k(F,” + F,’ cotg 0) + 2Fi(pcpF, — k) = —p(6F. + Fo”) (4 
i.e., the exact solution of the energy equation for axisymmetrical 
analogous to the Jeffery-Hamel flow. 
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